AD-752  558 


V/STOL  AIRCRAFT  AERODYNAMIC  PREDICTION 
METHODS  INVESTIGATION.  VOLUME  II.  APPLI 
CATION  OF  PREDICTION  METHODS 


Peter  T.  Wooler,  et  al 
Northrop  Corporation 


V _ _ 

I 

r - -  ' 

Prepared  for: 

Air  Force  Flight  Dynamics  Laboratory 
January  1972 


DISTRIBUTED  BY: 


National  Technical  Information  Service 
U.  S.  DEPARTMENT  OF  COMMERCE 

5285  Port  Royal  Road,  Springfield  Va.  22151 


AD  752558 


AFFDL-TR-72-2# 
Volume  II 


V/STOL  AIRCRAFT  AERODYNAMIC 
PREDICTION  METHODS  INVESTIGATION 


Volume  II.  Application  of  Prediction  Methods 


P.T.  Wooler 
H.C.  Kao 

M.F.  Schwendemann 
H.R.  Wasson 
H.  Ziegler 


Northrop  Corporation 
Aircraft  Division 


TECHNICAL  REPORT  AFFDL-TR-72-26.  Volume  II 
January  1972 


Approved  for  public  release,  distribution  unlimited. 


Repr  *  ri  r  ;  f, , 

NATIONAL  TECHNICAL 
INFORMATION  SERVICE 


D  D  C 


i  K 

U; 

n 


n  1372 


—  B  ^ 


AIR  FORCE  FLIGHT  DYNAMICS  LABORATORY 
AIR  FORCE  SYSTEMS  COMMAND 
WRIGHT  PATTERSON  AIR  FORCE  BASE,  OHIO  45433 


NOTICE 


When  Government  drawings,  specifications,  or  other  data  are  used  for  any 
purpose  other  than  in  connection  with  a  definitely  related  Government  procurement 
operation,  the  United  States  Government  thereby  incurs  no  responsibility  nor  any 
obligation  whatsoever;  and  the  fact  that  the  Government  may  have  formulated, 
furnished,  or  in  any  way  supplied  the  said  drawings,  specifications,  or  other  data, 
is  not  to  be  regarded  fay  implication  or  otherwise  as  in  any  manner  licensing  the 
holder  or  any  other  person  or  corporation,  or  conveying  any  rights  or  permission  to 
manufacture,  use,  or  sell  any  patented  invention  that  may  in  any  way  be  related 
thereto. 


lussm  a r 


PTIS  Wilia  SkUn 

d:C  Bat;  fectlu  a' 

«na;.  □ 

Jl&IU  iw/.  IlliJ 


IT  .  .  . . .  ... 

tmmm/mmrn 


1  list  AWL  mi/w  KMl  1 

4 

Copies  of  this  report  should  not  be  returned  unless  return  is  required  by  security  con¬ 
siderations,  contractual  obligations,  or  notice  on  a  specific  document. 


A1  R  FORCE/56780/30  November  1972  —  300 


/ 

»  -. 


Unclassified 


••emit*  Classification 

|  DOCUMENT  CONTROL  DATA  •  R  &  D 

1  (Security  c  las  al  tic  alien  of  title,  Mr  et  Mulrcci  M  annetatlon  aw«i  M  an'M  when  the  overall  report  la  claaaltled) 

ORIOINATINO  ACTIVITY 

a rn*  REPORT  SECURITY  CLASSlPlC  ITtON 

I  Northrop  Corporation 

Unclassified 

1  Aircraft  Division 
|  Hawthorne,  California 

li  CROUP 

[*  REPORT  TITLE 


V/STOL  AIRCRAFT  AERODYNAMIC  PREDICTION  METHODS  INVESTIGATION  - 
Volume  n 

«■  DCSCRIPTIVC  notcs  (Trpt  mi  rsparf  mnd  JncJusIys  dsfss)  ~ 

_  Final  Report  -  1  May  1969  to  31  January  1972 _ 

••  authorisi  ;rintn—f,  mid*,  InltUI.  Iialamj 


Peter  T.  Wooler  Howard  R.  Wasson 

Hsiao  C .  Kao  Henry  Ziegler 

Myles  F.  Schwendemann  _  _ 


••  REPORT  OATE  ’  '  " 

January  1972 

•S.  CONTRACT  OR  GRANT  NO.  ""  " 

am.  ORIGINATOR'S  REPORT  NUMIERIH 

F33615-69-C-1602 

*.  PROJCCT  NO. 

NOR  72-9 ,  Volume  n 

698  BT 

*.  T ask  No.,  698  BT  01 

•6.  OTHER  REPORT  NOtSI  (Any  ether  numbere  that  may  be  aealpned 
tflta  report) 

A 

AFFDL-TR-72-26,  Vo’ume  n 

10  DISTRIBUTION  STATtMINT 

Approved  for  public  release;  distribution  unlimited. 

II.  SUPPL CMSNT ARY  NOTCS 

•  ^  Ai«f»ACT  “ 

U.  SPONSORING  MILITARY  *  CTIV1TY 

Air  Force  Flight  Dynamics  Laboratory 
Wright-Patterson  Air  Force  Base 

Ohio  45433 

Analytical  engineering  methods  are  developed  for  use  in  predicting  the  static  and  dynamic 
stability  and  control  derivatives  and  force  and  moment  coefficients  of  lift-jet,  lift-far , 
and  vectored  thrust  V/STOL  aircraft  in  the  hover  and  transition  flight  regimes.  The 
methods  take  into  account  the  strong  power  effects,  large  variations  ir  angle  of  attack 
and  sideslip,  and  changes  in  aircraft  geometry  that  are  associated  with  high  disk  loaded 
V/STOL  aircraft  operating  in  the  aforementioned  flight  regimes.  The  aircraft 
configurations  studied  have  a  conventional  wing,  fuselage  and  empennage.  The  pre¬ 
diction  methods  are  suitable  for  use  by  design  personnel  during  the  preliminary  design 
and  evaluation  of  V/STOL  aircraft  of  the  type  previously  mentioned. 

This  report  consists  of  four  volumes..  The  prediction  methods  are  applied  to  a  number  of 
V/STOL  configurations  in  this  volume..  The  theoretical  development  of  the  prediction 
methods  may  be  founa  in  Volume  I.  Details  of  the  computer  programs  associated  with  the 
prediction  methods  are  given  in  Volume  HI.,  The  lesults  of  a  literature  survey  are  presented 
in  Volume  IV. 

Details  of  illustrations  in 
this  document  may  be  better 
studied  or  microfiche 


DD 


FORM 

I  NOV  • 8 


1473 


Unclassified 


!  :'c. 


Security  Classification 


V/STOL  Aircraft 
XV-5A 

V/&TOL  Stability  and  Control 

Aerodynamic  Characteristics  (Subsonic) 

Forces  (Normal  Force,  Lift,  Side  Force) 

Moments  (Pitching,  Soiling,  Yawing) 

Rotary  Derivatives 

V/STOL  Aerodynamics 

Nonlinear  Aerodynamics 

Wing  Stall 

Separated  Flows 

V/STOL  Transitional  Flight 

Engine  Wake  Effects 

Jet  Exhaust  Fields 

Jet  interference  Effects 

Jet  in  Cross  Flow 

Jet  Path 

Multiple  Jets 

Jet  Interaction 

Power  Effects 

Engine  Inlet  Effects 

Vortices 

Computer  Programs 


«U  S  Government  Printing  Office  1972  -  7E9-483/195 


V/STOL  AIRCRAFT  AERODYNAMIC 
PREDICTION  METHODS  INVESTIGATION 


Volume  II 

Application  of  Prediction  Methods 


P.T.  Wooler 
H.C.  Kao 

M.F.  Schwendemann 
H.R.  Wasson 
H.  Ziegler 


Approved  for  public  release,  distribution  unlimited. 


/  ^ 


*  t* 


-warn 


FOREWORD 

This  report  summarizes  the  work  accomplished  by  the  Aircraft  Division  of 
Northrop  Corporation,  Hawthorne,  California  for  the  Air  Force  Flight  Dynamics 
Laboratory,  AFSC,  Wright  Patterson  Air  Force  Base,  Ohio,  and  USAF  Contract 
No.  F33615-69-C-i602  (Project  698  BT).  This  document  constitutes  the  Final  Report 
under  the  contract. 

This  work  was  accomplished  during  the  period  1  May  1969  to  31  January  1972, 
and  this  report  was  released  by  the  authors  in  January  1972.  The  Air  Force  Project 
Engineers  were  Mr.  Robert  Nicholson  ar.d  Mr.  Henry  W.  Woolard  of  the  Control 
Criteria  Branch,  Flight  Control  Division,  AFFDL.  Their  assistance  in  monitoring 
the  work  and  providing  data  is  greatly  appreciated. 

The  authors  gratefully  acknowledge  the  assistance  and  '  .  operation  of  NASA 
Langley  Research  Center  personnel  during  the  wind  tunnel  Try.--  iyi  testing  in  the  NASA 
Langley  V/STOL  tunnel. 

Special  recognition  is  due  Mr.  Richard  J.  Margasor,  of  NASA  Langley  Research 
Center  who,  besides  being  ac  lively  involved  in  the  test;  ng  at  Langley,  has  made  valu¬ 
able  contributions  to  other  areas  of  the  investigation. 

Various  people  at  Northrop's  Aircraft  Division  contributed  to  the  investigation, 
particularly  the  following  persons  in  the  areas  des  gnated: 

Lynn  B.  Fricke  Developed  empirical  methods  for  the  wing.  Was  Test 

Engineer  for  the  wind  tunnel  testing  of  the  component 
model. 

Hsiao  C.  Kao  Developed  the  transformation  method  for  estimating 

power  effects  on  wings  and  fuselages.  Developed  the 
empirical  method  for  (he  body. 

Myles  F.  Schwendemann  Developed  the  method  for  estimating  engine  inlet  effects. 

Provided  prediction  method  wind  tunnel  testing  interlace 
for  the  configuration  model.  Pardcipa'od  in  the  hover 
analysis. 

.  i 


Martin  F.  Silady 

Howard  R.  Wasson 


Peter  T.  Wooler 

Henry  Ziegler 


Assembled  the  V/STOL  bibliography  and  was  responsible 
for  the  literature  survey. 

Developed  the  method  for  mapping  general  sections. 
Developed  the  nonlinear  body  prediction  method.  Assisted 
with  the  development  of  the  perturbation  method. 

Directed  the  technical  effort  and  developed  the  nonlinear 
wing  prediction  method. 

Developed  the  jet  flow  field  prediction  method.  Per¬ 
formed  the  analysis  of  wing  power  effects  employing 
lifting  surface  theory. 


Contributions  have  also  been  made  by  U.  A.  G.  Brynjestad  to  this  study  in  a  num¬ 
ber  of  areas,  particularly  the  literature  search  and  perturbation  method  developmer.  ; 
by  M.  S.  Cahn  in  the  development  of  the  method  for  mapping  general  sections  and  per¬ 
turbation  method  development;  by  members  of  the  Northrop  Aerosciences  Laboratories 
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ABSTRACT 


Analytical  engineering  methods  are  developed  for  use  in  predicting  the  static 
anti  dynamic  stability  and  control  darivatives  and  force  and  moment  coefficients  of 
lift-jet,  lift-fan,  and  vectored  thrust  V/STOL  aircraft  in  the  hover  and  transition 
flight  regimes.  The  mathodt  take  into  account  the  strong  power  effects,  large 
variations  in  angle  of  attack  and  sideslip,  and  changes  in  aircraft  geometry  that  are 
associated  with  high  disk  loaded  V/STOL  aircraft  operating  in  the  aforementioned 
flight  regimes  The  aircraft  configurations  studied  have  a  conventional  wing, 
fuselage  and  empe.mage.  The  prediction  methods  are  suitable  for  use  by  design 
personnel  during  the  preliminary  design  and  evaluation  of  V/STOL  aircraft  of  the 
type  previously  mentioned. 

This  report  consists  of  four  volumes.  The  prediction  methods  are  applied  to 
a  number  of  V/STOL  configurations  in  this  volume.  The  theoretical  development  of 
the  prediction  methods  may  be  found  in  Volume  I.  Details  of  the  computer  programs 
associated  with  the  prediction  methods  are  given  in  Volume  III.  The  results  of  a 
literature  survey  are  presented  in  Volume  IV. 
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SECTION  I 


INTRODUCTION 

This  volume  is  the  second  of  three  volumes  treating  V/STOL  Aerodynamic 
prediction  methods.  This  rolanr*  is  directed  toward  presenting  applications  of  the 
methods  developed  during  the  study  program.  These  applications  represent  the 
complete  prediction  techniques  developed  except  for  empirical  methods  for  treating 
power  effects  on  wings  and  bodies  which  are  presented  in  Volume  I. 

1.  PURPOSE 

The  purpose  of  this  volume  is  to  demonstrate  the  use  of  the  prediction 
methods  in  computing  the  aerodynamic  coefficients  and  derivatives  of  V/STOL  aircraft. 
To  accomplish  this  purpose  example  problems  are  worked  out  using  the  prediction 
techniques.  The  accuracy  of  the  methods  can  be  assessed  by  comparing  the  predicted 
aerodynamic  coefficients  against  test  results.  Limitations  of  the  methods  are  described 
and  the  necessary  assumptions  and  simplifications  which  must  be  made  are  pointed  out. 
It  should  be  possible  for  the  user  to  evaluate  the  described  methods  and  to  assess  their 
limitations  for  use  in  his  problems  by  examining  the  included  sample  problems,. 

2.  SCOPE 

The  methods  which  are  applied  in  this  volume  represent  the  theoretical 
prediction  techniques  which  have  been  developed  during  the  study  contract.  The  one 
theoretical  thechnique  which  is  not  presented  is  the  vortex  tracking  method  which  is 
not  applied  to  c.  given  example  problem  since  the  method  was  not  operable.  Also  not 
included  in  this  volume  are  the  empirical  methods  for  treating  wing  or  body.  These 
methods  are  treated  completely  in  Volume  I. 

In  presenting  samples  of  the  use  of  the  method  the  specific  nature  of  the 
chosen  problem  often  does  not  demonstrate  the  full  capability  of  the  methods  being 
employed.  To  further  demonstrate  the  method  capabilities  it  is  frequently  pointed  out 
where  the  methods  are  more  general  and  where  different  problems  can  be  treated  by 
the  same  methods. 
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3.  TECHNICAL  APPROACH 


The  method  of  demonstrate^  die  application  of  Hie  individuil  methods  which 
was  adopted  for  this  volume  was  to  select  a  sample  problem  and  to  apply  Hie  methods 
to  this  problem  where  applicable.  Hie  details  involved  in  using  each  method  are 
explained  by  presenting  a  complete  treatment  of  foe  sample  problem.  Where  computer 
programs  are  involved  complete  inputs,  outputs  and  comparison  of  results  with  experi¬ 
ments  are  presented.  Where  computer  programs  are  not  involved  the  equations  used 
are  presented  and  comparison  with  experiment  is  presented. 

Sample  problems  have  been  chosen  where  test  data  is  available  for  comparison 
with  theoretical  results.  This  permits  foe  accuracy  of  the  prediction  techniques  to  be 
evaluated  and  permits  a  better  understanding  of  the  difficulties  to  be  encountered  in 
predicting  the  aerodynamics  of  V/STOL  aircraft. 

4.  ORGANIZATION  OF  VOLUME  H 

The  mata  body  of  Volume  H  is  devoted  to  sample  problems  treating 
V/STOL  aircraft  configurations.  The  treatment  is  divided  into  sections  considering 
the  power  induced  effect*,  and  the  power  off  effects. 

Section  n  presents  two  methods  of  treating  power  effects  on  the  wing  of  a 
V/STOL  configuration.  Each  of  these  methods  has  certain  advantages  for  the  user  and 
both  will  be  useful  in  treating  power  induced  effects  on  wings. 

Section  HI  presents  the  method  of  predicting  power  induced  effects  on  the 
fuselage  and  illustrates  the  accuracy  which  can  be  obtained. 

Section  IV  demonstrates  the  use  of  the  prediction  methods  in  obtaining  down- 
wash  and  sidewash  effects  at  the  empennage  location  and  illustrates  how  these  results 
can  be  used  to  estimate  power  effects  on  the  tail  surfaces . 

Section  V  presents  the  method  applicable  to  predicting  the  power  effects  of 
the  inlets.  Comparisons  are  made  with  test  data  to  illustrate  the  accuracy  of  the 
method. 

Sections  VI  and  VII  present  the  methods  for  treating  the  unpowered  effects  of 
nonlinear  body  and  wing  aerodynamics.  Section  VI  presents  the  nonlinear  body  method 
and  shows  comparisons  with  test  data.  Section  VII  presents  similar  comparisons  for 
wing  aerodynamics. 
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Four  appendices  are  added  to  this  volume  to  complete  the  documentation 
of  work  performed  under  the  contract.  The  first  three  appendices  describe  the  wind 
tunnel  test  program  undertaken  during  the  study.  This  description  i3  presented  in 
thi.3  volume  to  permit  the  details  of  the  test  program  and  model  description  to  be 
available  to  the  user.  This  is  desirable  since  much  of  the  test  data  used  for  compari¬ 
son  with  the  p  ediction  method  came  from  the  test  program.  A  method  for  estimating 
normal  force  and  pitching  moment  in  the  lift  jet  wake  region  is  presented  in  Appendix 
IV. 
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SECTION  II 


POWER  EFFECTS  ON  THE  WING 

There  are  two  alternative  procedures  for  computing  the  power  induced  effects  on 
the  wing.  The  fi^st  of  these  combines  the  mapping  method,  the  jet  flow  field 
program  and  the  transformation  method,  hi  this  method  the  mapping  program  is  used 
to  describe  the  wing,  the  jet  program  calculates  velocities  induced  by  the  jet  at  the  wing 
surface  and  the  transformation  method  calculates  pressures,  forces  and  moments 
induced  by  the  jet  on  the  wing. 

In  the  second  procedure  the  jet  program  is  used  in  conjunction  with  the  lifting 
surface  theory  to  predict  power  induced  wing  effects.  The  jet  program  computes  a 
downwash  field  at  the  wing  plane  which  gives  an  effective  camber  distribution  for  the 
wing.  The  lifting  surface  theory  then  utilizes  this  camber  distribution  to  compute  lift 
and  pitching  moment  on  the  \/ing.  This  method  does  not  give  the  pressure  distribution 
about  the  wing  but  is  simpler  and  easier  to  use  than  the  first  method. 

The  following  subsections  will  describe  the  use  of  both  of  these  methods  as  applied 
to  a  given  wing.  The  use  of  the  methods  will  be  described  in  detail  so  that  a  complete 
understanding  of  the  methods  can  be  obtained. 
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1.  SAMPLE  PROBLEM 


To  demonstrate  the  method  of  predicting  the  power  induced  aerodynamic  effects 
on  a  wing,  a  sample  problem  is  given.  This  problem  is  for  a  single  jet  in  the  presence 
of  an  isolated  wing.  The  wing  of  the  sami  problem  is  the  ore  tested  in  the  configura¬ 
tion  wind  tunnel  test  program  discussed  in  Appendix  I. 


Section  coordinates: 
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Jet  location  and  description: 


A  sketch  of  the  coordinate  system  defining  the  jet  location  relative  to  the  wing 
is  shown  below. 


For  the  sample  problem: 


yj 

z. 

J 

Jet  diameter  d. 

J 


=  .  6  inches 
=  0 

=  -6.635  inches  (positive  upward) 
=  2. 25  inches 


The  jet  exhausts  into  the  mainstream  along  the  negative  z-direction,  that 
is  perpendicular  to  the  wing  planform. 


00 

Velocity  rrtio  -jj - =  .  2 

Jo 

Wing  attitude: 

Angle  of  attack  q  =  0  degrees 
Sideslip  angle  0  ~  0  degrees 


2.  APPLICATION  OF  MAPPING  METHOD  TO  THE  WING  ,f 

lii  the  sample  problem  all  the  wing  sections  are  geometrically  similar  so  that  it  is 
only  necessary  to  map  a  single  section  and  then  to  scale  the  coefficients  to  provide  a 
mapping  for  different  spanwise  stations  Figure  1  shows  the  inputs  to  the  mapping 
program  used  to  obtain  the  initial  mapping. 

The  first  card  lists  in  order  the  number  of  coordinates  being  input,  the  number  of 
corners  (and  pseudo-corners)  being  input,  the  number  of  terms  to  be  taken  in  the  ex- 

I  i 

pansion  for  the  potential  and  the  mapping,  and  a  zero  to  indicate  that  the  section  being 
mapped  is  symmetrical  about  the  x-axis. 

Cards  2  through  5  give  the  x-coordinates  selected  as  input  points  starting  at  the 
section  trailing  edge  and  proceeding  arcund  to  tfye  nose.  Cards  6  through  9  are  the 
z-coordinates  of  the  airfoil  at  the  same  points.  Stece  the  section  is  symmetrical  only 
the  upper  half  plane  coordinates  are  input, 

Gard  number  10  specifies  that  the  airfoil  is  to  be  shifted  .  5  units  in  the  negative 
x-direction.  This  is  necessary  since  to  obtain  the  mapping  it  is  necessary  that  the  air¬ 
foil  be  centered  with  respect  to  the  origin  to  some  degree.  It  is  not  necessary  to  center 
the  section  exactly  but  it  should  be  centered  somewhere  near  the  centroid  of  area. 

Card  11  specifies  which  input  points  are  comer  points.  This  card  indicates  that 
the  first  input  point  (the  trailing  edge)  is  a  comer  point.  The  second  number,  the 
zero,  indicates  that  the  second  comer  is  a  pseudo-corner;  that  is,  is  pot  a  true  comer 
but  merely  indicates  a  region  of  large  curvature.  This  second  number  refers  to  the 
nose  of  the  airfoil. 

Card  12  specifies  the  x  and  z  coordinates  of  the  first  comer  and  the  angle  turned 
through  at  the  comer  in  radians. 

Card  13  similarly  describes  the  second  corner  specifying  the  location  of  the  center 
of  the  leading  edge  radius  and  specifying  that  a  comer  of  n  radians  is  turned.  This  angle 
turned  is  only  an  approximation;  it  is  not  necessary  to  specify  the  angle  exactly. 

Card  14  specifies  parameters  needed  to  tell  the  program  how  to  write  out  the 
mapped  section  with  the  comers  included.  The  first  two  numbers  are  the  x-  and  z- 
coordinates  of  the  initial  poiit  to  be  mapped,  in  this  case  the  leading  edge  point.  The 
third  and  fourth  numbers  specify  the  first  and  last  points  to  be  mapped  specified  as 
angular  distances  around  the  mapping  circle.  In  the  the  sample  section  the  mapping 
is  to  start  at  the  nose  (180  degrees)  and  proceed  around  the  lower  surface  until  just 
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FIGURE  1.  MAPPING  PROGRAM  INPUT  DATA  FOR  SAMPLE  PROBLEM  (Wing) 


ahead  of  the  trailing  edge  (355  degrees).  The  last  number  specifies  that  mapped  points 
are  to  be  obtained  at  increments  about  5  degrees  apart.  In  this  example,  it  was  neces¬ 
sary  to  choose  these  parameters  to  avoid  specifying  the  trailing  edge  as  one  of  the  end 
points.  This  is  because  the  trailing  edge  is  a  corner  and  a  corner  point  cannot  be 
specified  as  one  of  the  end  points. 

Card  15  specifies  the  necessary  parameters  needed  to  map  out  points  on  the 
section  with  the  comers  removed.  This  card  soecifies  that  37  points  are  to  be  printed 
out  with  a  spacing  increment  of  5  degrees  about  the  mapping  circle  and  that  the  mapping 
is  to  start  at  9  =  0  degrees  on  the  mapping  circle. 

The  outputs  of  the  mapping  functions  for  the  inputs  of  Figure  1  are  shown  in 
Figure  2.  The  first  page  of  outputs.  Figure  2(a) ,  relates  to  the  computations  made  in 
calculating  the  potential  and  the  point-to-point  correspondence  between  points  on  the 
section  and  points  on  the  mapping  circle. 

The  first  two  qolumns  represent  x-  and  y-coordinates  of  the  input  points  with  the 
x-coordinate  shifted  by  the  incremental  value  of  AX  input  (in  this  case,  -.  5).  The 
third  column  represents  the  distance,  R,  from  the  new  origin  to  the  point  on  the  section. 
The  fourth  column  gives  the  computed  perimeter,  S,  of  the  section  from  the  positive  real 
axis  to  the  point  in  question.  The  fifth  column  gives  the  velocity,  V,  calculated  due  to 
the  unit  vortex  about  the  body.  The  velocities  written  out  for  comer  points  are  mean¬ 
ingless.  ALPHA  is  the  angle  of  the  section  slope  at  each  point  as  calculated  in  the 
program.  OMEGA  specifics  the  angular  distance  in  degrees  around  the  section 
measured  about  the  new  origin.  THETA  is  the  predicted  angular  distance  or  the  points 
around  the  mapping  circle  in  degrees. 

The  second  page  of  printouts,  Figure  2(b),  represents  the  results  obtained  by 
computing  the  derivative  of  the  mapping  function  with  the  comers  contained  explicitly 
and  integrating  the  resulting  expression  numerically.  The  location  of  the  first  point 
was  specified  as  explained  above  (this  point  is  not  printed  out).  The  three  columns 
give  the  x-coordinate,  the  y-coordinate  and  the  angle  around  the  mapping  circle  0 
for  a  series  of  points  as  specified  by  the  parameters  specified  in  the  input  cards.  The 
degree  to  which  these  points  agree  with  the  original  section  represents  the  degree  of 
accuracy  obtained  by  using  the  mapping  method. 

The  third  page  of  outputs,  Figure  2(c),  represents  the  results  of  multiplying  out 
factors  representing  any  comer  singularities  and  integrating  the  resulting  expression 
analytically. 
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FIGURE  2,  MAPPING  PROGRAM  OUTPUT  DATA  FOR  SAMPLE  PROBLEM  (Wing) 


SECTITN  V^PPINC  ev  NU^EMICAl  INTEGRATION 
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(b)  SECTION  MAPPING  BY  NUMERICAL  INTEGRATION 


MAPPING  Of  SECTION  *1Th  CC«NfcfiS  «£MnVtO 
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<d)  MAPPING  OF  SECTION  WITH  COHN  MRS  REMOVED 


First,  the  radius  of  tbe  mapping  circle  is  printed  oat  as  compoted.  Next,  tbe 
coefficients  of  tbe  mapping  function  are  written  oat  -  first  tbe  real  parts  aQ  —  files 
fie  imaginary  parts  bQ,  Tbe  mapping  ftnc&oc  is  written  in  fbe  form 

-  *  .dt+ibt  .  + 

. *  <» 


The  coefficients  are  written  out  in  order  a^,  a2>  a3  *  *  *  *  an  *°r  Partica^ar 
case,  all  the  imaginary  parts  of  tbe  coefficients  are  zero  since  tbe  section  is  symmetrical. 

Page  four.  Figure  2(d) ,  prints  out  the  x-  and  y -coordinates  obtained  with  tbe 
analytically  integrated  mapping  function.  It  is  not  possible  to  obtain  the  location  of  tbe 
body  directly  by  this  method  of  computation;  i.  e. ,  the  constant  term  in  the  above 
mapping  is  missing.  This  prevails  the  location  of  the  section  being  mapped  from 
being  specified.  This  is  readily  remedied  by  plotting  the  original  section  and  the 
mapped  section  and  the  displacement  required  to  obtain  a  good  fit  between  tbe  two 
sections  represents  the  constant  term  of  the  mapping. 

Figure  3  shows  a  comparison  between  the  mapped  output  and  the  original  input 
section,  the  lower  surface  also  being  shown  since  the  section  is  symmetrical  and  the 
mapping  retains  the  proper  y  of  symmetry. 

To  obtain  coefficients  for  the  wing  of  the  sample  problem  it  is  necessaiy  to  change 
the  coefficients  to  account  for  the  size  of  the  actual  wing.  The  coefficients  of  Figure  2 
are  based  on  a  unit  chord  and  the  mapping  for  this  section  can  be  written: 


where 


a' 

=  .60815  x  10  1 

1 

-3 

a  2 

=  .  75775  x  10  ° 

-4 

a'3 

=  .69368  x  10 

-6 

*4 

=  .12786  x  10 

a'_ 

=  -.23231 x  10" 

0 

-6 

a’e 

=  15*1-04  x  10 

14 


Original  Airfoil 
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15 


FIGURE  3.  COMPARISON  OF  ORIGINAL  AND  MAPPED  AIRFOIL 


and 


a'?  =  .  11822  xio"7 

a'g  =  -.46527  x  l(f10 
a'g  =  -.44950  1 10~S 


*.26916 


a*  has  not  been  determined  and  it  le  best  to  leave  this  coefficient  undefined  until  file 
o 

coefficients  have  been  ratioed  up  to  t  rue  size. 

S  is  desired  to  reexpress  Equation  (.)  in  the  form: 

|U . (2) 

where  flie  coefficients  a^,  a 2,  .  .  .  .  ag  reflect  the  true  living  dimensions  at  given 
stations  along  the  wing.  Rewriting  Eqtrion  (1)  in  the  form 

. 


and  equating  cj  to  f  we  obtain 


2  =  C*ca^ 


+ 


c,0a'9 

*  V 


so  that  it  can  be  shown  that 


an-cn*’a'n  . ,9 

i  0 

The  radius  of  the  new  mapping  circle  defined  by  £  =  r  e  is  now  r  =  cr' 

c  c  c 

To  obtain  the  coefficient  a  it  is  sufficient  to  note  that  from  the  first  and  last  numbers 

o 

of  the  last  page  of  Figure  2,  the  mapping  without  the  constant  coefficient  maps  the 
section  about  a  chordwise  point  of  .  4891,  so  that  for  the  origin  of  the  wing  located  at  the 

nose  of  the  root  chord  and  a  leading  edge  sweep  of  9.76  degrees,  the  coefficient 
aQ  can  be  found  as 

aQ  =  ,4891c  +  y  •  tan  9.76  degrees, 

where  c  -  c  -  (c  -  c.)  y/,  , 
r  '  r  t' J  b/, 
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For  use  in.  **»e  transformation  method,  it  is  also  necessary  to  define 
This  can  be  done  by  noting  that 


for  the  wing. 


dr 
_ c 

dy 


r  =  .26916  c 
c 


dr 
_ c 

dv 


. 26916 


dc 

dy 


K  is  now  possible  to  compute  all  of  the  numbers  needed  for  the  transformation  method. 
These  numbers  are  tabulated  in  Table  1. 
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3.  APPLICATION  OF  JET  FLOW  FIELD  THEORY  TO  WING 


a.  Transformation  Method 

Tbe  purpose  of  the  Jet  Flow  Field  theory,  when  used  in  conjunction  with  the 
Transformation  Method,  is  to  predict  jet-induced  velocity  components  at  those  control 
points  on  the  wing  at  which  tbe  Transformation  Method  requires  them  to  evaluate 
power  effects.  This  is  accomplished  by  executing  the  Jet  Flow  Field  computer  program 
to  generate  the  required  data  for  the  Transformation  Method  in  the  form  of  punched 
data  cards.  To  insure  compatibility  with  the  Transformation  Method,  the  control 
points  on  the  wing  are  specified  by  utilizing  the  mapping  coefficients  for  the  wing  cross 
sections  obtained  by  the  procedure  described  in  Section  II.  2.  The  punched  output  is 
generated  in  a  manner  to  provide  a  continuous  block  of  input  data  to  the  Transformation 
Method  computer  program.  Both  of  the  above  features  will  be  described  in  greater 
detail  in  the  discussion  below. 


(1)  Sample  Problem  Computation 

For  the  sample  problem  being  considered,  the  Jet  Flow  Field  program  is  now 
used  to  compute  the  jet-induced  velocities  at  the  eight  spanwise  stations  on  the  wing 
described  in  Section  II.  2.  Figure  4(a)  shows  a  sketch  of  the  wing  and  the  location  of 
the  jet  with  respect  to  the  input/output  coordinate  system.  Figure  4(b)  defines  the 
jet  exhaust  angles  <b  and  V  .  It  should  be  noted  that  the  input/output  coordinate  system 
shown  below  differs  from  the  general  coordinate  system  utilized  in  Section  n,  Volume  I. 


FIGFRE  4.  COORDINATE  SYSTEM  FOR  TYPICAL  WING 
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(a)  Input  for  Sample  Problem 

The  input  cards  required  for  the  sample  problem  are  tabulated  in  Figure  5. 

Card  1  lists  three  control  indices.  The  fiist  one,  MULT  =  1,  indicates  that  a 
single  jet  configuration  is  being  treated.  The  second  one,  IGE$M  =  1,  specifies  that 
control  points  on  wing  cross  sections  will  be  generated,  utilizing  the  mapping  coef¬ 
ficients  obtained  in  Section  n.  2.  The  third  control  index,  1PUNCH  =  1,  generates  the 
punched  output  for  the  Transformation  Method  program. 

Card  2  specifies  angle  of  attack  a  =  0,  and  angle  of  sideslip  /?  =  0. 

Card  3  controls  the  number  of  steps  and  the  step  size  in  the  numerical  integra¬ 
tion  of  the  equations  of  motion  for  the  jet  path.  For  the  sample  problem,  90  steps  with 
a  step  size  of  0.4  jet  exit  diameters  are  chosen. 

Cards  4  and  5  contain  information  about  the  jet.  The  jet  location  in  the  coor¬ 
dinate  system  of  Figure  4  is  X  =  0. 6,  Y  =  0. ,  Z=-6. 63.  The  jet  exhaust  angles  0  and 
defined  in  Figure  4(b)  are  180  and  0  degrees  respectively.  The  jet  exit  diameter, 
do  =  2.25  and  the  velocity  ratio  U^/U^  =  0,  2. 

Card  6  may  be  left  blank  for  computations  r.volving  a  single  jet.  For  a  multiple- 
jet  configuration  it  would  be  used  to  control  the  geometry  of  the  jet  resulting  from  the 
intersection  of  two  other  jets. 

The  remaining  input  cards  provide  data  to  generate  the  control  points  at  which 
jet-induced  velocities  are  to  be  evaluated.  These  control  points  are  generated  by 
utilizing  the  mapping  coefficients  and  mapping  radii  obtained  in  Section  II.  2  for  the 
eight  wing  stations  of  the  sample  problem.  The  numbej  of  control  points  generated 
at  each  spanwise  station  is  governed  by  the  input  variable  NTHT,  which  is  the  number 
of  increments  4  &  into  which  the  mapping  circle  is  to  be  divided  in  the  Transformation 
Method  computer  program. 

Card  7  specifies  that  NTHT,  the  number  of  control  points  at  each  spanwise 
station  or  the  number  of  46  increments  around  the  mapping  circle  is  36,  and  that 
the  number  of  spanwise  stations  NS  =8  .  It  also  defines  the  number  of  terms  used 
in  the  mapping  expansion  of  Section  II.  2,  NC0EF  =  11,  and  through  the  control  index 
IRECT  =  1  indicates  that  the  wing  is  nonrectangular. 

Cards  8  —  12  provide  che  da«,a  from  which  the  wing  cross  section  at  the  first 
spanwise  station  can  be  generated  by  the  computer  program.  Card  8  specifies  the 
location  of  the  station,  Y  =  0, ,  the  mapping  circle  radius  R  =  2.8889  and  the  rate  of 
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FIGURE  5.  JET  FLOW  FIELD  PROGRAM  INPUT  DATA  FOR  SAMPLE  PROBLEM 

(Wing;  Transformation  Method) 


change  of  R  with  Y,  DRDY  =  0.  Cards  9  —  12  list  the  real  and  imaginary  parts  of  the 
coefficients  to  be  used  in  the  mapping  expansion.  Cards  13  —  47  are  similar  data 
blocks  for  the  spanwise  stations  Y  =2.5,  5.04,  7.7925,  10.545,  13.,  16.05  and  20.125. 
The  data  in  cards  8—47  are  taken  from  the  Table  I  in  Section  II.  2. 

Note:  The  rate  of  change  of  the  mapping  circle  radius  with  spanwise  distance, 
DRDY,  is  not  required  for  any  of  the  computations  performed  by  the  Jet  Flow 
Field  program.  It  will,  however,  be  required  by  the  Transformation  Method 
program  and  is  read  as  part  of  the  input  that  it  may  be  punched  out  in  the 
proper  sequence  in  the  data  package  to  be  provided  to  the  Transformation 
Method  program. 

(b)  Output  for  Sample  Problem 

For  the  problem  being  considered  both  printed  and  punched  outputs  are  obtained. 

•  Printed  Output: 

Figure  6(a)  shows  the  first  page  of  printed  output.  The  jet  configuration 
being  treated  is  identified  both  by  appropriate  heading  and  by  other  pertinent 
input  information.  Input  controlling  the  numerical  integration  procedure  is 
also  displayed.  Figure  6(b)  shows  the  location  of  the  computed  jet  centerline. 
The  nondimensionalized  jet  speed  U^/U^  and  the  nondimensionalized  major 
diameter  of  the  ellipse  representing  the  cross  section  of  the  jet,  d/do  are 
also  out.  'Hiese  properties  are  printed  out  at  each  integration  inter¬ 

val.  Output  shown  in  Figure  6(b)  shows  only  the  initial  portion  of  the  printed 
output  generated  for  this  example.  The  centerline  was  computed  to  Z  =  87.63, 
which  represents  integration  of  the  jet  equations  over  the  range  Z  =  90  x  0.4 
x  2.25  =  81.  The  variables  XOOORD  and  DLA  show  a  monotonic  increase 
over  this  region,  while  UJ  =  Uj/Uj0  «.  2  as  the  mean  jet  speed  approaches 
the  freestream  speed  l)®.  Figure  6(c)  shows  the  printout  for  the  jet- 
induced  velocity  components  at  the  first  spanwise  station  specified,  Y  =  0. 

The  coordinates  of  the  36  control  points  at  the  station  are  identified.  The 
induced  velocity  components  U,  V,  W  all  nondimensionalized  by  U  are 
printed  out  for  each  control  point.  Figure  6(d)  shows  the  output  for  the  last 
station  considered  in  the  sample  problem,  Y  =  20. 125.  Similar  printouts 
are  obtained  for  the  other  intermediate  stations  specified  as  part  of  the  input. 
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***  SINGLE  JET 

CONFIGURATION  *** 

XJET  VJET 

ZJET 

“HI 

PS! 

U/UJC 

0.6000  0.0 

-6.6300 

180.0000 

0.0 

0.2m  0 

ANGLE  OF  ATTACK 

0.0 

ANGLF  OF  SIDESLIP  * 

O.C 

NUMBER  OF  STEPS  IN  INTEGRATION  =  90 

integration  interval 

*  C.40  JET  EXIT 

diameters 

FIGURE  6(a).  INPUT  PARAMETERS  FOR  SAMPLE  PROBLEM 


** 

SINGLE  JE1 

CENTFPL INF  ** 

******************************** 

XCOORC 

VCOORO 

2COORD 

UJ 

D1A 

0.6C 

0.0 

-6.63 

1.000 

l.OC 

0.61 

0.0 

-7.53 

0.948 

1.18 

0.64 

0.0 

-8.43 

0.893 

1.45 

0.T1 

0.0 

-9.33 

0.833 

1.90 

0.83 

0.0 

-10.23 

0.760 

2.64 

1.01 

0.0 

-11.13 

0.688 

2.93 

1.27 

0.0 

-12.03 

0.626 

3.23 

1.62 

0.0 

-12.93 

0.573 

3.55 

2.06 

0.0 

-13.83 

0.528 

3.89 

2.61 

0.0 

-14,73 

0.489 

4.23 

3.28 

0.0 

-15.63 

0.456 

4.60 

4.09 

0.0 

-16.53 

0.427 

4.98 

5.06 

0.0 

-17.43 

0.402 

5.37 

6.21 

0.0 

-18.33 

0.381 

5.77 

7.56 

0.0 

-19.23 

0.363 

6.18 

9.15 

0.0 

-20.13 

0.347 

6.60 

11.02 

0.0 

-21.03 

0.333 

7.03 

13.21 

0.0 

-21.93 

0.321 

7.47 

15.78 

0.0 

-22.83 

0.311 

7.92 

18.78 

0.0 

-23.73 

0.301 

8.38 

22.29 

0.0 

-2*t.63 

0.293 

8.86 

26.39 

c.o 

-25.53 

0.285 

9.34 

31.18 

0.0 

-26.43 

0.278 

9.85 

36.78 

0.0 

-27.33 

C.272 

10.37 

43.34 

0.0 

-28.23 

0.266 

10.91 

51.00 

0.0 

-29.13 

0.261 

11.47 

FIGURE  6(b).  JET  CENTERLINE  FOR  SAMPLE  PROBLEM 


•**  INOUCEO  VFIOCITIFS  ON  WING  **• 


i 

r  i 

f  ! 

f  ; 

j 

******************************** 

l  ! 

1 

X 

¥ 

l 

U 

V 

M 

! 

I 

10.717 

0.0 

0.0 

0.1136?E-01 

0.0 

-0.47648F-01 

'• 

I 

10.623 

0.0 

0.020 

0. 1 1 582E-01 

0.0 

-0.47592F-01 

l 

10.346 

0.0 

0.050 

0. 12220F-01 

0.0 

-0.47586F-01 

i 

i 

9.004 

0.0 

0.097 

0.13282E-01 

0.0 

-0.475356-01 

i 

1 

0.319 

0.0 

0.159 

0. 14771E-01 

0.0 

-0.47360E-01 

1 

8.616 

0.0 

0.229 

0. 16691E-01 

0.0 

-0. 469856-01 

1 

7.817 

0.0 

0.302 

0. 19044F-01 

0.0 

-0.46294E-01 

1 

6.948 

0.0 

0.377 

0.21786F-01 

0.0 

-0.45121E-01 

| 

6.043 

0.0 

0.449 

0.24793F-01 

0.0 

-0.43336F-01 

5.133 

0.0 

0.506 

0.27913E-01 

0.0 

-0.40918F-01 

\ 

4.242 

0.0 

0.533 

0.3C975F-01 

0.0 

-O.37918E-01 

3.394 

0.0 

0.529 

0. 33823E-01 

0.0 

-0. 343946-01 

f 

2.578 

0.0 

0.496 

0.36272E-01 

0.0 

-0.304616-01 

1.847 

0.0 

0.442 

0. 381 73F-C1 

0.0 

-0.26360E-01 

I 

1 

1.214 

0.0 

0.373 

0.39493F-01 

0.0 

-0.22435F-01 

1 

0.699 

0.0 

0.292 

0.40340E-01 

0.0 

-0. 190396-01 

1 

| 

0.321 

0.0 

0.20? 

0.40931 E-01 

0.0 

-0.1 646 1F-0 1 

1 

0.092 

0.0 

0.105 

0.41508E-01 

o.o 

-0. 14885E-C 1 

1 

i 

0.015 

0.0 

-0.000 

0.42236F-01 

0.0 

-0.1 44 10F-01 

m 

0.092 

0.0 

-0.105 

0.43145E-01 

0.0 

-O.15127E-01 

fl 

0.321 

0.0 

-0.203 

0.44135E-01 

0.0 

-0.17164E-01 

m 

0.699 

0.0 

-C. 292 

0.4497 lc-01 

0.0 

-0.206196-01 

1 

1.214 

0.0 

-0.373 

0.45273E-01 

0.0 

-0.254216-01 

I 

1.94? 

0.0 

-0.442 

0.44584E-01 

0.0 

-0.312C2E-01 

m 

I 

2.578 

0.0 

-0.496 

0.42576F-01 

0.0 

-0. 37266E-01 

m 

i 

3.384 

0.0 

-0.529 

0. 39280E-01 

0.0 

-0.42753E-01 

F 

4.242 

0.0 

-0.533 

0.350996-01 

0.0 

-C.46970E-01 

K 

1 

5.133 

0.0 

-0.506 

0. 30601E-01 

0.0 

-0.49642F-01 

i 

6.043 

0.0 

-0.449 

0.262766-01 

0.0 

-0.50929E-01 

i 

6.948 

0.0 

-0.377 

0.22435F-01 

0.0 

-C. 512236-01 

i  & 

7.81? 

0.0 

-0.30? 

0.1922OF-01 

0.0 

-0 . 50° 1 RP-0 1 

I 

j 

9.616 

0.0 

-0.229 

0. 16644E-01 

o.o 

-0.50293F-01 

: 

f 

9.319 

0.0 

-0.159 

0. 1 4654F-01 

o.o 

-0.49538F-01 

1 

\ 

9.904 

0.0 

-0.097 

0.131 77F-01 

0.0 

-0.488096-01 

; 

10.346 

0.0 

-o.oso 

0. 12 1 54E-01 

0.0 

-0.492?  3E-01 

) 

i 

10.623 

0.0 

-0.020 

0.11553E-01 

0.0 

-0 . 47839T- 0 1 

1 

******************** 

************ 

FIGURE  6(c).  INDUCED  VELOCITY  COMPONENTS  AT  STATION  Y  =  0. 
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*••  INOUCFO  VELOCITIES  OH  WING  *♦* 


** 


X 

y 

2 

U 

V 

W 

a.e?i 

20.125 

0.0 

0.111 50F-01 

-0.13150E-01 

-0.43097F-02 

8.774 

20.125 

o.oio 

0.1ilo3F-01 

-0. 131 15E-01 

-0.429035-07 

8.635 

20.125 

0.025 

0.11206E-01 

-0.13024E-01 

-0.42270E-02 

8.414 

20.125 

0.049 

0.11273E-01 

-0.12878E-01 

-0.412495-07 

8.122 

20.125 

0.079 

0.1 1357E-01 

-0.12682E-01 

-0. 39885F-0? 

7.770 

20.125 

0.114 

0. 11454E-01 

-0.12443E-01 

-0. 38213E-02 

7.37C 

20.125 

0.151 

0.11558E-01 

-0. 17167C-01 

-0.36773F-02 

6.936 

20.125 

0.189 

0. 1 I661F-01 

-0.H862E-01 

-0. 3*1 76E-0? 

6.483 

20.125 

0.225 

0.1175BF-01 

-0.11538E-01 

-0.31849E-02 

6.029 

20.125 

0.253 

0.11 847F-01 

-0.1121 1 E-01 

-0.29503F-0? 

5.583 

20.125 

0.267 

0.11929E-01 

-0.1C892E-01 

-0.27121E-0? 

5.154 

20.125 

0.264 

0. 12004E-01 

-0.10588E-01 

-0.247436-0? 

4.751 

20.125 

0.248 

0. 1 2070F-01 

-0. 10303F-01 

-0.22434E-02 

4.385 

20.125 

0.271 

0. 12127F-C1 

-0.10046E-01 

-0.20272F-02 

4.069 

20.125 

0.186 

0. I2174E-01 

-0.9P750F-02 

-0.18339E-02 

3.811 

20.125 

0.146 

0. 12214F-01 

-0.96477E-02 

-0.16705E-02 

3.62? 

20.125 

0.101 

0. 12249F-01 

-0.95725E-02 

-0.15429E-02 

3.508 

20.125 

0.052 

0. 127805-01 

-0.94554F-07 

-0.14541E-02 

3.469 

20.125 

-0.000 

0. 1731  IE-01 

-0.944B7F-02 

-0.14057E-02 

3.508 

20.125 

-0.057 

0.12339E-01 

-0.95020E-02 

-0.13995E-02 

3.622 

20.125 

-0.101 

0.12362F-01 

-0.96 1455-02 

-0.14375E-0? 

3.811 

20.125 

-0.146 

O.12377F-01 

-0 . 97  845E-02 

-0. 15195F-02 

4.069 

20.125 

-0.186 

0. 123B3E-01 

•0.100075-01 

-0. 16427F-02 

4.385 

20.125 

-0.221 

0.12374E-01 

-0.1O273F-01 

-0.18025E-92 

4.751 

20.125 

-0.248 

0.12^476-01 

-0. 10571F-01 

-0. 19943E-02 

5.154 

20.125 

-0.264 

0.177995-01 

-0.1OS89F-01 

-0 • 221 76F-02 

5.533 

20.125 

-0.267 

0.12725F-O1 

-0.11 212F-01 

-0 . 24522F-02 

6.029 

20.125 

-0.253 

0.171 26F-01 

-0 . 1 1 53 1F-01 

-0. 27081F-02 

6.433 

20.125 

-0.225 

0. 120046-01 

-0.11 836F-01 

-0. 29737E-0? 

6.936 

20.125 

-0. 189 

0.11 365E-01 

-0.17173F-01 

-0. 32386E-02 

7.370 

20.125 

-0. I5i 

0. 117196-01 

-0. 17385F-01 

-0. 34906E-02 

7.770 

20.125 

-0.114 

0.115745-01 

-0. 12614E-01 

-0.371985-07 

8.122 

20.125 

-0.079 

0.  11440E-01 

-0.12804E-01 

-0.39191E-0? 

8.414 

70.125 

-O.049 

0.11322E-01 

-0 . 12954E-01 

-0.40830E-02 

8.635 

20.125 

-0.025 

0.11731F-01 

-0. 13064E-01 

-0.42056F-02 

8.774 

20.125 

-0.010 

0.111 73C-01 

-0. 1 31 31E-01 

-0.47P13E-02 

******** •***♦***♦♦** ************ 


FIGURE  6(d).  INDUCED  VELOCITY  COMPONENTS  AT  STATION  Y  =  20. 125 


f-*~  Punched  Output  for  Station  Y 


0.0 

2.888900 

0.0 

1 

0.10000E  01 

0.0 

9.S2495E  01 

0.0 

0.70057E  01 

0.0 

l 

1 

0.93694E  00 

0.0 

0.92054E  00 

0.0 

0.18211E-01 

0.0 

1 

2 

-0.35514E  00 

0.0 

0.25274E  01 

0.0 

0.20819E  01 

0.0 

1 

3 

-0.87942E-01 

0.0 

-0.91189E  01 

0.0 

1 

4 

0.11362E— 01 

0. 11582E-01 

0.12220E-01 

0. 13282E- 01 

0.14771E-01 

0.16691E-01 

U 

1 

1 

0.19044E-01 

0.21786E— 01 

0.24798E-01 

0.27913E-01 

0.30975E-01 

0.33823E-01 

U 

1 

2 

0. 36272E-01 

0.38173E-01 

0.39493E-01 

0.40340E-01 

0.4093 IE-01 

0.41508E-01 

U 

1 

3 

0.42236E— 01 

0.43145E-01 

0.44135E-01 

0.44971E-01 

0.45273E-01 

0.44584E-01 

U 

1 

4 

0.42576E— 01 

0.39280E-01 

0.35099E-01 

0.30601E-01 

0.26276E-01 

0.22435E-01 

U 

1 

5 

0. 1922 OE -01 

0.16644E-01 

0.14654E-01 

0. 13177E-01 

0.12154E-01 

0. 11553E-01 

U 

1 

6 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

V 

1 

1 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

V 

1 

2 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

V 

1 

3 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

V 

1 

4 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

V 

1 

5 

0.0 

0.0 

0.0 

0.0 

0.0 

O.C 

V 

1 

6 

-0.47648E— 01— 0.47592 E— 01-0 .47586E-01- 

-0.47535E-01- 

-0.47360E-01- 

-0.46985E-01 

w 

1 

1 

-0.46294E— 01- 

-0.45121E-01 

-0.43336E-01-0.40918E-01 

-0.37918E-01- 

-0.34394E-01 

u 

1 

2 

-0.30461E-01-0.26360E-01-0.22435E-01-0.19039E-01-0.16461E-01-0.148B5E-01  W  1  3 
-0. 14410E— 01-0. 15127E-CI— 0. 17164E-0 1— 0.20619E  01-0.25421E-01-0.31202E-01  W  i  4 
-0.37266E-01-0.42753E-01-0.469/OE-01-0.49642E-01-0.50929E-01-0.51223E-01  MIS 
-0.50918E-01-0.50293E-01-0.49538E-01-0.48809E-01-0.48223E-01-0.47839E-01  M  1  6 


2.500000  2.709399  -0.071767  2 

0.10000E  01  0.0  0.S3S33E  01  0.0  0.61620E  01  0.0  2  1 

0.77289E  00  0.0  0.71218E  00  0.0  0.13214E-01  0.0  2  2 

-0.24166E  00  0.0  0.16130E  01  0.0  0.12461E  01  0.0  2  3 

-0.49365E-01  0.0  -0.48006E  01  0.0  2  4 


0. 12325E-01  0. 12S30E-01  0.13133E-01  0.14131E-01  0.15516E-01  0.17283E-01  U  2  1 

0.19419E-01  0.21872E-01  0.24526E-01  0.27241E-01  0.29896E-01  0.32376E-01  U  2  2 

0.34548E-0L  0.36300E-01  0.37600E-01  0.38516E-01  0.39194E-01  0.3980SE-01  U  2  3 

0.40465E-01  0.4U84E-01  0.41866E-01  0.42329E-01  0.42304E-01  0.41493E-01  U  2  4 

0.39697E-01  0.36945E-01  0.33499E-01  0.29741E-01  0.26035E-01  0.22652E-01  U  2  5 

0.19746E-01  0.17366E-01  0.15494E-01  0.14084E-01  0.13098E-01  0.12514E-01  U  2  6 

-0.95974E-02— 0.9591 1E-02-0.96194E-02-0.96592E-02-0.96937E-02-0.97102E-02  V  2  l 
-0.96874E-02-0.95914E-02— 0.93949E— 02-0. 9096 IE-02-0. 8 7079E-02- 0.82 380E-02  V  2  2 

-0. 76951 E-02-0. 71 089E-02-0. 653 13E-02-0.60244E-02-0. 5647 8E-02-0. 544 80E-02  V  2  3 

-0. 54534E-02-0.56833E-02-0.61559E-02-0.68779E-02-0. 78195E-02-0.88918E-02  V  2  4 

-0. 99505E-02-0. 10833E-01 -0.1 14 12E-0 1-0. I l648E-01-0.il 59 1E-01-0.11345E -01  V  2  5 

-0. 11011E-01-0.10654E-01-0. 10314E-01-0. 10021E-01-0.98005E-02-0.96614E-02  V  2  6 

-0.45398E -01-0.45 330E-01-0.45268E-0 1-0.4513 IE-01-0.6 485 1E-01-0.44363E-01  M  2  1 

-u.t3a»i>c-ui-u.*»^3ojt-ui-U.^0639fc-01-0.38415E-01-0.35756E-0l-0.32725E-01  W  2  2 

-0.29415E-0l-0.260HE-Ol-0.22772E-Oi-°.1997lE-Ol-0.17849E-Ol-0.16573E-01  W  2  3 
-0. 16235E-01-0. 16907E-01-0.18663E-01-0.21537E-01-0.25428E-01-0. 30030E-01  W  2  4 

-0.34843E-01-0.39282E-0l-0.42857E-01-0.45330E-01-0.46750E-01-0.47357E-01  W  2  5 
-0.47420E-01-0.47150E-01-0.46704E-01-0.46225E-01-0.45819E-01-0.45545E-01  M  2  6 


FIGURE  7.  JET  FLOW  FIELD  PROGRAM  PUNCHED  OUTPUT 
FOR  SAMPLE  PROBLEM 
(Wing;  Transformation  Method) 


5.0*0000  2.527200  >0.0717*7  3 

0.100006  01  0.0  0.54592E  01  0.0  0.536136  01  0.0  3  1 

0.6272*6  00  0.0  0.539116  00  0.0  0.932796-  02  0.0  3  S'. 

-0.159166  00  0.0  0.9*0916  00  0.0  0.714C4E  00  0.0  3  i 

-0.263856-01  0.0  -0.2393*6  01  0.0  3  * 


0.138066—01  0.139836-01  0.1*5126-0  1  0.1537*6-01  0.165616-01  0.180396-01  U  3 

0.197*16-01  0.2172*6-01  0.23762E-01  0.257906-01  0.277306-01  0.295186-01  U  3  t 

0.3107*6-01  0.3235*6-01  0.333326-01  0.3*0556-01  0.3*6096-01  0.350936-01  U  3  3 

0. 355736-01  0.3*0*56-01  0.36*666-01  0.366*96-01  0.365686-01  0.359856-01  U  3  * 

0. 36*036-01  0. 330096—01  0.3070*6-01  0.2*07*6-01  0.253*16-01  0.227136-01  U  3  5 

0.203*96-01  0.1*3356-01  0.166976-01  0.15*306-01  0.1*5286-01  0.139876-01  U  3  6 

-0.1693 56 -01-0. 169116-01-0. 11.91 36-01-0. 169056 -01-0. 1686  IE-01-0. 16 164E-01  V  3  I 
-0.165856— 01— 0.1628*6— 01 -0.1583 7E— 01-0.152596-01— 0.1*5826-01-0. 138296-01  V  3  2 

-0.130176-01— 0.121*66-01— 0.113976-01-0.107216-01— 0.102326-01-0.9986*6-02  V  3  3 

-0.1001*6-01— 0.103*36-01-0.109756-01— 0.1190*6— 01—0.131026— 01-0. 1*4506-01  V  3  * 

—3. 15*336— Cl— 0. 11 0606-01-0. 1799*6-0 1-0. 1*5566-01-0.187586-0 1-0. 186906-01  V  3  5 

-0.l**5*6-0l-0. 1*1316-01-0. 177776-0 1-0. 174*86-01-0.1 71886-01-0.170196-01  V  3  6 

—0.390006— 01 -0. 389166-01 -0.3S763E-0 1-0. 3<3*8*6-0l-O.38033E~0 1-0. 373666-01  U  3  1 

-0.36*266 -01-0. 351**E-01 -0.33*9*6-0 1-0. 3 1 5196-01-0.2929 16-0 1-0. 268736-01  W  3  2 

-  0.2* 3**E-01 -0.21 *296-01 -0.19*90 6-0 1-0. 17*956—0 1-0.159956-01-0.150956-01  *  3  3 

—0.1*85*6— 01-0. 153096-01 -0.16*9 36-01-0. i**086-01— 0.209816-01-0.2*9376-01  M  3  * 

-0.273166 -01-0. 305096-01-0. 333276-0 1-0. 35578E-01-0. 372106-0 1-0. 382 886-01  W  3  5 

-0.389236-01-0.3922*6-01-0.392966-01-0.392*26-01— 0.391*96-01-0. 39067E-01  *36 


7.792500  2.329599  -0.071767  * 

0.100006  01  0.0  0.557356  01  0.0  0.655566  01  0.0  *  1 

0.691306  00  0.0  0.389256  00  0.0  0.620976-02  0.0  *  2 

-0.976506-01  0.0  0.560*16  00  0.0  0.3722*6  00  0.0  *  3 

-0.126806-01  0.0  -0.106026  01  0.0  *  * 


0.150266-01  0. 151606-01  0.155706-01  0.162316-01  0.171166-01  0.181936-01  U  *  1 

0.19*266-01  0.20754E— 01  0.221006-01  0.23398E-01  0.2*6126-01  0.257126-01  U  *  2 

0.266666-0!  0.27650E-01  0.280626-01  0.28529E-01  0.288966-01  0.2921*6-01  U  *  3 

0.295166-01  0. 297966—01  0.30020E-0I  0.301336-01  0.300696-01  0.297*86-01  U  *  * 

0.291056-01  0.281166-01  0.26800E-01  0.252276-01  0.235006-01  0.217*86-01  U  *  5 

0.200896-01  0.186106-01  0.173576-01  0.16356E-01  0.15626E-01  0.151816-01  U  *  6 

-0.208176-01-0.207716-01— 0.20705E-0 1-0.20587E-C 1-0.203956-0 1-0. 201 166-01  V  *  l 
-0.197336-01-0.192196-01-0. 185716-01-0. 178216-01-0. l701*E-01-0. 161756-01  V  *  2 

-0.1532*E— 01-0. l**976-0l-0. 137*16-0 1-0. 1311 3E-01-0. 126686-0 1-0.12**96- 01  V  *  3 

-0.12*82E-01-0.12775E— 01-0.13331E-01-0. 1*1*16-01-0.1 5 17*6-01-0. 163616-01  V  *  * 

-0.17603E-01-0.18787E— 0I-0.19803E-01-0.20572E-01 -0.21073E-01-0.213*06-01  V  *  5 

-0.2l*30E— 01-0.21 390E-01-0.21268E-0 1-0.21 l l*E-0 1-0. 209766-0 1-0. 208 776-01  V  *  6 

-0.2993 7E-01-0.29868E— 01— 0.29637E-0 1-0 .2927*6-01-0. 281 38E-0 l- 0.280 136-01  W  *  1 

-0.27073E-01-0. 259016-01-0. 2*50 8E-0 1-0.229**6-01-0.212676-01-0. 1952*6-01  W  *  2 

-0.177696-01-0.160786-01-0. 1*5*16-01-0. 1 3247E-01- 0.l2276E-0i-0. 116856-01  M  *  ? 

-0.1 1 500E-01-0. 1 1 7-*0E— 01-0. 12*  18E-01-0. 135326-01-0.150*06-01-0.  1686’iE-Ol  8  4* 

-0.18875E-01 -0.209*16-01-0. 229156-01-0. 2*6816-01-0.2617*6-01-0. 27374E-01  8*5 
-0.28291E-01-0.28951E-01-0.29394E-01-0.29676E-01-0.29845E-01-0.29931E-01  8*6 


FIGURE  7.  (Continued) 


10.94499*  2.1)2000  *0.071747  5 

0.100006  01  0.C  0.940006  01  0.0  0.301576  01  0.0  S  1 

0.374416  00  0.0  0.273076  00  0.0  0.390496*02  0.0  5  2 

—0.573706*01  0.0  0.301376  CO  0.0  0.103206  00  0.0  5  3 

-0.571126-02  0.0  *0.437056  00  0.0  54 


0.152546-01  D.  153445*01  0.154276*0  1  0.140776-01  0.144476-01  0.173706*01  U  5  1 

0.101546*01  0.109756-01  C.197046-01  0.205476-01  0.212496-01  0.210796-01  052 
0.2 24256-01  0.220756-01  0.232326-01  0.235106-01  0.237336-01  0.239206-01  053 

0.241106-01  0.242746-01  0.244056-01  0.244746-01  0.244476-01  0.242006-01  O  5  4 

0.239406-01  0.234436-01  0.227356-01  0.210556-01  0.200506-01  0.197046-01  U  5  5 

0.107306-01  0.177476-01  0.149146-01  0.142166-01  0.156956-01  0.153716-01  U  5  6 

-0.211366-01-0.210776-01-0.209596-01-0.207596-01-0.204476-01-0.200926-01  V  5  1 

— 0. 19596C— (tt— 0. 190026-01-0. 103136-01-0. 175456-01-0.147946-0 1-0. 160296-01  V  5  2 

— 0. 152916—01*0. 145786-01-0. 139536-01-0. 1 34436-01-0. 1 30856—0 1-0.129096-01  V  5  3 

—0.129276—01—0.131446—01—0.135446-01-0.141746—01-0.149546-01-0.158616-01  V  5  4 

-0.148346-01-0.179056-01-0.187006-01-0.194636-01-0.200676-01-0.205186-01  V  5  5 

-0.208326-01-0.210286-01-0.211286-01-0.211636-01-0.211686-01-0.211606-01  V  5  6 

-0.212266-01— 0.211496— 01-0.209426-01-0. 205976— 01-0. 2011 16-01-0. 194816-01  U  5  l 
—0. 187056-01— 0. 177866-01-0. 16748 E—0 1-0. 156306—01— 0. 14468E-01-0. 132906-01  8  5  2 

-0.121326-01-0.110376-01-0.100556-01-0.923336-02-0.861486-02-0.822716-02  853 
-0.80838E-02-0.81932E-02-0.85626E-02-0.91900E-02-0.100546-01-0.111 156-01  854 
-0.123196-01-0.135996-01-0.148896-01-0. 161266-01-0.172686-01-0.182786-01  855 
-0. 191366-01-0. 198296-01-t. 203636-0 1-0. 207566-01-0.210266-01-0. 211826-01  8  5  6 


13.000000  1.955899  -0.071767  6 

0.100006  01  0.0  0.579036  01  0.0  0.321146  01  0.0  6  l 

0.290796  00  0.0  0.  93436  00  0.0  0.259096-02  0.0  6  2 

-0.342086-01  0.0  8.164836  00  0.0  0.91926E-01  0.0  6  3 

-0.262906-02  0.0  -0.184576  00  0.0  6  4 


0.146696-01  0.147286-01  0.14914E-01  0.152096-01  0.15591E-01  0.160396-01  U  6  1 

0.165316-01  0.170366-01  0.175256-01  0*179806-01  0.183946-01  0.187676-01  U  6  2 

0.190896-01  0.193576-01  0.195726-01  0.197436-01  0.198836-01  0.200066-01  U  6  3 

0.201216-01  0.202256-01  0.203076-01  0,203536-01  0.203466-01  0.202656-01  U  6  4 

0.200896-01  0.198046-01  0.194046-01  0.188946-01  0.182966-01  0.176456-01  U  6  5 

0.169876-01  0.163636-01  0.158026-01  0.153316  01  0.149756-01  0.147516-01  U  6  6 

-0.196496-01-0.195906-01-0.194556-01-0.192326-01-0.189206-01-0.185236-01  V  6  1 

-0.180416-01-0.174796-01-0.168536-01-0.161956-01-0.155356-01-0.148906-01  V  6  2 

-0. 14275E-01-0. 137066-01-0. 132086-0 1-0.128066— 01:-0.l2524E-0 1-0. 123826-01  963 
-0.123906-01-0.125486-01-0.128546-01-0.133056-01-0.138856-01-0.145656-01  V  6  4 

-0.153076—01-0. 160676-01-0.167966-01-0. 174576-0 1-0. 1802 8E-0 1-0. 185046-01  V  6  5 

-0.188876-01-0.191776-01-0. 193816-01-0. 195176-01-0.196036-01-0.196496-01  V  6  6 

-0.148666-01-0. 148066-01 -0.146356-0 1-0. 143546-01-0. 139666-01-0.1 34 76E-01  W  6  l 
-0.128896-01-0.122136-01-0.114726-01-0.106906-01-0.988976-02-0.908876-02  to  6  2 

-0.830896-02-0.757826-02-0.692656-02-0.638226-02-0.596936-02-0.570326-02  M  6  3 

- 0.559126-02-0. 563836-02-0. 534936-0 2-0. 62230E-02-0. 674736-02-0 • 74004E- 02  M  6  4 

-0.81 5416-02-0. 897526-02 -0.982876-0 2-0. 10682F-01-0. I 15086-0 1-0. 122 776-01  W  6  5 

-0.129636-01-0.135486-01-0.140246-01-0.143936-01-0.146576-01-0. 148156-01  H  6  6 

FIGURE  7.  (Continued) 


28 


U.IIMW  1.737200  -0.0717*7  7 

0.100806  01  0.0  0.5917*6  01  C.O  0.253326  01  0.0  7  l 

0.203726  00  0.0  0.1203*6  00  0.0  0.1*3106-02  0.0  7  2 

“0.1*7906—01  0.0  0.710526-01  0.0  0.355906-01  0.0  f  3 

—Or 903996  03  0.0  -0.5*3*06-01  0.0  7  4 


0.133126-01  0.133*36-01  0.13*476-01  0.13*096-01  0.130176-01  0.1*0506-01  U  7  1 

0.1*3106-01  0.1*5016-01  0.1*0316-01  0.150*26-01  0.152726-01  0.154*16-01  U  7  2 

0.15*2*6-01  0.157*56-01  0.150706-01  0.159706-01  0.160*76-01  0.1*1166-01  U  7  3 

0.1*1016-01  0.1*2*06-01  0.1*2076-01  0.1*3176-01  0.1*3206-01  0.162076-01  U  7  * 

0.1*2076-01  0.1*0726-01  0.150706-01  0.15*236-01  0.153166-01  0.1*9756-01  U  7  5 

0.1*6226-01  0.1*2006-01  0.1396*6-01  0.136906-01  0.13*916-01  0.133616-01  U  7  6 

— 0.160076-01— 0.160366— 01— 0.1671 16-0 1-0. 1650 06— 01-0.162316-01-0.150666-01  V  7  1 

-0.15*016-01 -0.150236-01-0.1*5276-01-0. 1*0176-01-0.135156-01-0. 130326-01  7  7  2 

—0.125706— 01-0.12 163E— 01-0.1100*6-01-0.115156— 01-0.113136—01-0.112006-01  V  7  3 

-0.112066-01-0. 113076-01-0. 1150W-0 1-0. 118106-01-0. 121996-0 1-0. 126606-01  V  7  * 

— 0. 131716-01-0. 137056-01-0. 1*2356-0 1-0. 1*7376— 01— 0. 151996-01-0. 156126-01  V  7  5 

-0.159716-01—8. 1*2606-01-0.1*5026-01—0. 16*776—01—0. 168006-01-0.160726-01  V  7  6 

-0.907006-02— 0.903996-02-0. 092156-0 2-0. 87200E— 02-0.0*60 06-02-0. 01***€-02  W  7  1 

—0.776306— 02-0.733616-02-0.607626-02-0.639066— 02— 0.591396-02-0.5*3136-02  U  7  2 

—0**96376— 02— 0.*52726-02-0.*138*6-0 2-0. 301286— 02-0. 356286— 02-0.339616-02  U  7  3 

-0.331636*02— 0.332*56-02-0.3*2976-02-0.362566— 02— 0.390766-02-0. *26516-02  W  7  * 
—0.*60526— 02-0.515296— 02-0.565296-02— 0.617076— 02-0.669106-02-0.719*06-02  8  7  5 
-0.766136-02-0.807*06-02-0.0*2336-02-0.070356-02-0.090926-02-0.903526-02  H  7  6 


20.125000  l.***599  -0.071767  0 

0.100006  01  0.0  0.600676  01  0.0  0.175186  01  0.0  8  1 

0.117156  00  0.0  0.575556-01  0.0  0.569376-03  0.0  0  2 

-0.555216-02  0.0  0.197596-01  0.0  0.813056-02  0.0  0  3 

-0.171916-03  0.0  -0.8913*6-02  0.0  8  * 


0.111506-01  0.111636-01  0.112066-01  0.112736-01  0.113576-01  0.11*5*6-01  U  8  l 

0.115586-01  0.116616—01  0.117586-01  0.110*76-01  0.119296-01  0.1230*6-01  U  8  2 

0.120706-01  0.121276-01  0.1217*6-01  0.1221*6-0;  0.122*96-01  0.122806-01  U  8  3 

0.123116-01  0.123396-0*  0.123626-01  0.123776-01  0.123036-01  0.1237*6-01  U  8  * 

0.123*76-01  0.122996-01  0.122256-01  0.121266-01  0.1200*6-01  0.118656-01  0  8  5 

0.117196-01  0.1157*6-01  0.11**06-01  0.113226-01  0.112316-01  0.111736-01  U  8  6 

— 0.131506-01— 0.131156— 01-0.1302*6— 01-0.128786-01-0. 126826-0 1-0. 12** 26-01  V  8  1 

-0.121676-01-0.110626-01-0.115386-01-0.112116-01-0.108926-01-0.105886-01  V  8  2 

-0.103036-01— 0.100*6 6-01-0. 982506-02-0.96*776-02-0.9522 56-02-0. 9*55*6-02  V  8  3 

-0.^4*876-02— 0.950206-02 -0.961*56— 02-0.978*56-02— 0.100076-01-0.102 736-01  V  8  ♦ 

-0.105716-01-0.108896-01-0.112126-01-0.115316-01-0.118366-01-0.121236-01  V  8  5 

-0.123856-01-0.1261*6-01-0.1280*6-01-0.1295*6-01-0.1306*6-01-0.131316-01  V  8  6 

—0.*30'?76— 02— 0.*29036-02—0. *22 706-0 2-0. *12*96-02-0. 398856-02-0. 382 136-02  8  8  1 

—0.36273£ -02—0. 3*1266-02-0. 3 18*96-02-0. 2950 36-02-0. 27 12 16-02-0. 2*7*36-02  W  8  2 

-0.22*3*6— 02— 0.202726— 02 -0. 183396-02-0. 167056-02-0 15*296-02-0.1*5*16-02  8  8  3 

-0.1*0576-02-0.139956-02-0.1*3756-02-0.151956-02-0.16*276-02-0.180256-02  88* 
-0.199*36-02-0.221266-02-0.2*5226-02-0.270816-02-0.297376-02-0.323866-02  885 
-0.3*9066-02-0.371986-02-0. 39191E-02-0.*08306-02-0.*20566-02-0.*28186-02  886 


FIGURE  7.  (Concluded) 


•  Punched  Output: 

The  punched  card  output  for  the  sample  problem  is  shown  in  tabulated  form  in 
Figure  7.  The  output  data  block  for  the  first  spanwise  station  is  identified. 
The  first  card  lists  the  spanwise  station  Y  =  0. ,  the  mapping  radius 
R  =  2.8889  2nd  the  rate  of  change  of  R  with  Y,  DRDY  =  0.  The  next  four 
cards  list  the  real  and  imaginary  parts  of  the  coefficients  used  in  the  map¬ 
ping  expansion.  Cards  6  —  11  list  the  induced  velocity  components  in  the 
X  direction  for  each  cf  the  36  control  points  at  Y  =  0.  The  induced  velocity 
components  in  the  Y-direction  are  listed  on  cards  12  —  17  and  cards  18-23 
specify  induced  velocity  components  in  the  Z -direction.  Data  blocks  of  this 
type,  each  consisting  of  23  cards,  follow  for  each  of  the  other  7  spanwise 
stations  specified  as  part  of  the  input.  The  punched  output  is  identified  in 
columns  73  —  80.  The  spanwise  station  number  is  shown  in  columns  75  —  77. 
Sequence  numbers  for  each  station  appear  in  columns  78  —  80.  The  letters 
U,  V,  W  in  column  74  identify  the  velocity  components  listed  on  the  data 
cards. 

Note:  From  the  tabulations  of  Figure  7  it  is  apparent  that  the  first  five 
cards  of  the  data  generated  for  each  spanwise  station  represent  an  exact 
duplication  of  input  cards  described  previously.  They  are  generated  as 
part  of  the  punched  output  so  that  a  more  complete  data  package  for  the 
Transformation  Method  program  may  be  obtained  and  additional  card 
handling  circumvented. 

(2)  Applicability  and  Limitations 

The  Jet  Flow  Field  program  may  be  utilized  to  evaluate  the  induced  flow  field  at 
given  control  points  due  to  one,  two  or  three  exhausting  jets.  For  a  single  jet  the 
initial  jet  exhaust  direction,  specified  by  <t>  and  4> ,  and  the  freestream  direction,  speci¬ 
fied  by  a  and  are  arbitrary.  For  a  two-jet  configuration  the  jet  exits  must  both  lie 
in  the  same  XY  plane  and  the  jet  exhaust  planes,  defined  by  the  freestream  vector  and 
the  initial  jet  exhaust  Vectors,  must  be  parallel.  The  same  restrictions  apply  to  a 
three-jet  configuration.  Additionally,  three- jet  configurations  must  be  colinear  and 
negative  angle -of-attack  cases  cannot  be  treated.  More  complex  configurations,  as 
shown  below,  may  be  treated  by  reduction  to  two-  or  three-jet  type  configurations,  as 
indicated,  and  adding  the  induced  velocities  at  each  control  point  due  to  configurations 
(ai  and  (b). 
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Four-Jet 

Configuration 


©  ®  @ 


Six-Jet 

Configuration 


Extensive  comparisons  between  computations  and  experimental  data  have  been  made 
for  velocity  ratios  0. 10  <  Uoo/U^o<  0.30  and  the  Jet  Flow  Field  program  may  be 
considered  most  applicable  in  this  range  of  velocity  ratios. 

The  choice  of  the  variables  governing  the  numerical  integration  for  the  jet  path 
is  related  to  the  velocity  ratio  of  the  problem  being  considered.  For  U^/U^  <  0. 125 
integration  in  the  direction  normal  to  the  freestream  over  an  extent  of  at  least  30  jet 
exit  diameters  has  been  found  desirable.  As  U^/U^  increases  this  may  be  decreased, 
as  the  jet  penetrates  less  at  the  higher  velocity  ratios.  For  the  above  range  of  velocity 
ratios  an  integration  step  size  of  <0.5  jet  exit  diameters  has  been  found  satisfactory. 

Control  points  at  which  jet-induced  velocity  components  are  to  be  evaluated  may 
not  lie  withi  n  the  jet  itself,  as  the  Jet  Flow  Field  theory  is  not  valid  in  this  region. 
Generally,  control  points  positioned  less  than  2  jet  exit  diameters  from  the  center  of 
the  jet  exit  should  be  avoided,  to  avoid  distortion  in  the  computed  velocity  distributions. 

b.  Lifting  Surface  Theory 

The  purpose  of  the  Jet  Flow  Field  theory,  when  used  in  conjunction  with  the 
Lifting  Surface  theory,  is  to  predict  jet-induced  downwash  distributions  on  the  wing 
to  be  utilized  by  the  Lifting  Surface  theory  in  evaluating  power  effects.  This  is  accom¬ 
plished  by  executing  the  Jet  Flow  Field  computer  program  to  generate  required  input 
data  for  the  Lifting  Surface  program  in  the  form  of  punched  data  cards.  These  data 
cards  will  then  constitute  the  downwash  matrix  [w] ,  which  forms  part  of  the  input  for 
the  Lifting  Surface  program  described  in  Section  II.  5. 

It  should  be  noted  that  the  manner  in  which  the  Jet  Flow  Field  program  is  utilized 
to  provide  dat?  for  the  Lifting  Surface  program  is  almost  identical  to  its  application  in 
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conjunction  with  the  Transformation  Method  described  in  Section  n.  3.  a.  The  discussion 
below  will  treat  in  detail  only  those  areas  of  input  and  output  which  differ  from  Section  XL  3 .a 

(1)  Sample  Problem  Computation 

For  the  sample  problem  being  considered,  the  Jet  Flow  Field  program  is  now 
used  to  compute  jet-induced  downwash  distributions  at  10  spanwise  stations  on  the 
wing.  Figure  8  shows  details  of  the  planform  of  the  wing  and  indicates  the  network  of 
control  points  to  be  utilized  in  the  computations. 


FIGURE  8.  CONTROL  POINTS  ON  WING  FOR  SAMPLE  PROBLEM 

The  wing  is  treated  as  a  planar  surface,  Z  =  0.  Each  spanwise  station  has  10  control 
points  spaced  evenly  between  0. 05  and  0. 95  of  the  local  chord.  The  spanwise  stations 
are  distributed  evenly  between  0. 05  and  0. 95  of  the  semi-span. 

(a)  Input  for  Sample  Problem 

The  input  cards  required  for  the  sample  problem  are  tabulated  in  Figure  9. 

Card  1  again  lists  the  three  control  indices.  IGE0M  =  3  indicates  that  the 
coordinates  of  all  control  points  will  be  provided  directly  as  part  of  the  input.  The 
other  two  indices  remain  unchanged  from  Section  II.  3.  a.  The  punched  output  generated 
will,  of  course,  now  be  suitable  for  use  with  the  Lifting  Surface  program. 

Card  £  remains  unchanged  from  Section  II.  3.  a. 
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G  JRE  9.  JET  FLOW  FIELD  PROGRAM  INPUT  DATA  FOR  SAMPLE  PROBLEM 

(Wing;  Lifting  Surface  Theory) 


Card  3  controls  the  number  of  steps  and  the  step  size  in  the  numerical  integra¬ 
tion  of  the  equations  of  motion  for  the  jet  path.  For  this  case,  the  number  of  steps 
has  been  cut  down  to  50,  as  computations  in  Section  n.  3.  a  showed  that  after  a  pene¬ 
tration  of  20  jet  exit  diameters  into  the  crossflow,  the  nondimensionalized  jet  speed 
U/U  =  0.202,  i.e. ,  the  jet  has  virtually  slowed  to  the  speed  of  the  crossflow  and 
further  contributions  to  the  jet-induced  velocities  will  be  negligible. 

Cards  4  —  6  see  Section  n.  3.  a. 

Card  7  lists  the  number  of  spanwise  icntrol  stations,  NS  =  10  and  specifies  that 
the  number  of  control  points  at  each  station  NC  =  10. 

Cards  8  —  57  list  the  coordinates  for  the  control  points  of  the  grid  shown 
schematically  in  Figure  8.  The  coordinates  for  each  control  point  appear  in  the  order 
X,  Y,  Z.  Control  points  are  listed  from  leading  edge  to  trailing  edge  for  each  span- 
wise  station,  with  the  spanwise  stations  appearing  in  a  root-to-tip  sequence.  The  total 
number  of  control  points  is  NC  x  NS.  The  listing  is  continuous,  i.e. ,  no  new  record 
is  required  for  the  first  control  point  at  each  spanwise  station. 

(b)  Output  for  Sample  Problem 

Both  printed  and  punched  outputs  are  obtained. 

•  Printed  Output : 

The  initial  part  of  the  printout  dealing  with  configuration  identification  and 
jet  centerline  printout  will  be  identical  to  that  shown  in  Figures  6(a)  and  6(b) 
of  Section  II.  3.  The  centerline  computations  are  now  carried  out  to  Z  =51.63 
consistent  with  the  50  integration  steps  specified.  Figure  10  shows  a  portion 
of  uie  printout  for  the  jet-induced  velocity  components  at  the  control  points 
specified  as  part  of  the  input.  The  control  points  are  listed  in  the  order  in 
which  they  were  read  in.  All  three  velocitv  components  are  printed  out, 
although  only  the  downwai  h  component  W  needed  in  conjunction  with  the 
Lifting  Surface  theory. 

•  Punched  Output : 

The  punched  output  for  the  sample  problem  is  shown  in  tabulated  form  in 
Figure  11.  The  data  block  for  the  first  spanwise  station  is  identified.  The 
cards  list  the  downwash  component  -W,  nondimensionalized  by  U  ,  for  each 
control  point  at  the  first  spanwise  station  in  a  leading  to  trailing  edge 
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1.006 


*** 

INOUC60 

1 

V6LOC1T16S  AT  CONTROL  POINTS  *** 

X 

V 

2 

u 

V 

U 

0.696 

1.006 

0.0 

0.  421856-0^ 

-0.25052E-02 

-0.194706-01 

1.742 

1.006 

0.0 

0.410696-01 

-0.31934E-02  , 

-0.273646-01 

2.788 

1.006 

o.o 

0.38306E-01 

•0.375686-02 

-0.342896-01 

3,835 

1.006 

0.0 

0.344886-01 

-0.414906-02 

-0. 397426-01 

4.881 

1.006 

0.0 

0.302336-01 

-0.437116-02 

-0.436256-01 

5.927 

1.006 

0.0 

0.260046-01 

-0.445366-02 

-0.461 186-bl 

6.973 

1.006 

0.0 

0.220696-01 

-0. 44349 E-02 

-0.475076-01 

8.019 

1.006 

0.0 

0.18547E-01 

-0.434956-02 

-0.480746-01 

9.066 

1.006 

0.0 

0. 15463E-01 

-0.42236E-02 

-0.480556-01 

10.112 

1.006 

0.0 

0.127976-01 

-0.40755E-02 

-0.476306r0l 

1.015 

3.018 

0.0 

0.395486-01 

-0.737346-02 

-0.  19404E-01 

2.008 

3.018 

0.0 

0.383486-01 

-0.397896-02 

-0.257786-01 

3.000 

3.018 

o.o 

0.35964E-01 

-0.103026-01 

-0.314106-01 

3.993 

3.018 

0.0 

0.327956-01 

-0.112666-01 

-0.359956-01 

4.986 

3.018 

0.0 

0.29248E-01 

-0.1 1869E-0 l 

-0.394556-01 

5.978 

3.018 

0.0 

0.25651E-01 

-0.12161E-01 

-0.418756-01 

6.971 

3.018 

0.0 

0.22215E-01 

-0.122116-01 

-0.434186-01 

7.963 

3.018 

0.0 

0.190566-01 

-0 . 12087E-0 1 

-0.44267E-01 

8.956 

3.018 

0.0 

0.J6219E-01 

-0. 1 1845E-0 1 

-0.445866-01 

9.949 

3.018 

0.0 

0. 13710E-01 

-0.115296-01 

-0.445146-01 

FIGURE  10.  INDUCED  VELOCITY  COMPONENTS  AT  CONTROL  POINTS 
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0.19469976-01 
0. 46118346-01 
0.19404246-01 
0.41874816-01 
C. 172291 06-01 
0.35073546-01 
0.14135876-01 
0.27716086-01 
0.11017706-01 
0.21044606-01 
0. 82872476-02 
0.15538076-01 
0.60588686-02 
0.11229796-01 
0.43101016-02 
0.79637516-02 
0.29697406-02 
0.55347086-02 
0.195792 76-02 
0.37495196-02 


0.27363576-01 
0.4750663E-01 
0.2577788E-01 
0.43418066-01 
0.21926766-01 
0.,36746506-01 
0.17456146-01 
0.29336516-01 
0.1332559E-01 
0.22460606-01 
0.9881828E-02 
0.16687146-01 
0.7158797E-02 
0.12116696-01 
0.50688916-02 
0.86256646-02 
0.34934536-02 
0.60174546-02 
0. 23193876-02 
0.40957486-02 


0.34289076-01 
0.48073736-01 
0.31409716-01 
0.4426672E-01 
0.26182716-01 
0.37892026-01 
0.20547006-01 
0.30601866-01 
0.15523176-01 
0.23660596-01 
0. 1142634E-01 
0.1771446E-01 
0.82375266-02 
0.12938906-01 
0.58197046-02 
0.92548846-02 
0.40148536-02 
0.64844976-02 
0.26808266-02 
0.44349136-02 


0.39741876-01 

0.4805501E-01 

0.35995056-01 

0.44586296-01 

0.29826826-01 

0.38601386-01 

0.2331*986-01 

0.3155021E-01 

0.17562436-01 

0.2465511E-01 

0.12896096-01 

0.18620016-01 

0.92824776-02 

0.13693736-01 

0.6556150E-02 

0.98488636-02 

0.45308256-02 

0.69340326-02 

0.30406656-02 

0.4*659240-02 


0.4362537E-01 
0.47630486-01 
0. 3945525E-01 
0.44513896-01 
0.32787346-01 
Q. 38960046-01 
0.25712726-01 
0.32223506-01 
0. 1940951E-01 
0.25459556-01 
0.14271436-01 
0. 1940690E-01 
0.10282816-01 
0.14380106-01 
0.72724896-02 
0. 10405786-01 
0. 50383516-02 
0. 7364653E-02 
0.3397372E-02 
0.50878046-02 


FIGURE  11.  JET  FLOW  FIELD  PROGRAM  PUNCHED  OUTPUT 
FOR  SAMPLE  PROBLEM 
(Wing;  Lifting  Surface  Theory) 
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sequence.  The  sign  of  W  is  changed  to  provide  compatibility  with  the 
Lifting  Surface  theory  where  downwash  is  conventionally  considered  to 
be  positive.  Similar  data  blocks  are  generated  for  the  other  9  spanwise  j 

stations.  Each  spanwise  station  starts  on  a  new  record. 

(2)  Applicability  and  Limitations 

See  discussion  in  Section  II.  3.  a. 

Additionally,  since  the  punched  output  variable  -W  serves  as  an  approximation  i 

to  the  tangent  of  the  downwash  angle  when  it  is  utilized  as  input  to  the  Lifting  Surface 
program,  the  application  of  the  Jet  Flow  Field  program  must  be  restricted  to  small 
angles  of  attack. 


4.  APPLICATION  OF  TRANSFORMATION  METHOD  TO  WING 

The  transformation  method  uses  the  jet-induced  velocity  components  at  a  number 
of  stations  on  the  wing  to  determine  wing  power  effects  in  the  form  of  surface  pressure 
distributions  and  integrated  force  and  moment.  The  transformation  method  requires 
that  the  mapping  function  for  each  of  the  wing  sections  is  known.  The  jet  induced 
velocity  components  are  determined  using  the  jet  flow  field  program  described  in 
Volume  T.  The  mapping  function  is  determined  using  the  techniques  developed  in 
Section  III  of  Volume  I. 

The  generation  of  the  coefficients  of  the  mapping  function  for  the  sample  problem 
has  been  described  in  Section  II.  2.  A  description  of  the  application  of  the  jet  flow 
field  program  to  the  sample  problem  has  been  given  in  Section  II.  3.  The  punched  card 
output  of  the  jet  flow  field  program  is  compatible  with  the  transformation  method  and 
includes  the  mapping  coefficients  and  jet  induced  velocity  components  for  each  of  the 
wing  sections.  To  complete  the  input  to  the  transformation  method  program  various 
flow  indices  must  also  be  specified.  Most  of  these  indices  have  already  been  defined 
in  the  preceding  section. 


a.  Inputs  to  Transformation  Method  for  Sample  Problem 

Shown  in  Figure  12  are  the  input  data  for  this  sample  problem  in  which  the 
punched  outputs  from  the  jet  flow  field  program  constitute  the  major  portion.  However, 
to  activate  the  computation  two  cards  must  precede  this  basic  input  block.  There  may 
be  none,  two  or  three  cards  following  this  block,  depending  on  the  specified  options. 
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Card  1  lists  iii  order  the  classification  index  (1  specifies  a  wing),  the  modi¬ 
fication  index  (1  denotes  the  option  of  three-dimensional  modification  being  exercised), 
the  number  of  iterations,  the  number  of  layers  for  distributing  residual  sources  and 
sinks,  the  power  index  (0  indicates  the  power  effect),  the  configuration  index  (1  refers 
to  a  nonrectangular  wing),  and  the  force  index  (1  indicates  forces  and  moments  to  be 
computed). 

Card  2  lists  in  order  the  number  of  stations,  the  number  of  pairs  of  the 
mapping  coefficients,  the  number  of  coefficients  in  the  Fourier  series  expansion,  the 
computation  index,  the  number  of  angular  increments  on  the  mapping  circle,  the  free- 
stream  to  the  jet  velocity  ratio,  the  angle  of  attack  in  degrees,  and  the  sideslip  angle 
in  degrees. 

Cards  3  through  186  contain  the  punched  output  data  provided  by  the  jet 
flow  field  program,  which  include  the  y  coordinate,  the  mapping  coefficients,  and  the 
induced  velocity  components  for  stations  No.  1  through  No.  8.  There  are  23  cards 
for  each  station. 

Card  187  lists  in  order  the  option  index  (0  denotes  no  average  value  usf  d  for 
‘he  station  next  to  the  exhausting  jet)  and  the  station  number  where  the  jet  is  located. 

Card  188  lists  the  number  of  stations  on  which  the  downwash  modification  is 
to  be  applied. 

Card  189  lists  in  order  the  number  of  jets,  the  jet  exit  diameter,  X  coordinate  of 
the  moment  center,  Z  coordinate  of  the  moment  center,  and  the  reference  length  for 
making  the  computed  moments  dimensionless. 

b.  Outputs  from  Transformation  Method  for  Sample  Problem 

Figure  13  lists  directly  or  indirectly  a  portion  of  the  input  data  on  Card  1 
through  Card  186. 

Figure  14  establishes  the  correspondence  between  the  angular  increments  of  the 
mapping  circle  and  their  corresponding  locations  on  the  wing  section  at  every  station. 
The  first  column  states  the  angular  increments  m  degrees. 

Figure  15  gives  the  pressure  distributions  in  coefficient  form  (p-  p^  ,  at 
every  station  after  completion  of  the  segment  method.  These  coefficients  are  tabu¬ 
lated  against  the  angular  increments..  To  obtain  the  actua1  location,  reference  must 
be  made  to  the  previous  figure.  The  second  and  the  third  lines  in  this  table  list  the 
radius  of  the  mapping  circle  (RB)  and  the  gradient  of  this  radius  in  y-direetion  (DRDY). 
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Figure  16  lists  the  pressure  distributions  at  various  sections  after  imposing 
the  residual  sources  and  sinks  in  the  network.  Columns  7,  8,  9  in  this  table  remain 
the  same  as  those  in  Figure  15,  since  the  flow  properties  near  the  wingtip  are  not 
modified. 

Figure  17  lists  the  pressure  distributions  after  completion  of  a  three-dimensional 
modification  of  one  iteration. 

Figures  18  and  19  show  printout  of  the  pressure  coefficients  after  imposing  the 
residual  sources  and  sinks  for  the  second  time  and  the  completion  of  a  three-dimen¬ 
sional  modification  of  two  iterations. 

Figure  20  lists  the  parameters  used  in  the  three-dimensional  modification  and  in 
the  force  and  moment  computations,  originally  read  in  as  input  data  on  Cards  187,  188, 
and  189.  Also  tabulated  are  the  computed  forces  (normalized  to  the  thrust)  and  moments 
(normalized  by  the  thrust  and  reference  length)  on  this  wing  after  two  iterations. 

Figures  21(a)  through  21(c)  show  the  comparison  between  the  computed  pressure 
coefficients  and  wind  tunnel  test  data  at  stations  y  =  5.04",  7.2925",  10.545",  and  16.0". 
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FIGURE  13.  PARTIAL  LIST  OF  INPUT  DATA  FOR  SAMPLE  WING 
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FIGURE  14.  CORRESPONDENCE  BETWEEN  ANGULAR  INCREMENTS  OF  MAPPING 
CIRCLES  AND  CARTESIAN  COORDINATES  OF  WING  SECTIONS 


FIGURE  14.  (Concluded) 
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FIGURE  15.  POWER-EFFECT  PRESSURE  COEFFICIENTS  ON  SAMPLE  WING 
AFTER  APPLICATION  OF  SEGMENT  METHOD 
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FIGURE  17.  POWER-EFFECT  PRESSURE  COEFFICIENTS  ON 
SAMPLE  WING  AFTER  ONE  ITERATION 


FIGURE  18,  POWER-EFFECT  PRESSURE  COEFFICIENTS  ON  SAMPLE  WING  AFTER 
RESIDUAL  SOURCE  AND  SINK  MODIFICATION  (Second  Time) 
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WING  AFTER  TWO  ITERATIONS 


FIGURE  20.  TORCES  AND  MOMENTS  ON  SAMPLE  WING  BY  TRANSFORMATION  METHOD 


POWER-EFFECT  PRESSURE  COEFFICIENTS  ON  SAMPLE 
WING  AT  STATION  Y  =  5.04 


FIGURE  21c.  POWER-EFFECT  PRESSURE  COEFFICIENTS  ON  SAMPLE 
WING  AT  STATIONS  Y  =  10.545  AND  Y  =  16.05 


Comparison  of  forces  for  sample  problem  with  wind  tunnel  test  date  of  Appendix 
I  is  shown  in  Figure  22,  together  with  further  calculations  and  test  data.  The  calcula¬ 
ted  power-induced  lift  follows  the  trend  of  the  test  data.  The  calculated  values  show  a 
greater  loss  in  lift  than  the  test  data.  This  is  consistent  with  the  surface  pressure 
results  shown  in  Figure  21.  The  reasons  for  these  differences  are  not  known  at  the 
present  time. 

c.  Method  Applicability  arid  Limitations 

This  method  is  generally  applicable  to  power  effects  on  the  wing.  In  addition  to 
the  present  configuration,  fairly  extensive  calculations  on  a  rectangular  wing  of  aspect 
ratio  equal  to  3  with  a  modified  NACA  65-010  section  have  also  been  performed 
(Figure  23(b)).  The  jet  (or  jets)  was  (were)  situated  at  the  midspan  of  the  wing  and 
exhausted  directly  from  the  lower  surface.  But  the  chordwise  location  of  the  jet  was 
allowed  to  vary  (three  positions:  20  percent,  50  percent  and  80  percent  of  the  chord 
length  from  the  leading  edge)  and  the  number  cf  the  jets  could  be  one  or  two.  Most 
of  those  calculations  have  been  compared  with  the  wind  tunnel  test  data  which  were 
obtained  at  Northrop  prior  to  the  present  study..  Some  of  these  comparisons  are 
shown  in  Figures  24  and  25. 

When  the  jet  exhausts  directly  from  the  wing  surface,  the  induced  velocity 
distribution  generally  shows  large  and  abrupt  changes  across  the  jet  station.  If  two 
iterations  are  planned  these  iterations  should  be  smoothed  out.  Otherwise,  the 
computed  results  following  the  second  iteration  may  exhibit  unacceptable  oscillations. 

If  the  calculation  is  limited  to  one  iteration  no  smoothing  was  found  to  be  necessary. 

Since  the  lift  jet  did  not  exhaust  directly  from  the  wing  surface  in  this  sample  problem, 
the  input  data  was  not  smoothed,  even  though  two  iterations  were  computed. 

Under  the  present  scheme  all  the  vertical  velocity  components  at  any  given 
station,  regardless  of  the  chordwise  position,  are  reduced  equally  to  a  magnitude  of 
one  third  of  the  lifting  line  downwash  value  (see  pages  101-102  of  Volume  I  for  details).; 
Since  the  aforementioned  procedure  is  somewhat  arbitrary  some  limitations  on  the 
method  presumably  exist.  For  flow  conditions  radically  different  from  the  present 
one,  the  approach  given  here  may  have  to  be  modified. 

The  present  method,  even  when  the  three-dimensional  modification  is  used, 
does  not  include  all  the  three-dimensional  effects..  The  method  is,  in  effect,  a  quasi- 
two-dimensional  one.  Like  the  widely  used  quasi-one-dimensional  approximation  for 
diffusers  and  nozzles,  Us  applicability  is  not  r.s  restricted  as  it  may  appear.  In  every 
example  considered,  the  agreement  between  calculation  and  test  data  is  fairly  satisfactory 
Because  of  the  quasi-two -dimensional  nature,  however,  computations  beyond  two 
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Segment  Plane 
y— Forward  Jet  Station 


FIGURE  23.  CONFIGURATIONS  TESTED  AT  NORTHROP 


iterations  may  not  be  warranted.  In  practice,  one  iteration  is  what  is  usually  needed. 
Two  iterations  have,  nevertheless,  been  carried  out  in  some  selected  problems  and 
also  in  the  sample  problem  here.  This  is  more  for  demonstration  purpose  than  for 
utility. 

The  present  computer  program  is,  in  a  formal  sense,  capable  of  treating  both 
power-on  and  power-off  problems.  However,  for  the  power-off  case,  the  wing  tips 
exert  a  much  larger  influence  upon  the  flow  property  than  in  the  case  when  power  effects 
are  being  calculated.  This  important  three-dimensional  effect  has  not  been  adequately 
accounted  for  by  the  present  method  and  the  calculations  from  it  are  likely  to  be  less 
accurate.  Therefore,  the  lifting  surface  theory  is  recommended  under  such 
circumstances. 
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5.  APPLICATION  OF  LIFTING  SURFACE  THEORY  TO  WING 


Lifting  Surface  theory  may  be  utilized  to  determine  the  load  distribution  and 
aerodynamic  coefficients  for  a  given  arbitrary  planform  and  specified  downwash  dis¬ 
tribution.  The  Lifting  Surface  computer  program  evaluates  power  effects  on  the  wing 
by  considering  a  known,  jet-induced  downwash  distribution. 

There  are  three  main  components  to  the  program  which  may  be  used  together  in 
one  continuous  operation  or  independently.  The  downwash  control  point  matrix  [d]  is 
generated  in  the  first  part  of  the  program,  its  least  squares  inverse  |DJ  is  generated 
in  the  second  part  of  the  program.  The  pressure  distribution  produced  by  the  speci¬ 
fied  downwash  matrix  [w]  is  computed  by  the  third  component  of  the  program.  The 
downwash  control  point  matrix  and  its  inverse  depend  only  on  the  planform,  the  loca¬ 
tion  of  the  downwash  control  points  and  the  number  of  terms  in  the  loading  series. 

Both  matrices  are  independent  of  the  downwash  distribution.  Once  the  inverse  is 
computed  it  forms  an  input  to  the  third  component  of  the  program,  where  the  pressure 
distribution  is  computed.  Thus  [DJ  may  be  retained  in  punched  card  form  and  then 
used  as  input  for  computing  the  pressure  distribution  due  to  specified  downwash  distri¬ 
butions.  The  inverse  need  not  be  recomputed  as  long  as  the  planform,  location  of 
downwash  control  points  and  the  size  of  the  pressure  loading  series  remain  unchanged. 

The  Lifting  Surface  program  is  a  modified  version  of  the  computer  program 
for  designing  and  analyzing  subsonic  lifting  surfaces  documented  in  Reference  12. 

The  design  options  have  been  eliminated.  The  capability  to  calculate  pressure  dis¬ 
tributions  produced  by  a  specified  cambered  surface  has  also  been  deleted. 

a.  Sample  Problem  Computation 

For  the  sample  problem  being  considered  the  Lifting  Surface  program  is  now 
used  to  determine  the  load  distribution  and  aerodynamic  coefficients  on  the  wing, 
produced  by  the  jet-induced  downwash  computed  in  Section  11,3. 

Figure  26  shows  details  of  the  planform  of  the  wing  and  indicates  the  location 
of  downwash  control  points.  Figure  26  is  identical  to  Figure  8  of  Section  II. 3,  except 
that  all  dimensions  are  now  based  on  a  semi-span  of  unity. 
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FIGURE  26.  DOWNWASII  CONTROL  POINTS  ON  WING 


(1)  Input  for  Sample  Problem 

The  input  cards  required  for  the  sample  problem  are  shown  in  Figure  27.  Since 
all  three  main  components  of  the  program,  discussed  in  detail  above,  are  being  exe¬ 
cuted  in  one  continuous  operation,  some  duplication  of  input  data  occurs. 

Card  1  lists  two  control  indices,  specifying  which  of  the  three  major  components 
of  the  program  are  to  be  executed.  The  combination  of  1ST  ART  =  1,  ISTOP  -  3  will 
execute  all  three  major  components.  Consequently,  the  program  will  start  by  com¬ 
puting  the  downwash  control  matrix  [d],  will  find  the  inverse  [p]^and  will  compute  the 
load  distribution  and  aerodynamic  coefficients. 

Card  2  is  a  title  card. 

Card  3  lists  the  number  of  spanwise  stations  on  semi  span  where  downwash  con¬ 
trol  points  are  to  be  located,  NS  =  10.  It  specifies  the  number  of  spanwise  modes  to 
be  used  in  the  pressure  loading  series,  M  =  6  and  the  number  of  chordwise  modes, 

N  -  8.  The  input  control  index  NEED  -  1  indicates  that  the  first  chordwise  mode, 
i.e.,  cot  9/2  mode,  is  to  be  used  in  the  computations.  The  number  of  leading  and 
trailing  edge  flaps  are  specified  with  NFLAP  -  0.  The  next  two  integers  are  print  and 
punch  controls  for  the  downwash  control  matrix  [dJ.  With  NPR  -  0,  NPU  =  0,  no 
printed  or  punched  output  on  [d]  will  be  generated.  The  print  control  NAY  -  0  speci¬ 
fies  that  no  intermediate  print  is  to  be  generated  during  the  computation  of  [d].  The 
number  of  leading  edge  discontinuities  i:  peoified  as  N0LED  -  2,  and  the  number  of 
trailing  edge  discontinuities,  N$TED  -  2,  is  indicated. 

59 


f 


GO 


FIGURE  27.  LIFTING  SURFACE  THEORY  PROGRAM  INPUT  DATA  FOR  SAMPLE  PROBLEM 


O  \S AMPLE 


FIGURE  27.  (Concluded) 


Card  4  indicates  that  the  chord1,  .ise  locations  of  the  downwash  control  points 
must  be  specified  through  input  cards  at  each  spanwise  station  by  listing  SPACE  =  0. 

It  also  lists  the  Mach  number,  FMACH  =  0.,  and  defines  the  root-semichord  F  =  .2667. 

Cards  5  and  6  list  the  spanwise  locations  of  the  downwash  control  points. 

Card  7  may  be  left  blank  for  a  wing  with  no  leading  or  trailing  edge  flaps. 

Card  8  specifies  the  tangents  of  the  sweepback  angles  of  the  leading  edges  of  the 
geometric  regions. 

Card  9  specifies  the  tangents  of  the  sweepback  angles  of  the  trailing  edges  of  lhe 
geometric  regions. 

Card  10  lists  the  spanwise  locations  of  the  leading  edge  discontinuities. 

Card  11  lists  the  spanwise  locations  of  the  trailing  edge  discontinuities. 

Card  12  specifies  the  number  of  downwash  control  points  at  each  spanwise 
station. 

Cards  13-  22  list  the  chordwise  locations  of  the  downwash  control  points  at  each 
spanwise  station  (in  percent  of  local  chord). 

This  completes  the  input  required  for  the  first  main  component  of  the  program, 
which  computes  the  downwash  control  point  matrix  [d]. 

Card  23  is  a  title  card  for  the  next  main  component  of  the  program. 

Card  24  lists  the  number  of  rows  in  the  downwash  control  point  matrix  or  the 
number  of  control  points  contained  in  [d],  NRGfW  =  100.  It  specifies  the  number  of 
columns  in  [d],  NC0L  =  36.  This  is  the  product  of  the  chordwise  and  spanwise  pres¬ 
sure  modes.  The  control  index  NREAD  =  0  indicates  that  the  second  main  component 
of  the  program  is  being  executed  in  a  continuous  operation  and  hence  [d]  will  be  read 
from  a  scratch  tape,  rather  than  input  cards.  With  NPR  =  0,  NPU  =  0,  no  printed  or 
punched  output  will  be  obtained  for  [d]1^  the  inverse  of  the  downwash  control  point 

matrix.  The  print  control  NAY  =  0  specifies  that  no  intermediate  print  is  to  be  gen- 

r  -itf' 

crated  during  the  computation  of  'DJ  . 

This  completes  the  input  for  the  second  main  component  of  the  program,  which 
Inverts  the  downwash  control  point  mutrlx  [d]  to  obtain  [d]^ 

Card  25  is  a  title  card  for  the  third  main  component  of  the  program  which  com¬ 
putes  the  load  distribution  and  aerodynamic  coefficients  for  a  specified  downwash 
matrix. 


Card  26  lists  the  number  of  chordwise  modes  used  in  the  pressure  loading 
series,  N  -  6,  and  the  number  of  spanwise  modes,  M  =  6.  It  specifies  the  number  of 
spanwise 'Stations  where  downwash  control  points  are  located,  NS  =  10,  and  specifies 
the  number  of  rows  in  the  downwash  control  point  matrix,  NR0W  =  100.  It  specifies 

the  number  of  span  wise  stations  where  the  chordwise  pressure  loading  distribution  is 

( 

to  be  calculated,  NETA  =  6.  It  lists  the  number  of  wing  discontinuities,  NDISC  =  2 
and  the  number  of  leading  and  trailing  edge  flaps,  NFLAP  =  0.  The  intermediate  print 
control  is  again  NAY  =  0.  It  also  specifies  the  number  cf  c>r  rdwise  points  at  which 
the  pressure  loading  is  to  be  computed,  NPSI  =  10. 

Card  27  lists  the  number  of  angles  of  attack  to  be  treated,  NALFA  =  1.  It  also 

,  '  l 

specifies  the  number  of  EPSLITsto  be  read  later,  NEPSLN  =  1.  The  input  control 
NEED  -  1  indicates  that  the  first  chordwise  mode,  i.e.,  the  cot 9/2  mode,  is  to  be 
used  in  the  computations.  The  control  index  NREAD1  “0  again  indicates  a  continuous 
operation  and  thus  [d]  will  be  read  from  a  scratch  tape  rather  than  from  input  cards. 
The  next  control  Integer,  NREAD2  *  1  Indicates  that  the  downwash  matrix  DaG  Is  read 
from  input  cards.  The  number  of  downwash  distributions  to  be  analyzed  is  specified 
with  NW  =  1. 

Card  28  specifies  the  root  semi-chord,  F  =  .2667.  It  indicates  that  the  chord- 
wise  locations  of  the  downwash  control  points  are  specified  through  input  cards  at  each 
span  wise  station  by  listing  SPACE  =  0.  It  lists  the  spanwise  location  of  the 
edge  of  the  fuselage,  YF  =  0.  It  indicates  how  the  points  at  which  the  pressure  loading 
is  calculated  are  located  chordwise,  by  giving  the  chordwise  spacing  DPSI  =  .1. 

Card  29  lists  the  spanwise  Coordinates  of  the  downwash  control  point  stations. 

Card  30  specifies  the  spanwise  locations  where  the  pressure  loading  is  to  be 
computed.  , 

Card  31  lists  the  angle  of  incidence  between  the  centerline  of  the  fuselage  and 
wing  root  chord  in  degrees,  EPSLN  =  0. 

Card  32  specifies  the  angle  of  attack  in  degrees,  ALFA  =  0. 

j  i 

Card  33  may  be  left  blank  for  a  wing  with  no  leading'or  trailing  edge  flaps.  , 

Card  34  specifies  the  chord  at  each  spanwise  discontinuity. 

Card  35  gives  the  location  of  each  spanwise  discontinuity. 

Card  36  lists  the  distance  from  root  leading  edge  tc  the  leading  edge  at  each 
spanwise  discontinuity. 


Card  37  specifies  the  number  of  downwash  control  points  at  each  span  wise 
station. 

Cards  38  —  57  specify  the  tangent  of  the  downwash  angle  at  every  control  point. 
Cards  38  —  57  are  the  downwash  matrix  [w]  generated  in  Section  n.3  and  shown  In 
tabulated  form  in  Figure  11  of  Section  IL3. 

(2)  Output  for  Sample  Problem 

With  the  choice  of  the  punch  controls  described  above,  only  printed  output  is 
obtained. 

Figure  28(a)  shows  a  composite  of  the  printout  generated  by  the  first  main  com¬ 
ponent  of  the  program  (CHAIN  1,  8)  which  computes  the  downwash  control  point  matrix. 

Figure  28(b)  shows  the  printout  generated  by  the  second  main  component  of  the 
program  (CHAIN  6,  8)  which  inverts  the  downwash  control  point  matrix  [d].  The 
determinant  of  the  unit  matrix  is  printed  out  as  a  check  on  numerical  accuracy. 

Figure  28(c)  shows  the  output  from  the  third  main  component  of  the  program 
(CHAIN  7,  8)  which  calculates  the  pressure  loading  and  aerodynamic  coefficients. 

Geometric  parameters  of  the  wing  are  shown  in  Figure  28(c)  and  are  all  identi¬ 
fied.  Aerodynamic  coefficients  and  the  pressure  loading  calculated  at  the  span  wise 
stations  specified  are  shown  in  Figure  28(c).  Again  all  computed  variables  are 
Identified. 

b.  Applicability  and  Limitations 

The  program  is  applicable  to  continuous  surfaces  of  arbitrary  planform  and  no 
interference  effects  such  as  slots,  ground  effects,  large  dihedral  angles  or  end  plates 
are  included.  The  program  does  contain  provisions  for  body  effects. 

Downwash  control  points  must  not  be  located  at  or  near  the  leading  edge,  since 
the  cotangent  elements  of  [d]  would  become  excessively  large  and  dominate  the  solu¬ 
tion  for  the  pressure  coefficient  matrix  [a].  Due  to  the  computing  techniques  utilized, 
downwash  control  points  must  not  be  located  at  discontinuities  in  the  planform  and  at 
flap  hinge  lines. 
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CHAIN  11, PI 

CALCULATION  OF  DOWN- ASH 

control 

PCINT  MATRIX  FOR  SAMPLE 

PROBLEM 

NO.  OF  SPANWlSE 

NODES 

»  * 

NO.  OF  CHORD* I SE 

MOOES 

•  6 

NO.  OF  FLAP  MOOES 

*  0 

COTANGENT  MODE* 

NEEO 

«  1 

POSITION  OF  FLAP 

1  -  0 

.0 

100 

OOWNWASH 

control 

POINTS  MACH  NO.  ■O.i' 

OOWNWASH  CONTROL 

POINTS 

1 

TO 

10 

Y» 

0.49999397E- 

■01 

OOWNWASH  CONTROL 

POINTS 

11 

TO 

20 

Y  * 

0.  149999986 

00 

OOWNWASH  CONTROL 

POINTS 

21 

TO 

30 

Y- 

0.2 5000000c 

00 

OOWNWASH  CONTROL 

POINTS 

31 

TO 

40 

Y» 

0. 34999996E 

00 

OOWNWASH  CONTROL 

POINTS 

41 

TO 

SO 

Y* 

0.44999499E 

00 

OOWNWASH  CONTROL 

POINTS 

51 

TO 

60 

Y« 

0.54999O95E 

00 

OOWNWASH  CONTROL 

POINTS 

61 

TO 

TO 

Y* 

0.649996986 

30 

OOWNWASH  CONTROL 

POINTS 

7t 

TO 

PO 

V» 

0. 7500000 JE 

oo 

OOWNWASH  CONTROL 

POINTS 

ttl 

TO 

VO 

v« 

0.64)4999oE 

0  0 

OOWNWASH  CONTROL 

POINTS 

91 

TJ 

100 

Y« 

0.9494I9496 

00 

FIGURE  28(a).  LIFTING  SURFACE  THEORY  PROGRAM  PRINTED 
OUTPUT  I  OR  SAMPLE  PROBLEM 


CHAIN  (6.61 

INVERT  OOWNWASH  CONTROL  POINT  MATRIX  FOK  SAMPLE  PR03LE** 
DETERMINANT  OF  UNIT  MATRIX  ■  0.100000003  01 


FIGURE  28(b).  (Continued) 
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CHAl*  C7,«» 


CALCULATION  OF  PRFS5URE  LOAHNG  OISTR13UTIQN  FOR  5AVLE  «M06LEM 

NO  300V 

GEOMETRIC  PARAMETERS 
AVERAGE  CHORD,  CAVE  «  C.4J0100 

NEAN  AERODYNAMIC  CK3S !"» ,  CoAR  «  0.414331 

LOCATION  OF  1/4  C3AR,  XMAS  ■  0.18J13? 

SPANWlSE  LOCATION  OF  CHAR ,  YBAR  *  0.644314 

RESULTS  F JR  ALF4»  0.0  ,  AND  EPSILON-  0.0  DTGREIS 

LIFT  COEFFICIENT,  CL  »  *0.1101 4 

MOMENT  COEFFICIENT,  Cm  -  -0„0J0o9 

I  NOllCEi)  (WAG  CtlF.FHC  KNT  ,  CD!  «  J.0J095 

PRESSURE  LOADING  0 1  ST* 1 HUT I  UN,  PR 


SPAN  - 

0. 1000 

0.230 J 

3. 3RM0 

0.5230 

J.ROOO 

.0.5300 

FRACTION 

of  chord 
0.1000 

-0.2422 

-0.2237 

-O.  1913 

-0.1376 

-0.0881 

-0.0401 

0.2000 

-0.1884 

-C. 1703 

-0.  1437 

-0. 1 1 83 

-0.0649 

-0.0327 

0.3000 

-0.1439 

-0.1499 

-0. 12H4 

-0.1073 

-0.0405 

-0.C275 

0.4000 

-0. 1 340 

-0.1242 

-0. 1098 

-0.0922 

-0.0520 

-0 .0231 

0.3000 

-0.1072 

-0. 101 » 

-0. 0896 

-0.0756 

-0.0419 

-0.0190 

0.6000 

-0.0834 

-0. OH  1 1 

-0.0737 

-0.0629 

-0.0350 

-O.L’162 

0.7000 

-0 . 0639 

-0.0667 

-0.0630 

-0.0549 

-0.0319 

-0.0147 

0.8000 

-0.0493 

-0.0326 

-0.031  7 

-0.0461 

-0.0282 

-0.0129 

0.9000 

-0.0272 

-0.0301 

-'.'.0307 

-0.0279 

-0.0167 

-0.0084 

1.0000 

-0.0000 

-0.0000 

-0.0000 

-0.0000 

0.0000 

-0.0000 

LOCAL  SEMICH.JSO, 

C/2 

0.2333 

C. 23  33 

0.2130 

0. 1967 

0. 1601 

0.1401 

CL  C/CAVC 

-0.  1841 

-0.1644 

-0. *353 

-0. 1031 

-0.0503 

-  0 . 0  2  2  1 

CM  C**2/CAVt  CHAR 

-0.0132 

-0.003? 

O.OOOrt 

0.0044 

0.0050 

0 • 0028 

LOPC/CAVC 

0.0029 

0.0019 

0.0009 

0.0003 

-0.0000 

-0.0000 

FIGURE  28(c).  (Concluded) 
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c.  Additional  Calculations  and  Comparison  with  Test  Data 

Once  the  inverse  of  the  downwash  control  point  matrix  has  been  obtained  for  a 
particular  wing  planform,  it  is  a  relatively  easy  task  to  obtain  the  power  induced 
aerodynamics  for  a  range  of  downwash  distribution  corresponding  to  different  power 
conditions.  Figures  29a  through  29g  show  calculations  for  foe  test  model,  described 
in  Appendix  I,  compared  with  test  data  for  a  number  of  power  configurations, 
velocity  ratios  and  angles  of  attack. 

The  interference  lift  for  the  vectored  thrust,  forward  position,  90°  nozzle 
deflection  angle  is  shown  in  Figure  29a  for  a  range  of  velocity  ratios.  The  wing 
calculations  show  that,  for  all  velocity  ratios,  the  induced  lift  (lift  with  power  on 
minus  lift  with  power  off)  is  less  than  for  the  static  case  that  is,  UM  =0.  In  con¬ 
trast  the  test  data  for  wing  plus  body  indicate  that  there  is  a  lift  augmentation  at 
the  higher  velocity  ratios.  For  nozzle  deflection  angles  of  45°  (Figure  29b),  although 
there  is  now  no  positive  Jet  interference  lift  at  the  higher  velocity  ratios,  the  induced 
lift  from  the  test  data  is  nearly  constant  for  .  3 £  /U  5  ,  whereas  the  calcula- 

JO 

tions  indicate  the  induced  lift  decreases  as  the  velocity  ratio  increases. 

For  both  the  above  cases  foe  Induced  lift  was  determined  from  foe  test  data 
using  the  inlet  plugged  as  the  unpowered  case.  With  the  inlet  open  foe  maSnstream 
flow  through  the  ejector  produces  a  "power"  effect  which  can  be  quite  lf-rge  as 
indicated  in  Appendix  I. 

The  body  alone  lift,  with  power  on,  for  the  vectored  thrust,  forward  position, 
90°  nozzle  deflection  angle  is  shown  in  Figure  29c.  Since  this  configuration  was  not 
tested  with  the  inlets  plugged,  it  is  not  possible  to  determine  foe  induced  lift. 
However,  due  to  model  symmetry,  it  is  reasonable  to  assume  that  the  power  off 
plugged  inlet  lift  will  be  quite  small  and  so  that  the  induced  lift  graph  for  the  body 
alone  configuration  will  be  very  similar  to  Figure  29c. 

The  positive  lift  arises  due  to  power  induced  uploads  on  the  nacelles  and  body 
ahead  of  the  Jet  exits.  If  this  lift  increment  is  subtracted  from  the  wing-body  test 
data,  we  get  very  good  agreement  with  foe  wing  calculations. 

Figures  29d  and  29e  show  comparisons  between  the  calculations  and  test  data 
for  the  aft  nozzle  positions  with  two  vectoring  angles.  For  the  aft  position  the  body 
alone  power  effects  are  expected  to  be  smaller  than  for  the  forward  position  due  to  the 
jet  exhausts  being  farther  removed  from  the  nacelles.  The  agreement  between  test 
and  calculations  is  seen  to  be  extremely  good  for  all  velocity  ratios. 

Figure  29f  shows  calculations  of  the  induced  lift  on  the  wing  with  the,  lift  jet 
operating  at  a  velocity  ratio  of  .  20  for  a  range  of  or.  For  ffs8°,  Li / ^  and  hence 
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Cjj  (induceu  lift  coefficient),  is  effectively  independent  of  a.  For  larger  a  the 
calculated  induced  lift  is  still  approximately  constant  whereas  the  test  data  shows  a 
sudden  increase  in  induced  lift.  This  is  due  to  a  change  in  the  stalling  characteristics 
for  the  wing,  deduced  from  the  pressure  measurements,  brought  about  by  the  jet 
induced  flow  field  producing  a  downwash  over  the  wing  for  this  particular  jet 
arrangement  and  velocity  ratio. 

This  result  identifies  an  area  in  which  care  must  be  taken  in  using  the  prediction 
methods.  It  has  been  assumed  that  one  can  calculate  power  effects  and  add  these 
to  the  unpowered  aerodynamics  of  the  vehicle.  This  procedure  of  superposition  appears 
to  be  justified  for  the  linear  range  of  a,  but  must  not  be  used  for  nonlinear  a. 

Instead,  these  two  effects  must  be  considered  together  for  nonlinear  a.  It  is  possible 
that  the  induced  flow  field  due  to  the  power  could  be  included  in  the  nonlinear  wing 
aerodynamics  procedure  presented  in  Section  VU  of  Volume  I  but  this  has  not  been 
studied  under  the  present  investigation. 

Similar  observations  may  be  made  for  the  vectored  thrust  configuration. 
Calculations  and  test  data  are  3hown  in  Figure  29g. 

Calculations  of  pitching  moment  due  to  power  effects  show  that  there  is  a 
nose  up  pitching  moment  for  body  alone,  the  magnitude  of  which  does  not  change 
noticeably  with  the  addition  of  the  wing.  This  result  is  in  agreement  with  the 
calculations. 
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FIGURE  29a.  INTERFERENCE  LIFT  FOR  VECTORED  THRUST, 
FORWARD  POSITION,  90°  DEFLECTION  ANGLE ,o  =  (1  =  0 
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FIGURE  29b.  INTERFERENCE  LIFT  FOR  VECTORED  THRUST, 
FORWARD  POSITION,  45°  DEFLECTION  ANGLE, A'  =  /?=0 
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FIGURE  29c.  LIFT  FOR  VECTORED  THRUST,  FORWARD 
POSITION,  90“  DEFLECTION  ANGLE,cv  =  0=O 
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FIGURE  29e.  INTERFERENCE  LIFT  FOR  VECTORED  THRUST, 
AFT  POSITION,  45°  DEFLECTION  ANGLE ,a  =/?  =  0 
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FIGURE  29f.  INDUCED  LIFT  VERSUS  ANGLE  OF  ATTACK, 
LIFT  JET  (U.  /U.q  =  .  20) ,  ;3  =  0 
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FIGURE  29g.  INDUCED  LIFT  VERSUS  ANGLE  OF  ATTACK, 
VECTORED  THRUST,  FORWARD  POSITION,  90° 
DEFLECTION  ANGLE  (Uw  /U.  =.20), 0=0 


SECTION  m 

POWER  EFFECTS  ON  THE  FUSELAGE 

The  calculation  of  the  jet  induced  loads  on  a  fuselage  is  accomplished  by  using 
the  transformation  method  with  the  disturbance  velocities  at  the  surface  of  the  body 
calculated  by  the  jet  program.  To  use  the  transformation  method,  it  is  necessary  to 
map  the  body  at  different  body  stations.  This  section  describes  the  application  of 
these  methods  to  the  calculation  of  fuselage  loads. 

1.  SAMPLE  PROBLEM 

To  demonstrate  the  application  of  the  methods,  the  fuselage  of  the  wind  tunnel 
test  model  which  was  tested  during  this  investigation  will  be  used.  This  fuselage  is 
described  in  Appendix  I  of  this  volume.,  A  sketch  of  the  model  fuselage  with  coordi¬ 
nate  system  is  shown  below. 


U. 

jo 


The  power  and  flight  conditions  which  must  be  specified  to  complete  the  problem 
description  are  as  follows: 

Jet  location  (single  jet  issuing  from  the  fuselage): 

Xj  =  208  inches 

V° 


Zj  =  -30.8  inches 

Jet  diameter  =  22.5  inches 


Flight  Conditions: 
or  -  0 


U 

Jet  velocity  ratio  =  0.2  0  -  0 

jo 

All  the  above  dimensions  are  ten  times  wind  tunnel  model  dimensions,  as  will  be 
the  case  throughout  this  section. 

Jet  inclination  angles: 

The  jet  will  be  taken  to  be  exhausting  along  the  negative  Z  direction  for  the 
sample  problem. 


2.  APPLICATION  OF  THE  MAPPING  METHOD 

A  complete  description  of  the  fuselage  will  not  be  given  here  since  this  is  not 
necessary  to  describe  the  application  of  the  mapping  method.  Instead,  a  complete 
treatment  of  one  section  of  the  body  will  be  given,  this  being  sufficient  to  demonstrate 
application  of  the  method. 

Figure  29  shows  the  section  of  the  wind  cunnel  test  model  body  at  station  264.  25 
together  with  terminology  for  mapping  into  a  circle.  This  section  has  been  rotated 
90°  counterclockwise  so  that  the  axis  of  symmetry  is  along  the  X-axis  with  the  bottom 
of  the  section  cutting  the  positive  X-axis.  This  coordinate  system  is  not  related  to  the 
original  fuselage  coordinate  system  but  is  in  the  terminology  used  in  the  mapping 
program.  This  location  of  the  section  is  the  proper  one  for  the  mapping  method,  and 
the  mapping  coefficients  obtained  with  this  orientation  are  in  the  correct  form  for  the 
transformation  method.  Figure  30  shows  the  inputs  to  the  mapping  computer  program 
required  to  obtain  a  mapping  of  this  section. 

The  first  card  contains  four  integers.  The  first  of  these  number;  specifies  the 
number  of  points  specified  about  the  section.  The  second  number  b,  eciiies  the  number 
of  corners  around  the  section  (for  a  symmetrical  section  this  would  represent  the  num¬ 
ber  of  pair3  of  corners  except  for  corners  on  the  X-axis),  which  for  this  section  is 
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FIGURE  30.  STATION  264.  25  OF  TEST  MODEL  FUSELAGE 
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zero.  The  third  number  represents  the  number  of  terms  in  the  expansions  for  the 
potential  and  for  the  mapping  function.  The  last  number  on  this  card  specifies  that 
the  section  being  mapped  is  symmetrical. 

Cards  2  through  6  specify  the  X-coordinates  of  the  (in  this  case)  37  points  being 
taken  about  the  section  starting  on  the  positive  X-axis  and  ending  on  the  negative  X-axis. 
Cards  7  through  11  specify  the  Y-coordinates  of  these  same  points  on  the  section. 

Since  there  are  no  comers  to  be  specified  on  this  section,  the  next  number  (on 
card  12)  specifies  what  shift  is  desired  to  translate  the  body  along  the  X-axis.  In  this 
case  there  is  no  need  to  shift  the  section  location  as  it  is  sufficiently  well  centered. 

Card  13  specifies  inputs  needed  to  specify  parameters  for  the  numerical  inte¬ 
gration  of  the  mapping  function  obtained.  The  first  two  numbers  specify  the  X-  and 
Y-  locations  of  the  initial  point  of  the  mapping. 

The  next  three  numbers  specify  the  angular  range  (about  the  mapping  circle)  of 
points  to  be  obtained  on  the  section,  and  the  approximate  spacings  to  be  obtained.  In 
this  case  it  is  specified  that  points  from  0°  to  180°  around  the  mapping  circle  are  to 
be  calculated  and  are  at  an  interval  size  of  5°. 

Card  14  specifies  a  similar  set  of  parameters  for  the  analytically  integrated 
mapping.  It  specifies  that  37  points  are  to  be  obtained  with  a  spacing  of  5°  of  theta 
and  that  the  points  are  to  start  at  0  =  0°. 

Figure  31  shows  the  output  of  the  mapping  program.  These  outputs  have  been 
described  in  Section  II.2  and  so  will  not  be  further  explained  here. 

Figure  32  shows  a  plot  of  the  mapping  as  obtained  by  analytical  integration.  A 
comparison  of  Figure  32  and  Figure  29  shows  that  the  mapped  section  must  be  shifted 
4.42  inches  in  the  negative  X  direction  to  give  the  best  fit.  This,  then,  is  the  value  of 
the  constant  term  in  the  mapping. 

The  complete  tabulation  of  coefficients  for  the  wind  tunnel  test  model  body  with 
canopy  off  are  shown  in  Table  n.  The  values  of  drc/dx,  the  rate  of  change  cf  the  map¬ 
ping  circle  radius  with  body  station,  have  been  obtained  by  graphical  differentiation  of 
rc  plotted  versus  x.  In  addition  to  the  coefficients  tabulated  in  Table  n,  the  initial 
mapping  coefficient  a_p  required  as  input  to  the  jet  program,  is  equal  to  1.0  for  all 
the  fuselage  stations. 
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FIGURE  32.  OUTPUT  OF  MAPPING  FUNCTION  PROGRAM  AT  TEST  FUSELAGE  STATION  264. 
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FIGURE  32.  (Continued) 


MAPPING  OF  SECTION  WITH  CORNERS  REMOVED 
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FIGURE  32.  (Concluded) 


TABLE  n.  COEFFICIENTS  OF  MAPPING  FOR  FUSELAGE 


3.  APPLICATION  OF  JET  FLOW  FIELD  THEORY  TO  FUSELAGE 

The  purpose  of  the  Jet  Flow  Field  theory,  when  used  in  conjunction  with  the 
Transformation  Method,  is  to  predict  jet-induced  velocity  components  at  the  control 
points  on  the  fuselage  required  by  the  Transformation  Method  to  evaluate  power 
effects.  This  is  accomplished  by  executing  the  Jet  Flow  Field  program  to  generate 
required  data  for  the  Transformation  Method  in  the  form  of  punched  data  cards.  To 
insure  compatibility  with  the  Transformation  Method,  the  control  points  on  the  fuse¬ 
lage  where  induced  velocity  components  are  to  be  computed  are  specified  by  utilizing 
the  mapping  coefficients  for  the  fuselage  cross  sections  obtained  in  Section  III.2.  The 
punched  output  is  generated  in  a  manner  providing  a  continuous  block  of  input  data  to 
the  Transformation  Method  computer  program.  Both  of  the  above  points  will  be 
described  in  greater  detail  in  the  discussion  of  the  sample  problem  computations. 

It  should  also  be  noted  that  the  application  of  the  Jet  Flow  Field  program  to  pro¬ 
vide  data  to  the  Transformation  Method  program  for  the  computation  of  power  effects 
on  the  fuselage,  differs  only  slightly  from  its  application  to  computing  power  effects 
on  the  wing,  discussed  in  Section  II.3.  Consequently,  much  of  the  discussion  below 
will  parallel  that  of  Section  If.3. 


a.  Sample  Problem  Computation 

For  the  sample  problem  being  considered,  the  Jet  Flow  Field  program  is  now 
used  to  compute  the  jet-induced  velocities  at  the  16  fuselage  stations  described  in 
Section  A  sketch  of  the  fuselage  and  th^  location  of  the  jet  with  respect  to  the 

input/output  coordinate  system  is  shown  belcw.  The  jet  exhaust  angles  0  and  are 
also  defined. 
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(1)  Input  for  Sample  Problem 

The  input  cards  required  for  the  sample  problem  are  tabulated  in  Figure  34. 

Card  1  lists  three  control  indices.  The  first  one,  MULT  =  1,  indicates  that  a 
single  jet  configuration  is  being  treated.  The  second  one,  IGE0M  =  2,  specifies  that 
control  points  on  fuselage  cross  sections  will  be  generated  utilizing  the  mapping  coef¬ 
ficients  obtained.  The  third  control  .ndex,  IPUNCH  =  1,  generates  the  punched  output 
for  the  Transformation  Method  program. 

Card  2  specifies  angle  of  attack,  a  -  0  and  angle  of  sideslip,  0=0. 

Card  3  controls  the  number  of  steps  and  the  step  size  in  the  numerical  integra¬ 
tion  of  the  equations  of  motion  for  the  jet  path.  For  the  sample  problem,  90  steps  with 
a  step  size  of  .4  (jet  exit  diameters)  are  chosen. 

Cards  4  and  5  contain  information  about  the  jet.  The  jet  location,  in  the  coordi¬ 
nate  system  of  Figure  33,  is  X  =  208.,  Y  =  0,  Z  =  -30.8.  The  jet  exhaust  angles  <t>  and 
\l> ,  defined  in  Figure  33(b),  are  180  and  0  degrees,  respectively.  The  jet  exit  diame¬ 
ter,  do  =  22.5  and  the  velocity  ratio,  U^/Ujo  =  .2. 

Card  6  may  be  left  blank  for  single  jet  computations.  For  a  multiple -jet 
configuration  it  would  be  used  to  control  the  geometry  of  the  jet  resulting  from  the 
intersection  of  two  other  jets. 

The  remaining  input  cards  provide  data  to  generate  the  control  points  at  which 
jet-induced  velocities  are  to  be  evaluated.  These  control  points,  in  order  to  insure 
compatibility  v/ith  the  Transformation  Method,  are  generated  by  utilizing  the  mapping 
coefficients  and  mapping  circle  radii  obtained  for  the  16  fuselage  stations  of  the  sam¬ 
ple  problem.  The  number  of  control  points  generated  at  each  fuselage  station  is  gov¬ 
erned  by  the  input  variable,  NTKT,  which  is  the  number  of  increments  A0  into  which 
the  mapping  circle  (or  mapping  semicircle  if  only  half  the  fuselage  is  to  be  mapped  at 
each  station)  is  to  be  divided  in  the  Transformation  Method  computer  orogram.  Since 
the  flow  is  symmetric  (0  -  0),  computations  will  be  carried  out  for  only  half  the  body. 
Consequently,  the  number  of  control  paints  at  each  fuselage  station  will  be  NTHT  +  1. 

Card  7  specifies  that  NTHT,  the  number  of  equal  increments  A0  into  which  the 
mapping  semicircle  is  divided,  is  18,  which  will  generate  19  control  points  at  each 
fuselage  station.  It  defines  the  number  of  fuselage  stations  NS  -  16.  It  also  defines 
the  number  of  terms  used  in  the  mapping  expansion,  NCOEF  -  11,  and  through  the 
control  index  NSYM  -  0  indicates  that  computations  are  to  be  effected  for  only  half 
the  iuselagc  at  each  station. 
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FIGURE  34.  JET  FLOW  FIELD  PROGRAM  INPUT  DATA 
FOR  SAMPLE  PROBLEM  (Fuselage) 


The  data  in  cards  8-55  are  taken  from  Table  E,  which  is  determined  in  Section 

ffi.2. 

Cards  f1— 10  provide  the  data  from  which  the  fuselage  c.xis :  section  at  the  first 
station  can  be  generated  by  the  computer  program.  Card  8  specii.es  tie  location  of 
the  station,  X  =  23.7,  the  mapping  circle  radius  R  =  10.2633,  and  the  rate  of  change  of 
R  with  X,  DRDX  =  .285. 

Cards  9  and  10  list  the  real  parts  of  the  coefficients  to  be  used  in  the  mapping 
i  expansion  (thus  only  symmetrical  fuselages  may  be  treated). 

Cards  11-55  are  similar  data  blocks  for  fuselage  stations  X  =  41.,  73.,  94.,  118., 
143.5,  162.5,  185.5,  221.5,  264.5,  316.,  343.,  374.,  411.,  450.,  and  497. 

Note:  The  rate  of  change  of  the  mapping  circle  radius  with  distance  along  the 
fuselage,  DRDX,  is  not  required  for  any  of  the  computations  performed  by  the  Jet 
Flow  Field  program.  It  will,  however,  be  required  by  the  Transformation  Method 
program,  and  is  read  as  part  of  the  input  so  that  it  may  be  punched  out  in  the  proper 
sequence  in  the  data  package  to  be  provided  to  the  Transformation  Method  program. 

(2)  Output  for  Sample  Problem 

For  the  sample  problem  being  considered,  both  printed  and  punched  outputs  are 
obtained. 

•  Printed  Output 

Figure  35(a)  shows  the  first  page  of  printout.  The  jet  configuration  being  treated 
is  identified  both  by  appropriate  heading  and  by  printout  of  pertinent  input  infor¬ 
mation.  Input  controlling  the  numerical  integration  procedure  is  also  displayed. 
Figure  35(b)  shows  a  partial  printout  of  computed  jet  centerline  information. 

The  coordinates  of  the  jet  centerline  are  identified.  The  nondimensionalized  jet 
speed  Uj/UjQ  and  the  nondimensionalized  major  diameter  of  the  ellipse  repre¬ 
senting  the  cross  section  of  the  jet,  d/do,  are  also  printed  out.  These  proper¬ 
ties  are  printed  out  at  each  integration  interval  as  specified  on  card  3  of  the 
input.  The  variables  XCOORD  and  DIA  show  a  monotonic  increase  over  this 
region,  while  UJ  -  Uj/U^o  -  ,?  once  the  jet  speed  Uj  approaches  the  freestream 
speed  U^.  Figure  25(c)  shows  the  printout  for  the  jet-induced  velocity  compo¬ 
nents  at  the  first  fuselage  station  specified,  X  -  23.7.  The  coordinates  of  the 
19  control  points  at  the  station  are  identified.  For  this  symmetric  flow  sample 
pro  ilem,  only  the  positive  half  of  the  fuselage  station  is  generated.  The  induced 
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•••  SIMPLE  JET  CONFIGURATION  ••• 


UET  VJET  Z  JET  PHI  PS  I 

2QN.0000  0.0  -30.0000  100.0000  0.0 


ANCLE  OF  ATTACK  *  0.0 
ANCLE  OF  SIOESLIi'  *  0.0 


NUNOEft  OF  STEPS  IN  INTEGRATION  *  90 

INTEGRATION  INTERVAL  *  0.40  JET  EXIT  01 AHETERS 


FIGURE  35(a).  INPUT  PARAMETERS  FOR  SAMPLE  PROBLEM 


••  SINGLE  JET  CENTERLINE  ** 
*••**•*••••*•••****••**••*••«*•* 


XCOORO 

YCOORD 

Z COORD 

UJ 

DIA 

200.00 

0.0 

-30.8  0 

1.000 

1.00 

208.10 

0.0 

-39.80 

0.948 

1.18 

208.44 

0.0 

-48.80 

0.893 

1.45 

209.10 

0.0 

-57.80 

0.833 

1.90 

210.26 

0.0 

-66.80 

0.760 

2.64 

212.10 

0.0 

-75.80 

0.688 

2.93 

214.71 

0.0 

-84.  80 

0.626 

3.23 

218.18 

0.0 

-93.80 

0.573 

3.55 

222.60 

0.0 

-102.80 

0.528 

3.88 

228.11 

0.0 

-111.80 

0.489 

4.23 

234.83 

0.0 

-  120.80 

0.456 

4.  60 

242.92 

0.0 

-129.80 

0.427 

4.98 

252.59 

0.0 

-138.80 

0.402 

5.37 

264.06 

0.0 

-147.80 

0.381 

5.7  7 

277.60 

0.0 

-156.80 

0.363 

6.18 

293.52 

0.0 

-165.80 

0.347 

6.60 

312.22 

0.0 

-174.80 

0.333 

7.  03 

334.13 

0.0 

-183.80 

0.321 

7.47 

359.77 

0.0 

-192.80 

0.311 

7.92 

389.78 

0.0 

-201 '80 

0.301 

8.38 

424.86 

0.0 

-210.80 

0.293 

8.86 

465.87 

0.0 

-219.80 

0.285 

9.34 

513.81 

0.0 

-228.80 

0.278 

9.85 

569.84 

0.0 

-237.80 

0.272 

10.37 

635.36 

0.0 

-246.80 

0.266 

10.91 

711.98 

0.0 

-255.80 

0.261 

11.47 

FIGURE  35(b).  JET  CENTERLINE  FOR  SAMPLE  PROBLEM 


•*»  IXH-CEQ  VELOCITIES  ON  BOOT  ••• 


••••••••••••••• •••**•* 


X 

Y 

2 

U 

V 

V 

23.700 

0.0 

-18.691 

-0.4473  IE-04 

0.0 

0.93589E-02 

23.700 

1.731 

-18.452 

—0.2 9299 £-06 

0.74342E-04 

0.93522E-02 

23.700 

3.392 

-18.003 

0.17127E-04 

0. 14471E— 03 

0.93320E-02 

23.700 

4.929 

-17.252 

0.93755E-04 

0.20799E-03 

0.92982E— 02 

23.700 

6.300 

-16.217 

0.19770E-03 

0.26181E-03 

0.92517E—  02 

23.700 

7.462 

-14.928 

0. 32468E-93 

0 . 30423E-O3 

0.91937E-02 

23. 700 

8.375 

-13.417 

0.47016E  -03 

0.33399E— 03 

0.51253E-02 

23.700 

9.013 

-11.720 

0.62943E-03 

0. 35036E— 03 

C.90481E-02 

23. 700 

9.359 

-9.878 

C.79768E-03 

0.35398E— 03 

0.89638E-02 

23.700 

9.414 

-7.  930 

0.97076E— 03 

9. 34585E— 03 

0.88739E-02 

23.700 

9.211 

-5.925 

O.II440E-02 

0.32832E-03 

0.87802E-C2 

23.700 

8.797 
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FIGURE  35(C)..  INDUCED  VELOCITY  COMPONENTS  AT  STATIONS  23,7  AND  437. 
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velocity  components  (J,  v,  W,  all  nondimensionalized  fay  U^,  are  printed  out  for 
each  centred  point.  Figure  35(ci  also  shows  the  printout  for  the  last  fuselage 
station  considered  in  this  problem,  X  -  497.  Similar  printouts  are  obtained  for 
the  other  intermediate  stations  specified  as  part  of  the  input. 

•  Punched  Output 

The  punched  output  for  the  sample  problem  is  shown  in  tabulated  form  in  Figure 
36.  The  output  data  block  for  the  first  fuselage  station  is  identified.  The  first 
card  lists  foe  fuselage  station  X  =  23.7,  the  mapping  radius  R  =  10.2633  and  the 
rate  of  change  of  R  with  X,  DRDX  =  .265.  The  next  two  cards  list  the  real  parts 
of  the  coefficients  used  in  the  mapping  expansion.  Cards  4—7  list  the  induced 
velocity  components  in  the  X-direction  for  each  of  the  19  control  points  at  fuse¬ 
lage  station  X  =  23.7.  The  induced  velocity  components  in  the  Y-direction  are 
listed  on  cards  8—11  and  cards  12—15  specify  the  induced  velocity  components 
in  file  Z-direction.  Data  blocks  of  this  type,  each  consisting  of  15  cards,  follow 
for  each  of  the  other  15  fuselage  stations  specified  as  part  of  the  input.  The 
punched  output  is  identified  in  columns  73—80.  The  fuselage  station  number  is 
shown  in  columns  75-77.  Sequence  numbers  for  each  station  appear  in  columns 
78-80.  The  letters  U,  V,  W  and  column  74  identify  the  velocity  components 
listed  on  the  data  cards. 

Note:  From  the  tabulations  of  Figure  36  it  is  apparent  that  the  first  three  data 
cards  of  'he  data  generated  for  each  fuselage  station  represent  an  exact  duplication  of 
input  cards  described  previously.  They  are  generated  as  part  of  the  punched  output 
so  that  a  more  complete  data  block  for  the  Transformation  Method  program  may  be 
obtained  without  additional  card  handling. 

b.  Applicability  and  Limitations 


See  discussion  on  applicability  and  limitations  in  Section  II.  3. 
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FIGURE  36.  (Continued) 


374.000000 
0.100006  01-0. 
0.150036  07  0. 
-0.7*79 IE -02-0. 
-0. 32 79 7E -02-0. 
0.193236-02  0. 
0.206066-02 
0.0  -0. 
-0.121 1*6-01-0. 
-0.  760*36-"*?-  >. 
0.693056-07 
-0.600676-01-0. 
-0.502616-01-0. 
-0.629976-01-0. 
-0.395076-01 
611.000000 
0.100066  01-0. 
0.631076  06  0. 
-0.616026-02-0. 
-0.380666-02-0. 
-0.370286—03-0. 
0.601276-03 
0.0  -0. 
-0.812126-02-0. 
-0.576666-C2-0. 
0.615586-07 
-0.553356-01-0. 
-0.692026-01-0. 
-0.391096-01-0. 
-0. 366696-01 
650.000 COO 
0.10000E  01-0. 
0.792696  05  O. 
-0.506656-02-0. 
-0. 305696-02-0. 
-0.170506-02-0. 
-0.117606-02 
0.0  -0. 
-0.660706-02-0. 
-0.363206-02-0. 
0. 2603C6— 07 
-0.6536 5E -01-0. 
-0.618806-01-0. 
-0.358166-01-0. 
-0.339626-01 
697.000000 
0.10000E  01-0. 
0.76206E  03  0. 
-0. 3866 9E -02-0. 
-0. 36666E-02— 0. 
-0.26666E-02-0. 
-0.236006-02 
0.0  -0. 
-0. 13892E-02-0. 
-0.130166-02-0. 
0.91033E-08 
-0.36283E-01-0. 
-0.35053E-01-0. 
-0.32593E-01-0. 
-0. 3165CE-01 


27.511200  -0.076000 

110766  02  0.67503E  02  0.273316  03-0.267586  05-0.130196  06 
57690E  07-0.937206  09  C. 133666  10  0.237966  12 
733636-02-0.692266-02-0.626616-02-0.539936-02-0.637706- 
21 715E-02-0.10C68E-C2-0. 750366-06  0.787176-03  0.165116- 
227906-02  0.253026-02  0.269806-02  0.279966-02  G.285C8E- 


356 75E-02-0.673166-0 2-0.926906-02-0. 109536-01-0. 118736- 
110066-01-0. 1111 16-01-0.102716—01 -0.9665 i 6-02-0. 861 556- 
661306— 02 -0.5*3676-0 2-0. *18506— 07— 0.28538E— 02-0. 1*4786- 

606936-01-0.376666-01-0.658366-01-0.636656-01-0.610776- 
55296F-Ci-0„52301E-0 1-0.69376E-01-0.66 768E-01-0.666 166- 
61 7786*01-0. *08676-0 1-0. *02356— 01-0.398506-01-0.396*96- 


02  U 
02  U 
02  U 
U 

V 

V 

V 

V 
M 
M 
W 
M 


01 

02 

02 

01 

01 

01 


23.677185  -0.135000 

166396  02  0.25525E  02  0.133966  03-0.161226  05-0.692196  05 
375996  07-C.26662E  09-0.960706  96  0.31622E  11 
60936F-02-0.58935E-02-0.55605E— 02— 0.5099 36-0 2-0. 653016- 
31 9996-02 -0.26978E-02-0. 182186-02-0. 1 22606-02-0 . 7*55*6- 
100296-03  3.107276-03  0.251956-03  0.361196-03  0.387286- 


22508E-02-0.*2863E-02-0.59365E-02-0.711*1E-02-0.783C1E- 
805956-02-0. 77556E-0 2-0. 73*176-02-0.6892 16-02-0.638616- 
501526-02— 0. 61590E-02-0. 3222 16— 02-0. 22086E-02-0  112506- 

5515*E-01-0.5*6liE-01-0.537i36-01-0.52*7*E-01-G.509**E- 
47322E-01-0. *536*6-0 1-0. *3* 34E-01-0.41685E-0 1-0.402 386- 
38261E-01-0. 37585E-01-0. 37126E-01-0.36839E-0 1-0.366936- 


02  U 
03  U 
03  U 
U 

02 
02 
02 


01 

01 

01 


17.579895  -0.185000 

176556  02  0.266976  02  0.802066  02-0.306986  06-0.1556B6  05 
*29786  06-0.207396  08-3.562666  08  0.229706  10 
SO 323E-02-0.49296E-0 2-0. *75636-02-0.45 12 36-02-0. *206 1£- 
3*7396-02-0. 3075  76-02-0. 268246-02-0.23256E-02-0. 202686- 
1S886E-02-0. 1*3196-02-0. 131506-02-0. 12358E-02-0. 118956- 


1 1 3996-02-0.218796-02-0.3068*6-02-0. 3  7*386-02-0. *21336- 
45798E— 02— 0.4533 1E-02—0.44080E-02— 0.423* IE-02-0. 398 726- 
3172*6-02-0.263*26-02-0.203736-02-0.138886-02-0.702776- 

*52*7E-01-0.**951E-01-0.***55E-31-0.*3T586-01-0.*288*E- 
40786E-01— 0.396 35E-0 1-0. 38486E-0 1—0. 37*3 2E-0 1-0. 365*06- 
352216— 01-0. 34746E-0 1—0. 3*3916-01-0.3*1506-01-0.340096- 


02  U 
02  U 
02  U 
U 

V 

V 

V 

V 
w 
w 
u 

M 


02 

02 

03 

01 

01 

01 


8.212999  -0.206000 

212526  02  0.131166  02  0.80*016  01-0.11905E  03-0.329606  03 
16588E  04-0.34593E  05  0.90360E  03  0.11027E  07 
38559E-02 -0.382 29 E-0 2-0. 37675E-02-0. 36882E-02-0. 358586- 
33288E— 07-0. 31835E-02-0. 30356E-02-0.28954E-02-0.27713E- 
25732E-02 -0. 24969E-02-0. 2*3716-02-0.239*56-020.236876- 


327776-03-0.63*626-03-0.900606-03-0.11 1**6-02-0.127626- 
1*5526-02-0. 1*81 5E-0 2-0. 1*81 3E-02-0. 1*5926-02-0. 14038E- 
11519E-02-0.96575E-03-0. 751*96-03-0. 51400E-0 3-0. 26051E- 

362*96-01-0.361*86-0 1-0. 35979E-01-0. 357366-01-0.35*2*6- 
3*6386-01-0.341936-0 1-0. 3 3738F-0 1-0. 33  3066-01-0 . 329236- 
32 3106-01 -0. 3207*6-0 1-0.31889E-01-0. 31 757E-01-0. 316776- 


02  U 
02  U 
02  U 
U 

V 

V 

V 

V 

w 
w 
w 
w 


02 

02 

03 

01 

01 

01 


13 

13 

13 

13 

13 

13 

13 

13 

13 

13 

13 

13 

13 

13 

13 
1* 
1* 
1* 
1* 
1* 

14 
1* 
1* 
1* 
1* 
1* 
1* 
1* 
!* 
1* 
>5 

15 
15 
15 
15 
15 
13 
15 
15 
15 
15 
15 
15 
15 

15 

16 
16 
16 
16 
16 
16 
16 
16 
16 
lb 
16 
16 
16 
16 
16 


1 

2 

1 

2 
3 
* 
1 
2 

3 

4 
1 
2 
3 
* 

1 

2 

1 

2 

3 

* 

1 

2 

3 

* 

1 

2 
3 
* 

1 

2 

1 

2 

3 

<> 

1 

2 
3 
* 
1 
2 
3 
* 

1 

2 

1 

2 

3 

* 

1 

2 

3 

* 

1 

2 

3 

* 


FIGURE  36..  (Concluded) 


4.  APPLICATION  OF  TRANSFORMATION  METHOD  TO  FUSELAGE 


This  section  bears  close  resemblance  to  the  d  icussion  in  Section  II.  4  of  "Appli¬ 
cation  of  Transformation  Method  to  Wing."  A  finite  wing  is  also  &  three-dimensional 
body  bat  it  has  a  distinctive  property  of  generating  lift  in  a  uniform  stream.  Thus,  by 
ignoring  a  part  of  the  computer  program  which  is  concerned  with  circulation,  the  power 
effects  on  a  fuselage  are  calculated  in  a  manner  similar  to  that  for  a  wing. 

a.  Inputs  to  Transformation  Method  for  Sample  Problem 

All  the  input  data  for  this  sample  problem  are  shown  in  Figure  37.  The  punched 
outputs  from  the  jet  flow  field  program  again  make  up  the  major  portion  of  the  input. 

Two  cards  must  also  precede  this  basic  Input  block  to  specify  various  flow  indices, 
and  none,  two  or  four  cards  may  follow  this  block,  depending  on  the  options. 

Card  1  and  Card  2  denote  the  same  number  of  flow  indices  as  those  on  Cards  1 
and  2  in  the  previous  section  on  application  to  a  wing.  The  meaning  of  each  input  is 
also  the  same,  except  that  the  /turner  J  quantities  of  the  following  indices  are  differ¬ 
ent.  The  classification  index  (IGECM)  is  now  equal  to  2  to  denote  a  fuselage,  the 
number  of  iterations  (JSTOP)  i*  1,  the  number  of  stations  ,'NSTA)  is  16,  the  computation 
index  (NSYM)  is  0  to  indicate  the  existence  of  a  plane  of  symmetry,  and  ithe  number 
of  angular  increments  from  0  to  *  (MTHET)  is  18. 

Cards  3  through  242  contain  the  punched  output  data  furnished  by  the  jet  flow 
field  program,  which  include  the  X  coordinate,  the  mapping  coefficients,  and  the 
induced  velocity  component  for  stations  1  through  16.  There  are  15  cards  for  every 
station. 

Card  243  specifies  th*  station  number  immediately  preceding  the  exhausting  jet. 

Card  244  refers  to  the  X  coordinate  of  the  fuselage  tail. 

Card  245  lists  in  order  the  number  of  jets,  the  jet  exit  diameter,  X  coordinate 
of  the  center  of  gravity,  and  the  reference  length  for  nondimensionalizing  computed 
moments. 

Card  246  denotes  in  order  the  Y  coordinate  of  the  nose,  Z  coordinate  of  the  nose; 

X  coordinate,  Y  coordinate  and  Z  coordinate  of  the  tail. 

u.  Outputs  fr  jm  Transformation  Method  for  Sample  Problem 

Figure  38  lists  diree'  ■  ir  indirectly  a  portion  of  the  input  uata  on  cards  1 
through  242. 
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Figure  39  establishes  the  correspondence  between  the  angular  increments  of  the 
mapping  circle  and  their  corresponding  geometric  locations  at  every  fuselage  station. 
The  first  column  gives  the  angular  increments  in  degrees. 

Figure  40  shows  the  pressure  coefficient  distribution  at  each  station  after  appli¬ 
cation  of  the  segment  method.  Since  these  coefficients  are  tabulated  against  angular 
increments  only,  cross  reference  to  Figure  j.s  needed  in  order  to  establish  pres¬ 
sure  distributions  in  the  physical  plane. 

Figure  41  is  also  a  table  for  the  pressure  distribution  similar  to  the  one  in  Fig¬ 
ure  40  but  after  the  three-dimensional  modification  has  been  introduced. 

Figure  42  gives  the  parameters  used  in  the  three-dimensional  modification  and  in 
the  computation  of  forces  and  moments,  originally  read  in  as  input  data  on  Cards  243 
through  246.  The  computed  forces  (normalized  to  the  thrust)  and  moments  (normalized 
by  the  thrust  and  reference  length)  on  the  fuselage  after  one  iteration  are  also  tabulated 
in  the  same  figure. 

Some  of  the  computed  results  of  this  sample  problem  have  been  compared  with 
the  available  wind  tunnel  test  data  and  are  shown  in  Figures  43(a)  through  43(d).  The 
power-effect  pressure  distributions  in  these  figures  are  representative  for  fuselages 
with  a  lift  jet  in  a  uniform  flow.  The  typical  feature  of  this  class  may  be  described 
as  follows:  The  power  effect  in  the  fore  part  of  the  fuselage  is  generally  very  small, 
as  seen  in  Figure  43(a).  As  we  approach  the  lift  jet,  there  is  a  small  region  in  the 
front  where  the  power  effect  is  positive  (Figure  43(b)).  The  negative  power  effect  then 
decreases  rapidly  and  its  magnitude  becomes  very  large  in  the  immediate  neighborhood 
of  the  jet  but  tapers  off  to  a  constant  level  in  the  distance  of  one  or  two  jet  diameters. 
This  constant  level  of  power  effect  is  in  general  quite  high  and  prevails  over  the  entire 
back  part  of  the  fuselage  (Figures  43(c)  and  43(d)).  Thus,  it  contributes  a  n  aior  por¬ 
tion  of  the  forces  and  moments.  The  origin  of  this  effect  is  due  mostly  to  the  wake 
formed  behind  the  lift  jet.  Its  prediction  lies  beyond  the  scope  of  the  preserl  method. 
Consequently,  the  discrepancy  between  the  calculation  and  the  test  data  in  this  region 
is  large. 

When  the  sideslip  angle  is  not  zero,  the  jet  wake  region  does  not  completely 
enfold  the  back  portion  of  the  fuselage.  Ihe  comparison  between  the  prediction  and 
the  test  becomes  more  favorable  as  may  be  seen  in  Figures  44(a)  to  44(c). 

Further  calculations  and  wind  tunnel  test  data  are  compared  in  Figures  45  and 
46.  The  body  referred  to  in  Figure  46  is  the  one  tested  at  Northrop  prior  to  the 
present  contract  (see  Figure  23  in  Section  II.  4). 
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FIGURE  38.  PARTIAL  LIST  OF  INPUT  DATA  FOR  SAMPLE  FUSELAGE 
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FIGURE  39.  (Continued) 
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FIGURE  39.  (Continued) 
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FIGURE  39.  (Concluded) 
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FIGURE  40.  POWER-EFFECT  PRESSURE  COEFFICIENTS  ON  SAMPLE  FUSEL  AO  E 
AFTER  APPLICATION  OF  SEGMENT  METHOD 
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FIGURE  40.  (Continued) 
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FIGURE  40.  (Concluded) 


PRESSURE  COEFFICIENTS  AT  FUSELAGE.  THREE  01  HENS IONAL  MODIFICATION  OF  1  ITERATION 
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FIGURE  41.  POWER-EFFECT  PRESSURE  COEFFICIENTS  ON 
SAMPLE  FUSELAGE  AFTER  ONE  ITERATION 
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FIGURE  41.  (Continued) 


PRESSURE  COEFFICIENTS  AT  FUSELAGE,  THREE  DIMENSIONAL  MODIFICATION  OF  1  ITERATION 


FIGURE  41.  (Concluded) 


FIGURE  42.  FORCES  AND  MOMENTS  ON  SAMPLE  FUSELAGE 
BY  TRANSFORMATION  METHOD 


FIGURE  43a  POWER-EFFECT  PRESSURE  COEFFICIENTS  ON  SAMPLE 
FUSELAGE  AT  STATIONS  X  =  73.0  AND  X  =  118.0 
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FUSELAGE  AT  STATION  X  =  185.  5 


Test  D 


L?ft  Jet 


FIGURE  43d.  POWER-EFFECT  PRESSURE  COEFFICIENTS  ON  SAMPLE 
FUSELAGE  AT  STATION  X  =  343.0 


Test  Da 


FIGURE  44a.  POWER-EFFECT  PRESSURE  COEFFICIENTS  ON  FUSELAGE 

AT  STATIONS  X  =  73.  0  AND  X  =  118. 0 
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AT  STATION  X  =  185.  5 


Segment  Method 
3-D  Modification 


FIGURE  44c.  POWER-EFFECT  PRESSURE  COEFFICIENTS  ON  FUSELAGE 
AT  STATIONS  X  =  264.  25  AND  X  =  343.0 


A  Test  Data 

_  O  Calculation  (Transformation) 

□  Calculation  +  Wake  Correction  (Appendix  IV) 


c.  Method  Applicability  and  Limitations 


As  pointed  oat  in  the  beginning  of  this  section,  the  method  for  predicting  power 
effect  on  a  fuselage  is  very  similar  to  that  of  a  wing.  However,  for  the  sample 
problems  considered  in  this  Volume,  there  are  two  differences  between  fuselage  and 
wing  computations.  The  first  one  is  the  formation  of  a  wake  behind  the  exhausting  jet 
which  engulfs  a  large  portion  of  the  fuselage  and  is  a  fundamental  problem.  The 
second  difference  is  computational,  and  may  be  briefly  described  as  follows:  The 
induced  velocities  usually  undergo  very  large  changes  across  the  jet,  as  discussed 
previously.  These  are  most  noticeable  in  the  mainstream  direction,  which  is  also 
the  longitudinal  direction  of  a  fuselage  without  sideslip.  When  we  invoke  numerical 
differentiations  in  the  longitudinal  direction  to  determine  the  strength  of  the  residual 
sources  and  sinks,  these  large  variations  generally  magnify  the  intensity  of  the 
sources  and  sinks  and  may  cause  erratic  behavior  in  the  final  results.  This 
irregularity  usually  gets  worse  when  higher  iterations  are  carried  out,  even  though 
some  smoothing  of  the  input  data  is  applied.  For  this  reason,  two  iterations  were  not 
performed.  To  compare  the  case  of  a  wing,  wa  notice  that  the  largest  variation  in  the 
induced  velocities  exists  in  the  mainstream  direction  and  does  not  traverse  any  of  the 
wing  stations.  In  addition,  the  computational  procedure  automatically  smooths  out 
some  of  the  large  variations  in  the  mainstream  direction  by  integration  for  the 
"boundary  functions"  (see  Volume  I,  Section  IV,  for  details)  and  the  subsequent 
expansion  in  a  Fourier  series. 

Since  the  power  effect  on  a  fuselage,  disregarding  the  wake,  contributes  only  a 
fraction  to  the  total  power-effect  forces  and  moments,  some  uncertainty  in  prediction 
may  be  tolerable.  Based  on  this  assumption,  the  segment  method  is  useful  since  it 
includes  some  of  the  three-dimensional  effects  that  are  already  present  in  the  jet  flow 
field  and  gives  fairly  reasonable  results  in  most  cases. 

As  in  the  case  for  a  wing,  the  present  computer  program  is  also  capable  of 
treating  the  power-on  and  power-off  problems  in  a  formal  sense.  However,  the 
three-dimensional  effects  due  to  the  nose  and  the  tail  of  the  fuselage  are  not  included. 

In  the  examples  calculated  so  far,  no  smoothing  procedure  was  used  for  the 
input  data. 
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SECTION  IV 


POWER  EFFECTS  ON  CONTROL  SURFACES 

To  predict  the  aerodynamics  oi  a  V/STOL  aircraft  it  is  necessary  to  be  able  to 
predict  the  effects  of  power  on  the  horizontal  and  vertical  tail  surfaces.  This  power 
effect  can  be  attributed  primarily  to  an  induced  flow  angle  at  the  tail  location.  Having 
a  method  of  predicting  the  power  induced  flow  angle  it  is  then  possible  to  estimate  the 
aerodynamic  forces  and  moments  induced  on  the  tail  surfaces  by  the  exiting  jet.  This 
section  will  describe  bow  to  obtain  the  jet  induced  flow  angles  at  the  location  of  the 
tail  surfaces. 

1.  SAMPLE  PROBLEM 

The  sample  problem  to  demonstrate  foe  calculation  of  foe  downwash  at  the  tail 

will  be  the  wind  tunnel  test  model  described  in  Appendix  I,  with  the  lift  jet  operating. 

The  jet  parameters  and  the  flight  conditions  are  the  same  as  specified  in  Section  III, 

except  that  here  three  values  of  U  /U.  will  be  considered  .1,  .2,  and  .3. 

oo  jo 

Power  induced  downwash  and  sideslip  will  be  computed  at  the  test  rake  location: 
x  =  44.23  inches 

y  =  -6.75  inches  Dimensions  in  model  scale 

z  =  .9,  2.9,  4.9,  6.9,  8.9,  10.9,  12.9  inches 
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2.  APPLICATION  OF  JET  FLOW  FIELD  PROGRAM  TO  CONTROL  SURFACES 


The  use  of  the  jet  flow  field  program  consists  of  specifying  the  points  at  which 
Jet  induced  velocity  components  are  to  be  computed,  details  of  jet  location  and  the 
flight  variables.  A  description  of  sample  input  and  output  for  the  jet  flow  field  program 
has  been  given  in  Section  IL3.  The  application  of  the  jet  flow  field  program  to  tbe  tail 
problem  differs  only  in  tbe  location  of  tbe  control  points. 

In  this  example  the  down-'ash  and  sidewash  angles  will  be  computed  at  the  rake 
location  instead  of  tbe  vertical  or  horizontal  tail  location  since  test  data  exists  for  the 
rake  location. 

Figures  47  and  48  show  a  comparison  of  power  induced  downwash  and  sidewash 
obtained  by  the  jet  program  and  compared  with  wind  tunnel  test  data  for  tbe  body  alone 
with  tbe  rake  attached.  Hie  theory  will  not  account  for  the  presence  of  a  body  or  a 
wing  so  it  must  be  assumed  that  these  components  can  be  ignored  in  calculating  the 
power  induced  angles.  Hie  agreement  shown  in  the  figure  can  be  considered  satis¬ 
factory  considering  the  scatter  in  the  test  data. 

3.  CALCULATION  OF  POWER  INDUCED  FORCES  AND  MOMENTS  ON  CONTROL 
SURFACFS 

To  calculate  the  incremental  force  and  moments  induced  on  the  tail  it  is  necessary 

to  estimate  the  CT  of  the  surface  in  the  presence  of  the  fuselage  (and  other  tail  panels) 
" a 

and  the  cent.oid  of  the  panel  load.  These  values  may  be  estimated  by  empirical  methods 
such  as  are  to  be  found  in  DATCOM.  When  these  values  are  known,  the  jet  induced 
downwash  and  sid'  wash  can  be  used  to  estimate  an  incremental  angle  of  attack  or  side¬ 
slip  on  the  surface  in  question.  From  these  values  it  is  possible  to  estimate  power 
induced  forces  and  moments  on  the  tail  surfaces.  Some  estimates  of  this  nature  have 
been  made  for  the  wind  tunnel  test  model  of  this  study  but  the  accuracy  with  which  the 
incremental  forces  and  moments  can  be  measured  from  the  test  data  precludes  any 
conclusions  as  to  method  accuracy.  The  comparisons  made  do  show,  however,  that 
such  estimates  are  of  the  right  order  of  magnitude. 
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SECTION  V 


APPLICATION  OF  INLET  METHOD 

The  method  of  force  and  moment  estimation  for  normal  inlets  described  in  Volume 
I,  Section  V,  consists  of  three  parts.  The  inlet  induced  forces  are  two  parts,  lip  forces 
which  act  in  the  immediate  vicinity  of  the  inlet,and  surface  forces  which  act  at  greater 
distances.  The  third  part  consists  of  a  description  of  the  net  thrust  of  the  propulsive 
device  causing  the  inlet  flow.  An  approximation  for  a  lift  fan  has  been  derived ,  but 
data  obtained  from  test  could  be  used  in  place  of  the  model. 

The  inlet  method  may  be  applied  to  a  wide  variety  of  configurations.  This  versa¬ 
tility  is  a  direct  result  of  the  empirical  nature  of  the  model.  Because  of 
this  empirical  nature,  some  comment  is  required  on  selection  of  the  parameters  used 
in  the  equations  summarised  below. 
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(  *rinle*+  )  =  /inlet  lip  \  /  surface  \  /netfan\ 

\net  fan  thrust  /  l  forces  /  \ forces  /  (forces/ 


Lip  Forces: 
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Fan  Flow  and  Force: 


&.5 


Summation:  = 

D  ' 

' 

1.  SELECTION  OF  PARAMETERS 

a.  Inlet  Method 

In  die  application  of  the  inlet  method,  judgement  must  be  used  in  determining 
the  boundaries  of  the  surface  upon  which  the  inlet  is  assumed  to  have  an  effect.  The 
magnitude  of  the  pitching  moment  due  to  forward  velocity  as  well  as  the  magnitude 
and  sign  of  the  static  moment  is  a  direct  function  of  the  size  and  distribution  of  the 
effective  area.  In  simple  cases  such  as  the  presence  of  an  inlet  in  the  surface  of  a 
rectangular  wing,  it  is  obvious  that  the  entire  surface  should  be  considered  in  the  cal¬ 
culations.  In  more  complicated  situations  it  is  presently  thought  that  the  area  con¬ 
sidered  for  each  inlet  be  defined  by  a  radius  measured  from  the  centroid  of  the  inlet 
to  the  nearest  edge  of  die  body  or  the  nearest  barrier  to  flow,  such  as  a  sharp  corner. 
Lift  forces  are  concentrated  in  the  immediate  vicinity  of  the  inlet  and  are  not  so  greatly 
affected  by  the  definition  of  the  effective  area. 

The  radius  R^,  used  to  separate  the  area  on  which  "lip"  forces  act  from  the  area 
on  v'hich  "surface"  forces  act,  is  completely  arbitrary.  The  value  of  R^  does  not 
affect  the  final  result  of  the  calculation  which  contains  the  sum  of  "lip"  and  "surface" 

forces.  In  a  particular  calculation  it  may  be  convenient  to  make  R^  correspond  to  the 
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actual  inlet  size  or  to  make  a  larger  value  and  minimize  the  value  of  die  surface 
force  integrals. 

Although  the  integrals  appearing  in  the  surface  force  equations  are  easily 
adapted  to  machine  computation,  .graphical  integration  or  a  finite  summation  using  a 
worksheet  such  as  that  shown  in  Table  3  has  been  found  to  be  satisfactory.  Because 
the  forces  are  concentrated  in  the  immediate  vicinity  of  the  inlet,  large  angular 
increments  may  often  be  used  with  little  loss  in  accuracy.  For  configurations  with 

more  than  one  inlet,  the  present  procedure  is  to  superimpose  the  effects  of  the 

\  ' 

individual  inlets  without  consideration  of  interaction  between  inlets. 

b.  Fan  Flow  Model 

In  application  of  the  fan  flow  model,  values  of  the  parameters  17  and  K  must  be 

2 

selected.  These  parameters  appear  only  in  the  combination  17  -K  and  this  combined 
parameter  should  be  in  the  range  ±1. 

For  deep  inlets  K  may  be  assumed  to  be  zero.  For  very  thin  inlets  the  value  of 
K  will  approach  unity.  Reference  (2)  shows  that  K  =  0. 7  for  a  relatively  deep  fan-in¬ 
wing  installation. 

The  dynamic  head  recovery  factor  17  varies  from  near  unity  for  deep  inlets 
to  zero  for  thin  inlets  Even  for  deep  inlets,  the  value  of  rj  will  fall  quite  rapidly 
with  forward  speed  if  no  flow  turning  device  is  present. 

The  parameter  Ct  is  obtained  by  fitting  static  thrust  versus  RPM  curves  for  a 
particular  fan  after  values  have  been  selected  for  17  -K  and  cDs£g- 

2.  EXAMPLES  OF  APPLICATION 

The  analytical  models  are  applied  to  several  configurations  and  indicate  charac¬ 
teristics  consistent  with  the  data.  It  is  not  possible  to  completely  verify  the  analysis 
as  it  is  not  presently  possible  to  separate  the  effect  of  the  propulsive  wake,  present 
in  the  data,  from  the  effect  of  the  inlet  flow,  and  it  is  not  yet  possible  to  calculate  the 
exit  effect  with  complete  confidence  for  these  configurations.  No  sample  calculations 
are  shown  because  of  the  relative  simplicity  of  the  model,  but  the  selection  of  the 
empirical  parameters  is  indicated. 

Although  a  large  amount  of  data  is  available  for  both  lift-jet  and  lift-fan  con¬ 
figurations,  only  lift  fan  data  is  considered.  Lift -jet  inlet  effects  are  relatively 
small  due  to  the  lower  mass  flows  and  this  may  result  in  inlet  force  variations  of  the 
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same  order  as  data  scatter  ^ .  The  aiuouru  of  usable  data  aval1  able  is  further  reduced 
by  the  presence  of  a  high  degree  of  wind  tunnel  wall  interference  in  some  tests. 

The  methods  have  been  applied  to  four  configurations; 

1.  fan-in-nacelle  ^ 

2.  fan-in-fuselage  ^ 

3.  fan-in-wing  ^ 

4.  1/6-scale  XV5A(7) 

a.  Fan  Flow  Calculations 

Reference  (3)  and  Reference  (4)  have  been  used  to  substantiate  the  propulsion 
model.  The  static  thrust  versus  RPM  curves  were  fitted  for  both  the  fan-in-fuselage 
and  fan-in-wing  configurations.  Using  the  value  of  Ct  so  obtained,  fan  flow  versus 
forward  velocity  was  calculated  and  '<  compared  to  experimental  values  in  Figure  50. 
Both  curves  are  reasonably  well  predicted. 

It  should  be  noted  that  in  matching  the  static  thrust  curve,  the  parameter 
actually  obtained  is  the  combination  (Ct  /l  +  C^S^/A^),  not  Ct .  Therefore,  varia¬ 
tion  of  the  value  of  CDSCB  does  not  affect  the  static  values  produced  by  the  equations 

but  only  varies  the  apparent  value  of  Ct  and  the  effects  of  translational  velocity 

o 

through  the  parameter  combination  ( tj  -K  /I  +  CqS^/A^-).  Thus,  it  may  be  possible 
to  fit  data  while  using  values  for  individual  parameters  which  are  in  error. 

b.  Fan-in-Nacelle 

The  fan-in-nacelle  model  of  Reference  (6)  is  indicated  in  Figure  49.  The  twelve 
inch  fan  was  electrically  driven  and  was  not  equipped  with  inlet  devices  or  exit  vanes 
to  aid  in  turning  the  flow.  The  model  was  reversed  in  the  wind  tunnel  and  the  fan  was 
tested  in  both  leading  and  trailing  positions. 

The  results  yielded  by  the  empirical  model  are  shown  in  Figures  51  and  52;  no 
exit  effects  are  included  in  the  calculated  curves,  which  include  inlet  effect,  net  fan 
thrust,  and  unpowered  aerodynamics.  Available  data  for  similar  inlets  indicate  that 
the  dynamic  head  recovery  factor  can  be  expected  to  decrease  rapidly  as  free  stream 
velocity  increases  from  the  static  conditions.  Therefore,  curves  for  both  total  loss 
and  complete  recovery  are  shown.  The  parameter  K  was  assumed  to  be  zero. 


DIRECTION  OF  FLOW 


x/1  =  0. 64  x/1  =  0. 36 


FIGURE  49.  FAN-IN-NACELLE  CONFIGURATION 

The  method  reflects  the  change  in  flow  direction  as  seen  in  the  lift  data.  Lift 
forces  are  reasonably  well  predicted  by  assuming  high  recovery  at  low  speeds  and  low 
recovery  at  high  speeds.  Drag  is  also  well  predicted.  The  lack  of  accuracy  in  the 
moment  prediction  may  be  due  to  uncertainties  in  the  unpowered  aerodynamics,  exit 
effects  not  included  in  the  predictions,  or  errors  in  the  selection  of  the  propulsion 
parameters.  Increased  fan  flow  rates  would  improve  overall  correlation,  but 
arbitrary  changes  in  propulsion  parameters  would  yield  little  additional  information. 

c.  Fan-in-Fuselage 

The  full-scale  fan-in-fuselage  model  of  Reference  (3)  was  a  shoulder  wing  con¬ 
figuration  of  aspect  ratio  five.  The  single  fan  was  mounted  vertically  in  the  fuselage 
with  the  fan  axis  passing  through  the  wing  quarter  chord  and  moment  center.  A  single 
semicircular  vane  was  placed  behind  the  leading  edge  of  the  inlet  to  improve  pressure 
recovery  and  inhibit  separation.  Only  tail-off,  flap-retracted  data  are  considered. 


13.3 


In  the  application  of  the  empirical  method,  the  depth  of  die  inlet  led  to  the  as¬ 
sumption  that  the  flow  is  axial  to  the  fan;  data  shows  that  pressure  recovery  of  the 

2 

inlet  is  nearly  complete.  Therefore,  the  parameter  tj  -  K  is  assumed  to  be  unity. 

The  diameter  of  the  assumed  circular  inlet  was  chosen  to  contain  the  inter¬ 
section  of  the  actual  inlet  and  the  fuselage.  The  projected  planform  of  the  wing-body 
was  used  in  calculating  induced  surface  forces.  Inlet  sealed  data  was  used  to  repre¬ 
sent  power-off  terms. 

The  results  of  the  calculations  and  data  are  shown  in  Figures  53  and  54.  Figure 
53  compares  fan  flow  rate  and  fan  rotor  thrust.  Figure  54  shows  lift,  drag  and  moment 
coefficients.  It  can  be  seen  that  drag  is  best  predicted  and  pitching  moment  least  well 
predicted.  This  may,  however,  be  due  to  the  lack  of  fan  exit  effects  in  the  analytical 
predictions. 

NASA  TN-D-2560  identifies  the  presence  of  wind  tunnel  wall  interference  in  this 
test  and  noting  that  uncorrected  data  is  presented  here,  removal  of  the  interference 
would  be  expected  to  improve  lift  correlation.  The  effect  on  the  moment  correlation  is 
not  known. 


d.  Fan-in-Wing 

Calculations  were  also  made  for  the  one-sixth  scale  model  of  the  XV-5A  reported 
in  Reference  (7).  The  configuration  of  the  test  model  was  gear  down,  flaps  down,  and 
tail  off.  The  model  moment  center  is  ahead  of  the  fan  axis. 

2 

In  the  application  of  the  model,  the  parameter  tj-K  is  assumed  to  be  zero. 

The  assumed  inlet  size  closely  matches  that  of  the  actual  inlets  which  are  nearly 
circular.  Each  fan  is  assumed  to  produce  induced  forces  only  on  the  wing  panel  in 
which  it  is  mounted.  Data  obtained  with  inlet  and  exit  sealed  are  used  to  predict 
power-off  effects. 

The  results  are  presented  in  Figures  55,  56  and  57.  The  empirical  prediction 
is  presented  both  with  and  without  power-off  aerodynamics.  Fan  exit  effects  present 
in  the  test  are  not  included  in  the  empirical  predictions.  The  forces  and  moments  for 
this  test  were  nondimensionalized  by  the  use  of  the  "slipstream"  dynamic  pressure 
qS,  where: 

qs  p«ui  +  vaf 

L  =  static  lift  at  constant  RPM 
o 

A„  =  total  fan  area 

r 
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Again,  urag  force  is  best  predicted  and  moment  least  well  predicted.  Hie  addition  of 
exit  effects  to  the  prediction  should  improve  the  moment  correlation  and  may  improve 
the  lift  prediction.  Hie  use  of  exit  open,  power-off  data  would  improve  the  tnome.it 
prediction  at  the  expense  of  lift  prediction,  due  to  disturbance  of  flow  on  the  lower 
wing  surface. 
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TABLE  3 

INLET- SURFACE  FORCE  AND  MOMENT  CALCULATIONS 
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FIGURE  50.  FAN  FLOW  PARAMETERS 
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FIGURE  51.  FAN-IN-NACELLE,  INLET  LEADING 
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FIGURE  53.  FAN-IN-FUSELAGE,  PROPULSION  PARAMETERS 


FIGURE  56.  1/6-SCALE  XY-5A,  FAN-IN-WING,  DRAG  FORCE 
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Section  vi 


NONLINEAR  BODY  AERODYNAMICS 

The  body  aerodynamics  are  computed  by  combining  slender  body  theory  with 
viscous  crossflow  effects  to  obtain  nonlinear  coefficients  at  high  angles  of  attack  and 
sideslip.  The  method,  described  in  Volume  I,  has  been  programmed  for  the  computer 
with  the  linear  and  the  viscous  components  of  the  aerodynamic  coefficients  printed  out 
separately.  The  computer  program  requires  that  the  mapping  of  the  body  cross  sec¬ 
tions  be  known.  This  mapping  is  obtainable  by  means  of  the  mapping  method  described 
in  Sections  n  and  III.  A  simplified  method  of  obtaining  inputs  to  the  computer  program 
is  also  described  which  permits  more  ready  calculation  of  body  aerodynamics  by  by- 

,  ,  I 

passing  the  mapping  of  the  body  cross  sections.  1 

1.  SAMPLE  PROBLEM  i 

To  illustrate  the  use  of  the  nonlinear  body  aerodynamics  program  the  wind  tunnel 
test  model  body  described  in  Appendix  I  will  be  used.  The  aerodynamic  coefficients 
and  rotary  derivatives  will  be  calculated  for  this  body  through  an  angle  of  attack  and 
angle  of  sideslip  range. 

a.  Description  of  Body  Coordinate  System 

The  axis  system  used  to  describe  the  body  is  shown  in  Figure  58.  The  Coordinates 
are  body  axes  with  the  x-axis  aligned  along  the  body.  The  exact  location  of  the  x-axis 
is  chosen  to  permit  the  body  cross  sections  to  be  obtained  in  planes  perpendicular  to 
this  axis.  The  exact  location  of  the  origin  is  not  restricted  to  be  at  the  body  nose  but 
may  be  chosen  to  suit  the  user.  The  axis  system  chosen  consists  of  a  right  hand  system 
with  the  x-axis  directed  aft,  the  y-axis  to  the  right  and  the  z-axis  upward. 

I 

The  flight  conditions  for  the  static  coefficients  are  specified  a3  a  resultant  angle 
of  attack  a  and  a  roll  angle  <p.  The  resultant  angle  of  attack  is  the  angle  between  the 
freestream  direction  and  the  x-axis  and  is  always  defined  as  positive.  The  roll  angle 
is  then  specified  as  the  angle  between  the  freestream  componeht  in  the  y-z  plane  and 
the  body  vertical  plane  (x-z  plane)  as  shown  in  Figure  58. 
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FIGURE  58.  SIGN  CONVENTION  FOR  BODY  AERODYNAMIC  COEFFICIENTS 


Of  tfae  rotary  velocities  p  is  not  considered  of  importance  for  a  body  and  so  is 
not  included  in  the  computation.  The  rotary  components  q  and  r  are  specified  in  body 
zyjz*  with  q  as  a  vector  in  the  direction  of  the  positive  y-axis  and  r  as  a  vector  in  the 
direction  of  the  negative  z  axis.  This  convention  was  chosen  to  be  consistent  with  the 
terminology  of  Reference  6. 

A  reference  length  1  reference  area  S  and  a  center  of  gravity  location  x _ 

r  r  eg 

are  specified.  All  of  the  for^e  coefficients  are  based  on  the  same  reference  area. 
Pitching,  yawing  and  rolling  moment  coefficients  are  defined  on  the  reference  area 
and  the  reference  length  specified.  This  differs  from  the  conventional  way  of  defining 
the  moments.  The  reference  point  about  which  the  moments  are  taken  is  the  specified 
center  of  gravity  location  which  is  located  on  the  x-axis.  The  effects  of  including 
rotary  velocities  assume  that  the  center  of  rotations  of  q  and  r  are  at  xc^. 

The  force  and  moment  coefficients  are  in  body  axes  as  shown  in  Figure  58, 
positive  along  the  positive  z-axis,  Cm  positive  for  a  moment  pitching  the  nose  up:. 

Cy  is  positive  in  the  positive  y  direction  and  CQ  positive  for  a  nose  right  moment. 

Rolling  moment  coefficient  Cj  is  defined  positive  counterclockwise  the  moment  being 
specified  about  the  x-axis.  No  attempt  has  been  made  to  incorporate  the  axial  force 
coefficient  as  the  method  used  is  not  suitable  for  that  purpose. 

b.  Body  Description  for  Nonlinear  Force  and  Moment  Program 

To  use  the  computer  program  for  aerodynamic  coefficients  the  body  must  be 
described  for  a  series  of  sections  taken  perpendicular  to  the  x-axis.  It  is  not  necessary 
to  take  the  sections  chosen  at  equal  intervals  but  the  spacing  should  be  relatively  uni¬ 
form  with  more  sections  being  taken  in  regions  where  the  cross  sectional  parameters 
are  changing  rapidly. 

The  section  inputs  include  both  the  geometrical  variables  S  and  dS/dx  and  the 
coefficients  of  the  mapping  function  and  their  derivatives.  It  is  therefore  necessary 
to  know  the  mapping  of  each  of  the  sections  being  inputed.  Although  it  is  necessary  to 
include  mapping  coefficients,  the  nature  of  the  slender  body  theory  is  such  that  only 
the  first  few  coefficients  are  of  primary  significance.  For  this  reason  it  is  possible 
to  approximate  the  coefficients  of  the  mapping  function  and  still  retain  reasonable 
accuracy.  Therefore,  methods  are  presented  for  obtaining  the  mapping  coefficients 
accurately  and  an  easier  approximate  method  which  retains  the  more  significant 
coefficients. 
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The  coefficients  as  obtained  by  the  mapping  program  must  be  modified  before 
they  are  suitable  for  use  with  the  nonlinear  body  aerodynamics  program  and  the 
method  of  modification  will  be  described.  The  simplified  method  of  obtaining  coefficients 
for  the  mapping  function  will  also  be  described  and  a  complete  set  of  inputs  to  die  pro¬ 
gram  will  be  given  for  this  simplified  method. 

c.  Modification  of  Mapping  Function  for  Body  Aerodynamics  Program 

Equation  (1)  of  Section  n  is  not  in  the  proper  form  for  use  with  the  body  aero¬ 
dynamics  program.  When  this  equation  is  rewritten  to  include  the  constant  term,  i.e., 
to  locate  the  section 


On  +  ibn 


(3) 


the  secH  -1  -mated  and  mapped  as  shown  in  Figure  59(a).  For  the  body  aerodynamics 

progran  c.uc  tion  is  rotated  as  shown  in  Figure  59(b)  and  the  mapping  is  commenced 
at  a  diii.ri.ut  point  jn  the  section.  In  addition,  instead  of  basing  the  mapping  on  a 
circle  of  radius  rc,  t  ^  mapping  is  rewritten  to  base  the  new  mapping  on  the  unit  circle 
coordinate  a=  e1^.  final  form  of  the  mapping  is  then: 


2  ~tc(T  +  c0 


+  i(jn  ,  . . *  + 

fia°  +  +  -g-  + 


(4) 


where  these  coefficients  are  related  to  the  coefficients  of  Equation  (3)  by  the  relation 


Cn  *  * dn 


(6n+  ibn) 


(5) 


For  a  symmetrical  shape  this  reduces  to  the  relation 


T7  <3«  i  a  p 

c,  id  p 
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a.  Original  Mapping  Relation 


Z  -  Plane 


cl+  idl 

Z  =  r  o  +  C-  + id-  +  - — —  + 

c  0  0  cr 


CN+  ldN 
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b.  Final  Mapping  Relation 


FIGURE  59.  CHANGE  OF  MAPPING  FUNCTION  FOR  BODY 
AERODYNAMICS  PROGRAM 
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d.  Simplified  Handbook  Method  for  Obtaining  Coefficients 

Since  it  is  desired  to  use  the  body  aerodynamics  program  for  preliminary  design 
type  work  it  is,  where  possible,  desirable  to  avoid  the  complexity  of  obtaining  the 
mapping  function.  It  is  possible  to  do  this  for  the  usual  fuselage  shapes  encountered 
and  still  retain  sufficient  accuracy  for  preliminary  design  purposes. 

The  actual  fuselage  is  replaced  by  an  equivalent  body  in  which  the  sections  an 
replaced  by  equivalent  ellipses  keeping  the  same  body  camber.  This  replaces  the  me.*, 
ping  of  Equation  (6)  above  with  the  truncated  mapping 

Zmtca  +  id0+£l-  (7) 


This  expression  retains  the  most  critical  terms  in  the  mapping  as  far  as  obtaining  the 
body  aerodynamics  and  it  is  possible  to  approximate  the  three  coefficients  rc,  dQ  and 

cr 

Defining  a  to  be  half  of  the  maximum  vertical  dimension  of  the  true  cross  sec¬ 
tion  and  b  to  be  half  the  maximum  lateral  dimension  of  the  section,  it  is  possible  to 
approximate  r  and  c1  by  the  expressions 

C  1 


the  results  being  reasonably  accurate  for  fuselages  which  do  not  depart  too  far  from 
the  elliptical.  The  coefficient  dQ  can  be  replaced  by  the  centroid  of  the  sectional  area. 

The  remainder  of  the  inputs  required  for  the  slender  body  portion  of  the  computer 
program  are  quantities  which  can  be  obtained  directly  from  the  body  geometry,  such 
as  the  cross  sectional  areas  and  its  derivative  with  respect  to  x. 

e.  Viscous  Cross  Flow  Input 

The  computer  program  also  requires  that  a  cross  flow  drag  coefficient  be  input 
for  both  the  components  of  flow  in  the  vertical  plane  and  the  lateral  plane. 
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The  viscous  crossflow  terms  can  be  obtained  by  using  the  drag  coefficients  of 
an  infinite  two-dimensional  elliptical  cylinder  at  each  body  station,  the  character  of 
the  ellipse  changing  according  to  the  maximum  dimensions  of  the  body  as  described 
above. 

These  drag  coefficients  can  be  obtained  using  the  drag  coefficient  given  in 
Reference  9  for  two  dimensional  subcritical  ellipses  and  given  by 

Cv-.0?5(l+f)i-t.1%  (10) 

where  t  =  maximum  dimension  perpendicular  to  crossflow 

c  =  maximum  dimension  parallel  to  crossflow 

A  coefficient  is  computed  for  a  crossflow  velocity  component  in  the  vertical 
c'irection  in  one  case  and  a  second  coefficient  is  computed  for  a  crossflow  component 
in  the  lateral  plane.  These  two  coefficients  are  then  multiplied  by  the  maximum 
dimensions  perpendicular  to  the  flow  direction  for  the  input  parameters. 
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2.  SAMPLE  COMPUTATIONS  FOR  TEST  BODY 


Table  IV  shows  a  set  of  computations  for  rc  and  for  die  wind  tunnel  test 
model  of  this  study.  The  computations  are  straightforward  involving  no  difficulties. 
The  major  inconvenience  in  developing  the  inputs  for  the  program  is  in  obtaining  S 
and  the  derivatives  of  each  of  the  coefficients  with  respect  to  body  station.  There  is 
no  conceptual  difficulty  involved  but  the  required  integration  for  the  body  cross  sec¬ 
tional  area  and  the  centroid  location  and  the  graphical  differentiations  needed  to  obtain 
dS/dx,  drc/dx,  ddQ/dx  and  dCj/dx  is  tedious.  If  necessary,  computer  programs  can 
be  written  to  do  the  necessary  integration  and  differentiations  but  such  programs  are 
not  included  here. 

a.  Sample  Inputs  for  Force  and  Moment  Program 

Figure  CO  shows  a  sample  set  of  inputs  for  the  body  aerodynamics  program.  The 
data  are  for  the  wind  tunnel  test  model  tody  of  this  study  contract  with  the  canopy  off. 

Card  1  specifies  the  maximum  number  of  mapping  coefficients  for  any  section 
input  (maximum  of  12)  and  the  number  of  stations  for  which  section  data  is  input 
(maximum  of  40) . 

Cards  2—4  give  the  station  locations  of  the  input  sections,  maximum  of  40.  The 
remainder  of  the  cards  must  be  in  units  consistent  with  these  numbers.  The  stations 
input  must  not  include  the  nose  location  nor  the  tail  section  if  these  stations  have  zero 
area  and  mapping  circle  radius. 

Cards  5—7  contain  the  radii  of  the  mapping  circle  at  the  input  stations. 

dr 

Cards  8—10  give  the  values  of  for  the  same  stations. 

Cards  11—13  are  the  cross  sectional  areas  S  of  the  sections. 

Cards  14—16  are  the  values  of  dS/ ,  . 

dx 

Cards  17—19  give  the  values  of  the  side  viscous  crossflow  drag  coefficient  per 
unit  length  times  the  maximum  vertical  dimension  at  the  section. 

Cards  20—22  give  the  vertical  viscous  crossflow  drag  coefficient  times  the 
maximum  lateral  section  dimension. 

Cards  23—76  consist  of  sets  of  three  cards  for  each  of  the  input  sections  (in  this 
case  18  sets).  The  first  card  of  this  set  contains  two  numbers.  The  first  of  these 
specifies  the  number  of  mapping  coefficients  of  the  given  section.  If  the  number 
specified  is  zero,  the  program  uses  the  number  given  on  Card  1.  The  number  given 
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TABLE  IV.  COMPUTATIONS  FOR  WIND  TUNNEL 
TEST  MODEL 


X 

y/2 

z- 

z  * 

r 

c 

C1 

7.25 

4.25 

-14.1 

-5.6 

4.25 

0 

23.7 

9.2 

-18.7 

3.1 

10.05 

-.85 

41.0 

13.1 

-22.2 

11.0 

14.85 

-1.75 

73.0 

17.4 

-26. 25 

24.0 

21.2625 

-3. 8625 

94.0 

19.25 

-27, 95 

31. 3b 

21.45 

-5.2 

118.0 

21.0 

-29.2 

37.75 

27.2375 

-6.2375 

143.5 

22.45 

-30.2 

40.7 

28.95 

-6.5 

163.5 

23.45 

-30.55 

40.7 

29.5375 

-6.0875 

185.5 

24.3 

-30.8 

40.7 

30.025 

-5. 725 

221.5 

25.05 

-30.6 

40.7 

30.35 

-5.3 

264.25 

25.1 

-29.3 

40.7 

30.05 

-4.95 

316.0 

24.8 

-25.0 

40.7 

28.825 

-4. 025 

343.0 

24.6 

-21.5 

40.4 

27.775 

-3.175 

374.0 

24.45 

-16.3 

39.75 

26.2375 

-1.7875 

411.0 

21.5 

-9.1 

38.3 

22.6 

-1.10 

450.0 

IS.  6 

-.85 

35.95 

17.0 

-1.4 

497.0 

6.4 

11.6 

30.8 

8.0 

-1.6 

512.0 

3.545 

16. 11 

27.5 

4.62 

-1.075 
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FIGURE  60.  (Continued) 


304JZ6£.  0  /1- •  53000E  J> 


157 


must  sot  exceed  the  number  on  the  first  card.  The  second  number  specifies  whether 
or  not  the  section  under  consideration  is  symmetrical  or  not.  If  the  section  is  sym¬ 
metrical,  the  number  0  is  input;  otherwise  1.  The  axis  of  symmetry  is  assumed  to  be 
the  vertical  plane. 

The  second  card  for  each  set  specifies  the  coefficients  of  the  mapping  function. 
For  a  symmetrical  section  the  coefficients  input  follow  the  order  d^,  Cj.  d2,  c3>  d4, 
c5,  ...  since  the  other  coefficients  are  zero.  For  an  unsymmetrical  section  the 
coefficients  are  specified  In  the  order  cQ,  dQ,  c^  dj,  c2,  c^,  ...  up  to  the  maxi¬ 
mum  number  of  pairs  specified.  The  subscript  0  specifies  the  second  coefficient  (or 
coefficient  pair),  the  subscript  1  the  third  coefficient,  etc. 

The  remaining  cards  specify  the  flight  conditions  under  which  the  coefficients 
are  to  be  found.  Card  77  is  a  comment  card  and  can  contain  any  pertinent  informa¬ 
tion  desired. 

Card  78  specifies  four  numbers:  reference  length,  1^,  reference  area  Sr,  center 
of  gravity  (and  moment  center)  location  x  and  the  incremental  step  size  along  the 
body  at  which  computations  are  made.  The  program  assumes  linearity  between  the 
incremental  steps  here  specified  so  that  a  reasonably  large  number  of  steps  are  re¬ 
quired  along  the  body. 

The  79th  card  specifies,  respectively,  the  number  of  angles  of  attack  to  be 
computed  (maximum  of  18),  the  number  of  roll  angles  (maximum  of  9),  number  of 
pitching  velocities  (maximum  of  9)  and  number  of  yawing  velocities  (maximum  of  9). 

Cards  80  and  81  specify  die  angle  of  attack  at  which  the  coefficients  are  to  be 
evaluated.  These  angles  are  in  degrees. 

Card  82  specifies  the  roll  angles  which  are  to  be  computed,  also  in  degrees. 

Card  83  specifies  the  desired  values  of  pitching  velocity  inputs.  The  number 
inputed  represents  (dimensionless). 

Similarly  Card  84  specifies  yawing  velocity  inputs  specified  as  rl  /OT. 

r  c\j  oo 

(dimensionless) . 

This  completes  the  input  cards  needed  to  compute  the  nonlinear  body  aero¬ 
dynamics.  Cards  are  added  or  subtracted  as  necessary  to  input  all  die  specified  data. 
That  is,  enough  cards  are  used  to  input  the  numbers  required  and  no  blank  cards  are 
to  be  inputed.  As  an  example,  the  number  of  stations  here  specified  is  .18  which 
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requires  3  cards  to  specify  each  of  the  first  set  of  parameters.  More  of  less  cards 
would  be  used  depending  on  the  number  of  stations  specified. 


b.  Sample  Outputs  from  Force  and  Moment  Program 

Figure  61  shows  the  output  from  the  computer  program  for  the  inputs  listed  in 
Figure  60.  The  first  line  written  out  in  the  information  input,  on  the  comment  card. 

/qlr  \  /rlr  \ 

The  second  line  lists  the  flight  conditions  PHI  (<p),  Q\ 20^/’  *or 

which  the  coefficients  are  calculated. 


After  this  are  tabulated  the  five  coefficients  C  .  C  ,  C  ,  C  and  CL  in  that 

n-  m  v  n  .  1 

order  as  functions  of  angle  of  attack.  The  coefficients  as  written  out  are  separated 
into  a  potential  component  (obtained  by  slender  body  theory)  and  a  viscous  component 
(using  viscous  crossflow).  To  obtain  the  coefficient  these  two  components  must  be 
added  together.  The  program  does  not  calculate  a  viscous  component  to  the  rolling 
moment  so  this  is  printed  out  as  zero. 

When  more  than  one  set  of  flight  conditions  (other  than  a)  are  input  this  tabula¬ 
tion  is  repeated. 
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FIGURE  81.  NONLINEAR  BODY  AERODYNAMICS  PROGRAM  OUTPUT 
DATA  FOR  SAMPLE  PROBLEM 


3.  COMPARISON  WITH  TEST  DATA 


Figure  62  shows  the  comparison  between  theory  and  test  for  normal  force  co¬ 
efficient  C^and  pitching  moment  coefficient,  Cm  at  zero  sideslip.  The  theoretically 
predicted  values  of  are  somewhat  low.  In  particular,  the  predicted  lift  curve 
slope  at  zero  angle  of  attack  is  predicted  to  be  zero  while  the  test  results  show  a  finite 
value.  The  value  of  at  zero  angle  of  attack  is  also  in  error.  The  agreement  obtained 
for  the  pitching  moment  can  be  considered  to  be  good. 

Figure  63  shows  a  comparison  between  test  and  theoretical  side  force  C  and 
yawing  moment  C^.  The  theoretical  side  force  tends  to  be  somewhat  low,  and  the 
agreement  in  yawing  moment  is  somewhat  worse  than  the  pitching  moment  agreement. 

The  computer  program  does  not  calculate  the  derivatives  but  it  is  possible  to 
compute  two  sets  of  coefficients  at  different  values  of  a  parameter  and  obtain  the 
derivative  by  dividing  the  difference  of  the  coefficients  by  the  incremental  change  in 
the  parameter. 

Figures  64  and  65  show  samples  of  rotary  derivatives  obtained  using  the  computer 
program.  No  test  data  is  available  with  which  to  compare  the  theoretical  results, 
so  it  is  not  possible  to  predict  what  accuracy  is  obtained  by  the  theoretical  method. 


Test  Points 


FIGURE  62.  THEORY  AND  TEST  CONDITIONS  OF  CXT  AND  C  FOR  WIND  TUNNEL 

N  m 

TEST  MODEL  FUSELAGE 
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TEST  MODEL  FUSELAGE 


SECTION  vn 


NONLINEAR  WING  AERODYNAMICS 

The  method  used  for  predicting  the  nonlinear  aerodynamics  of  wings  is  that 
described  in  Section  VII  of  Volume  I.  This  method  has  been  programmed  for  the 
computer  and  essentially  converts  the  known  section  characteristics  for  the  wing 
into  finite  wing  characteristics.  First  the  wing  section  characteristics  have  to  be 
identified  from  available  test  data.  This  information  is  then  used  to  obtain  the 
positions  and  relative  strengths  of  the  two  lifting  lines.  The  wing  planform  and 
flight  condition  are  included  with  the  section  model  into  the  computer  program  to 
enable  the  wing  aerodynamics  to  be  determined. 

1.  SAMPLE  PROBLEM 

/ 

To  illustrate  the  use  of  the  nonlinear  wing  aerodynamics  program  the  wing  for 
the  wind  tunnel  test  model,  described  in  Appendix  I,  will  be  used. 

The  wing  planform  is  shown  in  Figure  66  with  the  axes  system,  lifting  lines 
(discussed  later)  and  downwash  control  line  indicated.  The  wing  employed  a 
NACA  63 AO  10  section. 


FIGURE  66.  WING  PLANFORM  WITH  LIFTING  LINES  AND  DOWNWASH  CONTROL 

LINE  SHOWN 


The  axis  system  chosen  is  a  right  hand  system  with  the  x-axis  directed  aft  and  the 
z-axis  directed  upwards.  The  flight  conditions  for  the  sample  problem  will  be  a=  14° 

and  16°,  13=  0°  with  the  rotrry  variables  p,  q,  r  all  equal  to  zero. 

a.  Determination  of  Section  Parameters 

The  separation  characteristics  of  airfoil  sections  are  a  function  of  the 
Reynolds  Number,  R,  the  flow  tending  to  separate  at  a  lov/er  value  of  a , 
for  lower  values  of  R.  The  unpowered  runs  in  the  tests  described  in  Appendix 

c 

I  were  at  a  Reynolds  Number/foot  =  1. 3  x  10  .  No  data  could  be  found  for  the 
test  model  wing  section  (NACA  63A010)  for  Reynolds  Numbers  as  low  as  the 
test  value.  However,  data  for  section  NACA  64 ^ -012,  with  a  similar  section 
geometry,  for  the  required  values  of  R  has  been  determined  in  Reference  10  and 
is  tabulated  in  Table  V. 

The  drag  coefficient  for  the  airfoil  at  a  =  90°  is  taken  to  be  2. 08  (value  for 
NACA  0012  airfoil)  with  the  line  of  action  passing  through  the  chord  center.  The 
leading  lifting  line,  in  the  mathematical  model  for  the  airfoil  section,  is 
positioned  at  the  quarter  chord  position  ( Cmi =  0  in  linear  a  range).  The 
downwash  control  point  is  taken  to  be  at  .  75  chord.  The  aft  lifting  line  is  now 
chosen  to  give  a  good  fit  to  the  section  pitching  moment  while  satisfying  the 
boundary  condition  of  no  flow  through  the  section  at  the  downwash  control 
station  and  providing  an  exact  duplication  of  the  section  normal  force.  In  this 
case,  by  taking  the  aft  lifting  line  to  be  located  at  .  70  of  the  chord,  a  satis¬ 
factory  representation  to  the  section  characteristics  is  possible  as  shown  in 
Table  V.  The  weighting  function  Wt  (weighting  of  circulation  between  leading 
and  aft  lifting  lines)  is  also  shown  in  Table  V.  The  value  of  W^.  is  needed  for 
ot  =  o  and  this  is  determined  by  extrapolation  of  the  values  for  larger  a. 

b.  Inputs  to  Nonlinear  Wing  Aerodynamics  Program  for  Sample  Problem 
A  sample  set  of  inputs  for  the  nonlinear  wing  aerodynamics  program  is 

shown  in  Figure  67. 

Card  1  specifies  the  initial  value  of  wing  angle  of  attack  a,  angle  of  sideslip 
(3 ,  and  the  step  size  in  a  all  in  degrees. 

Card  2  provides  wing  planform  information  in  the  form  of  the  y-coordinates 
(relative  to  the  wing  root  chord),  of  the  wing  root  and  tip  chords,  wing  taper 
ratio,  and  the  tangent  of  the  leading  edge  sweep  angle. 

Card  3  gives  in  order,  the  rolling,  pitching  and  yawing  velocities. 

Card  4  contains  the  reference  length  (relative  to  the  wing  root  chord),  the  x- 
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TABLE  V.  SECTION  DATA  FOR  AIRFOIL  SECTION 


NACA  64  a  -012 


TEST 

LIFT 

TNG  LINE  METHOD 

a 

CN 

Cm  1/4 

CN 

Cm  1/4 

wt 

0 

0 

0.0 

0. 

0. 

Indet. 

2 

0.18 

0.0 

0.18 

0.974 

4 

0.36 

0.0 

0.36 

-0.007 

0.972 

6 

0.54 

0.0 

0.54 

-0.013 

0.970 

7 

0.63 

0.63 

-0.016 

0.969 

8 

0.71 

-0.008 

0.71 

1 

0.967 

9 

0.78 

-0.010 

0.78 

0.962 

10 

0.83 

-0.015 

0.83 

-0.032 

0.953 

12 

0.83 

1 

0.83 

-0.051 

0.914 

14 

0.68 

0.68 

-0.079 

0.807 
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FIGURE  67.  NONLINEAR  WING  AERODYNAMICS  PROGRAM  INPUT 
DATA  FOR  SAMPLE  PROBLEM 


and  z-  coordinates  of  the  center  of  gravity  location. 

Card  5  contains  the  wing  section  normal  force  coefficient  for  a  =  90°  and  the 
x -coordinate  of  the  intersection  of  its  line  of  action  with  the  airfoil  chord. 

Card  6  gives  the  number  of  circulation  and  downwash  control  points. 

Card  7  specifies  the  number  of  angle  of  attacks  to  be  computed  and  the 
number  of  iterations  for  each  a  . 

Card  8  indicates  whether  the  wing  loading  is  symmetrical  about  the  x-axis, 
in  this  case  NSYM  =  0  indicating  symmetry.  (NSYM  =  1  indicates  asymmetrical 
loading.) 

Cards  9  and  10  list  the  y-coordinates  of  the  circulation  control  points  and 
downv/ash  control  points . 

Cards  11,  12  and  13  specify  the  x-coordinates  (relative  to  chord)  and  the 
tangents  of  the  sweep  angles  for  the  leading  lifting  line,  the  aft  lifting  line  and 
the  line  connecting  the  downwash  control  points. 

Card  14  contains  the  effective  angles  of  attack  for  the  downwash  control 
point  stations.  In  this  case  the  values  have  been  determined  from  the  previous 
calculations  for  a  =  12°. 

Cards  15  and  16  list  the  values  of  oi  in  degrees  at  which  the  circulation 
weighting  function  is  tabulated. 

Cards  17  and  18  contain  the  tabulated  values  of  the  circulation  weighting 
function  determined  in  this  case  from  Table  V. 

c.  Outputs  to  Nonlinear  Wing  Aerodynamics  Program  for  Sample  Problem 

Figure  ;8  shows  the  output  from  the  computer  program  for  the  input  listed 
in  Figure  67.  The  flight  condition  is  shown  in  the  form  of  values  for  m  ,  •>  and 
p,  q,  r.  The  spanwise  loading  and  effective  angle  of  attack  at  the  spanwise 
stations  are  shown  next.  The  normal  force  coefficient  and  the  moments  about 
the  y-  and  x-axes  follow.  This  output,  excepting  the  flight  condition  variables, 
is  repeated  for  the  number  of  iterations  on  a.  (in  this  case  2), 

This  set  of  output  is  then  repeated  for  the  number  of  a’s  input  to  the  program 
which  in  this  case  is  two. 

d.  Method  Applicability  and  Limitations 

The  method,  as  presently  programmed,  Is  restricted  to  straight  tapered 
wings.  Because  the  approach  uses  lifting  lines  the  method  Is  really  applicable 
to  large  aspect  ratio  wings  and  the  accuracy  of  the  predictions  will  not  be  as 
good  for  lower  aspect  ratio  wings.  Flaps  may  be  accounted  for  by  changing  the 
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sectio  characteristics  (cr  the  wing.  Calculations  of  the  finite  wing  aerodyna¬ 
mic  characteristics  depend  on  a  knowledge  of  the  section  characteristics  so  that 
in  some  cases  this  may  be  a  limitation  on  the  method.  The  method  is  program¬ 
med  to  calculate  the  effects  of  sideslip  and  the  rotary  derivatives  although  no 
significant  attempt  has  been  made  to  validate  these  options. 

In  typical  calculations  two  iterations  are  used  in  the  linear  range  of  a  and 
five  iterations  for  nonlinear  or.  Estimates  of  the  effective  angle  of  attack  are 
determined  from  previous  computations  at  lower  values  of  a.  The  initial  calcu¬ 
lations  have  usually  started  with  a  =  4°  and  effective  alpha  equal  to  2°. 

2.  COMPARISON  WITH  TEST  DATA 

The  calculations  shown  in  Figure  68  are  for  the  test  model  described  in 
Appendix  I.  The  test  data  for  these  conditions  for  comparison  were 

l 

a  =  14°  :  CN  =  .  708  ,  Cj^  =  .  365 
<*  =  16°  :  CN  =  .665  ,  ^^  =  .375 

Further  calculations  for  an  aspect  ratio  6  wing  are  shown  in  Figure  69  compared  with 
the  test  data  of  NASA  TM  2-27-29A. 
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Revolts  for  Alpha  *  14.0000,  and  Beta  «  O.i  JO?  Ccgrecs 
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FIGURE  68.  SAMPLE  OUTPUTS  FOR  NONLINEAR  WING 
AERODYNAMICS  PROCRAM 
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- —Calculation 

®  Test  Data 

(NASA  TM  2 -27 -5 9 A) 


FIGURE  69.  WING  CALCULATIONS  FOR  AN 
ASPECT  RATIO  6  WING 
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APPENDIX  I 


WIND  TUNNEL  TESTING  OF  V/STOL 
CONFIGURATION  MODEL 

A  model  was  constructed  and  tested  to  supplement  the  analytical  investigat  ; 
and  gain  further  data  for  use  in  validation  and  improvement  of  the  analytical  predic¬ 
tion  techniques.  The  model  was  configured  to  resemble  a  feasible  military  airci  ift 
but  was  designed  to  operate  with  jets  of  varied  number  and  location  to  allow  a  v  riety 
of  data  to  be  generated.  The  test  contributes  new  data  because  of  extensive  pressure 
instrumentation  present  on  the  model.  These  data  facilitate  identification  of  the 
sources  of  the  various  induced  loads  measured  during  the  test. 

1.  MODEL  AND  APPARATUS 

The  model  is  a  shoulder  wing  configuration  equipped  with  an  external  airfoil  flap 
and  a  stabilator  mounted  on  the  vertical  tail  above  the  fuselage.  Two  vectored  thrust 
engines  are  contained  in  large  nacelles  mounted  beneath  the  wing  adjacent  to  the  fuse¬ 
lage;  a  single  lift  engine  is  mounted  within  the  body  and  forward  of  the  vectored  exits. 
The  general  arrangement  of  the  model  is  shown  in  the  three  view  drawing  of  Figure 
1-1  and  the  photographs  of  Figures  1-2,  1-3,  and  1-4.  Detailed  geometrical  data  are 
given  in  Table  I-I. 

The  vectored  thrust  engines  were  simulated  by  ejector  type  jet  engine  simula¬ 
tors.  Two  3.67-inch  diameter  nozzles  v/ere  tested  at  two  longitudinal  positions  (11 
and  111  percent  c)  at  nominal  deflections  of  0,  45,  or  90  degrees  at  the  forward  loca¬ 
tion,  and  45  or  90  degrees  at  ihe  aft  location.  Larger  4.5-inch  diameter  nozzles  were 
tested  in  the  aft  position  at  a  nominal  deflection  of  90  degrees.  Plugs  were  used  to 
seal  the  inlets  when  -nlet  flow  was  not  desired. 

Typical  operating  curves  for  the  ejectors  are  shown  in  Figure  1-5.  The  design 
of  the  ejector  nozzles  produced  a  relatively  nonuniform  jet  velocity  profile,  shown  in 
Figure  1-6. 


Preceding  page  blank 
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The  lift  jet  was  sinmlated  by  a  convergent  nozzle  supplied  through  a  perforated 
plate.  The  exit  diameter  was  2.25  inches.  No  inlet  simulation  was  provided.  The  lilt 
jet  possessed  a  relatively  uniform  exit  profile,  shown  in  Figure  I-G. 

Details  of  jet  exit  location  and  point  of  application  of  resultant  forces  are  shown 
in  Table  I-H  in  terms  of  location  along  the  mean  aerodynamic  chord  and  distance  below 
the  wing  plane. 

The  engine  simulators  were  driven  by  cold  dry  air  supplied  through  a  common 
plenum  fed  by  a  flexible  metal  tube  passing  through  the  sting. 

Model  forces  and  moments  were  recorded  using  a  six-component  internal  balance. 

Two  hundred  and  sixty-four  pressure  taps  are  present  on  the  left  half  of  the 
model.  They  are  placed  in  four  groups:  a  wing  pattern,  a  lower  fuselage  pattern, 
a  circumferential  pattern  at  five  fuselage  stations,  a  nacelle  centerline  pattern. 

Lower  fuselage  and  wing  pressure  patterns  are  described  in  Tables  I-DI  and  I-IV. 
Circumferential  fuselage  and  nacelle  centerline  patterns  are  shown  in  Figures  1-7 
and  1-8. 

A  seven-probe  flow  angularity  rake  was  sometimes  mounted  at  the  tail  station 
in  place  of  the  empennage.  Its  general  arrangement  may  be  seen  in  Figure  1-2. 

Testing  was  performed  in  the  NASA  Langley  V/STOL  tunnel  which  has  a  test 
section  of  14  X  22  feet. 

2.  TEST  PROCEDURE 

Testing  of  the  model  was  performed  in  two  phases:  calibration  of  balance  and 
engine  simulators,  and  aerodynamic  testing  within  the  wind  tunnel. 

a.  Balance  Calibration  and  Corrections 

The  testing  of  a  powered  model  requires  additional  balance  calibration  and  cor¬ 
rection  beyond  that  routinely  performed  during  unpowered  testing.  The  air  supply 
balance  arrangement  used  during  this  test  is  indicated  schematically  in  Figure  1-9. 

Although  the  air  supply  line  is  designed  to  be  highly  flexible  and  cause  minimum 
interference  with  force  measurements,  corrections  must  be  applied  to  balance  mea¬ 
surements  to  reflect  that  portion  of  the  total  load  which  is  carried  by  the  line.  A 
series  of  known  loads  was  applied  to  the  balance-air  supply  line  assembly,  and  the 
resultant  balance  measurements  were  used  to  obtain  a  matrix  of  linear  correction 
coefficients. 
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Pressurization  of  the  system  caused  forces  to  be  exerted  on  the  model  by  the 
flexible  supply  line.  These  forces  were  measured  and  calibrated  with  the  plenum 

sealed  arid  no  nongravitational  loads  applied  to  the  model.  The  calibrated  forces  have 
been  removed  from  final  force  data. 

The  final  form  of  the  corrections  is  shown  below.  The  momentum  input  due  to 
the  supply  mass  flow  was  found  to  be  negligible. 

'■’w  rw 

F  =  K  F.  .  -  F  (P  ) 
cor  bal  p'  sup 


where 

a  6x1  matrix  of  applied  loads 

a  6x6  matrix  of  correction  coefficients 

a  6x1  matrix  of  loads  measured  by  balance 

a  6x1  matrix  of  loads  applied  by  the  air  supply  line,  calibrated 
against  supply  pressure 

Additional  details  concerning  the  magnitude  and  accuracy  of  the  corrections  may 
be  found  in  Appendix  II. 

b.  Calibration  of  Engine  Simulators 

Each  engine  simulator-nozzle  combination  was  tested  individually  to  determine 
direct  thrust  force  and  moment  applied  to  the  model.  The  same  air  supply-balance 
assembly  used  during  tunnel  tests  was  used  for  this  calibration.  A  calibrated  bell- 
mouth  was  used  to  monitor  ejector  inlet  flow.  A  limited  survey  was  made  of  nozzle 
exit  profiles  to  determine  their  basic  character. 

These  tests  were  performed  with  the  units  mounted  on  the  bare  model  plenum  to 
minimize  static  interference  effects  by  reducing  surface  area  near  the  exits.  Some  sma’l 
interference  is  of  necessity  present  in  the  data  because  of  the  physical  proximity  of 
the  external  portions  of  the  nozzles  and  drive  system. 

Lift  jet  thrust  was  calibrated  on  the  basis  of  total  pressure  within  the  nozzle. 
Ejectors  were  calibrated  as  a  function  of  primary  nozzle  plenum  pressure.  No  cor¬ 
rections  were  made  to  ejector  calibrations  due  to  forward  speed  present  in  the  tunnei. 
Erior  due  to  this  approximation  is  examined  in  Appendix  III. 
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c.  Wall  Corrections 

No  corrections  were  made  to  the  data' for  possible  wall  effects.  The  dimensions 
of  the  14  X  27,  foot  test  section  are  quite  large  in  comparison  to  basic  model  dimen¬ 
sions,  and  computations,  made  prior  to  the  test  using  Heyson's  method  indicated 
that  wail  interference  effects  would  be  within  the  accuracy  of  the  data  acquistion 
system.  However,  the  pitch  mechanism  used  during  the  test  was  not  capable  of 
maintaining  the  model  at  a  fixed  position  within  the  test  section,  causing  the  model  to  be 
nearer  the  floor  at  the  lowest  angle  of  attack  tested.  An  investigation  was  made  to 
determine  if  the  ground  effect  of  the  floor,  the  nearest  surface,  was  causing  any 
significant  effect  on  the  model.  The  wing-body-nacelle-tail  was  equipped  with  the 
large  exits,  pitched  to  a  6  degree  angle,  and  tested  at  various  forward  speeds  at 
different  heights.  An  effect  was  observed  a  the  lowest  forward  speed  tested  but  not 
at  higher  speeds  nor  statically  (see  Figure  1-10).  The  minimum  rn.del  height  during 
norma1,  testing  was  42  inches.  The  nondimensionalized  minimun.  neight  shown  below 
indicates  that  the  large  exits  will  show  the  most  severe  interference  at  lower  angles 
of  attack  and  low  forward  speed.  The  order  of  magnitude  may  h>.'  t  to  three  percent 
of  total  load. 
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d.  Test  Parameters 


The  testing  of  powered  models  creates  a  requirement  for  a  parameter  relat’ig 
propulsive  and  aerodynamic  forces.  The  effective  \elocily  ratio  as  used  in  this  series 
of  tests  is  the  square  root  of  the  ratio  of  freestream  unit  momentum  flux  to  mean  jet 
unit  momentum  flux.  It  is  obtained  from  the  following  expression: 
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Another  parameter,  a  relative  measure  of  propulsive  and  aerodynamic  forces, 
is  the  thrust  coefficient  -  the  ratio  of  thrust  to  the  product  of  freestream  dynamic 
pressure  and  wing  area. 
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The  two  paramotorB  aro  rolatod  In  tho  manner  shown  below, 
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Tho  offootivo  volocity  ratio  Is  used  no  tho  primo  parameter  of  thiB  roport 
bocauso  it  Is  loss  rolatod  to  tho  goomotry  of  tho  particular  configuration  being  con¬ 
sidered,  and  it  bettor  doscribcs  the  state  of  tho  highly  dofloctod  Jot. 


o.  Tost  Program 

Tho  model  was  constructed  ouch  that  nacolloa,  lift  ongino,  wing,  and  omponnago 
could  bo  mounted  on  tho  fusolago  in  any  combination.  Comploto  modol  buildups  woro 
performed  for  tho  wing-body-nacollo-tuil  configuration  poworod  by  tho  two  ojoctors 
in  tho  forward  position,  and  for  tho  wing-body-tail  poworod  by  tho  lift  Jot,  Emphasis 
was  placed  on  tho  wing-body  nacollo  poworod  by  tho  two  ojoctors  with  nozzles  in  tho 
forward  or  aft  position,  and  the  wing-body  without  nacelles  poworod  by  tho  lift  Jot. 

A  standard  series  of  runs  was  used  throughout  tho  tost  although  not  all  configura¬ 
tions  woro  tostod  at  all  conditions.  Tho  series,  shown  in  Table  I-V,  consists  of  an 
unpoworod  anglo-of-attack  variation,  followed  by  anglo-of-attack  variation  at  sovoral 
velocity  ratios,  and  velocity  ratio  variation  at  fixod  angles  of  attack.  A  similar  pattern 
was  adoptod  for  runs  involving  sidoslip.  Both  force  and  pressure  data  woro  rocorded 
for  most  runs. 

Tho  range  of  variables  tostod  includod  anglo-of-attack  variations  of  0  to  20 
dogroos,  angles  of  sideslip  from  0  to  *  12  dogroos,  and  volocity  ratios  of  0  to  ,3  for 
lift  Jot  poworod  configurations  and  0  to  .5  for  ojcctor  poworod  configurations.  Very 
llmitod  testing  was  performed  at  combined  angles  of  attack  and  sideslip.  Testing  was 
accomplished  at  dynamic  pressures  of  0  to  71  psf,  resulting  in  Roynolds  numbers  of 
up  to  1.5xl0(}  por  foot, 

Dynumic  pressure  and  thrust  combinations  used  to  achlovo  those  volocity  ratios 
are  llstod  in  Table  I-VI.  Bocauso  of  tunnel  volocity  limits,  higher  volocity  ratios  woro 
obtainod  with  tho  use  of  roducod  thrust.  Noto  that  those  are  nominal  values.  In  gonoral, 
uctual  thrust  levols  wore  slightly  lowor  than  those  shown. 

A  summary  of  configurations  tostod  is  given  in  Table  J-VII. 
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3.  Results 


The  results  of  this  test  are  presented  in  variot  ts  sections  of  this  report  in  support 
of  the  development  and  verification  of  analytical  prediction  techniques.  Complete 
results  will  be  published  in  two  NASA  Technical  Memorandums  at  a  later  date. 

a.  Presentation  of  Results 

The  force  data  recorded  during  the  test  are  presented  in  several  forms: 

(1)  Body  axis  aerodynamic  coefficients . 

(2)  Stability  axis  aerodynamic  coefficients. 

(3)  Forces  and  moments  nondimensionalized  by  thrust. 

(4)  Forces  and  moments  with  direct  thrust  effects  removed,  nondimensionalized 
by  thrust. 

(5)  Forces  and  moments  with  direct  thrust  effects  removed,  nondimensionalized 
by  dynamic  pressure. 

In  the  reduction  of  aerodynamic  coefficients,  forces  were  nondimensionalized  by 
the  product  of  freestream  dynamic  pressure  and  wing  area.  Longitudinal  moments 
were  nondimensionalized  by  the  product  of  freestream  dynamic  pressure,  wing  area, 
and  wing  mean  aerodynamic  chord  length;  lateral  moments  were  nondimensionalized 
by  the  product  of  freestream  dynamic  pressure,  wing  area,  and  wing  span. 

In  the  reduction  of  thrust  coefficients,  forces  were  nondimensionalized  by  the 
total  calibrated  thrust  of  all  nozzles  operating.  Moments  were  nondimensionalized  by 
the  product  of  total  thrust  and  an  effective  diameter.  The  effective  diameter  is  defined 
as  the  diameter  of  a  circle  equivalent  in  area  to  the  sum  of  the  exit  areas  of  the  operat¬ 
ing  nozzles. 

Thrust  removed  coefficients  were  obtained  by  removing  the  forces  and  moments 
of  the  eperating  engine  simulators,  as  determined  during  static  calibration,  from  the 
balance  measured  loads  prior  to  nondimensionalization. 

b.  Selection  of  Unpowered  Baseline  Data 

Analysis  of  the  vectored  thrust  configuration  tests  has  revealed  that  signifi¬ 
cant  differences  exist  in  unpowered  data  taken  with  the  ejector  inlets  open  and  closed. 
These  differences  exceed  those  directly  attributable  to  vectoring  oi  freestream  flow 
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through  the  undriven  ejectors.  When  these  data  are  used  to  obtain  an  interference 
effect  due  to  jet  operation,  two  significantly  different  answers  result,  dependent  on 
which  "unpowered"  coefficient  is  used  in  the  calculations. 


interference 
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total 

direct 

unpowered 

effect 

load 

thrust 

aerodynamics 

The  effect  of  the  open  versus  the  closed  inlet  on  the  power  off  longitudinal 
aerodynamic  coefficients  is  shown  in  Figure  1-11  as  a  function  of  nozzle  deflection 
angle.  Opening  of  the  inlet  produces  a  lift  increment  of  0.216  with  the  forward  small 
nozzles  deflected  90  degrees.  The  portion  of  this  increment  attributable  to  momentum 
change  of  the  free  flow  passing  through  the  ejectors  and  nozzles  may  be  found  using 
Figures  f  -12  and  1-13.  The  unpowered  inlet  weight  flow  is  found  to  be  0.46  lbs/sec/ 
ejector  producing  an  estimated  vertical  force  of  2. 43  lbs.  This  is  equivalent  to  a 
lift  coefficient  of  only  0.018.  Thus,  the  large  part  of  the  inlet  lift  increment  cannot 
be  attributed  to  flow-through  momentum,  but  must  be  due  to  interference  with  the 
external  aerodynamics. 

Powered  lift  data,  nondimensionalized  by  thrust  and  shown  in  Figure  1-14,  do 
not  reflect  the  significant  differences  caused  by  the  open  inlet  in  the  unpowered  data. 
The  difference  in  powered  lift  is  seen  to  remain  approximately  constant  at  various 
velocity  ratios  and,  in  general,  it  is  small  in  comparison  to  the  lift  loss  equivalent 
to  the  aerodynamic  coefficient  increment  obtained  from  power  off  data.  In  contrast, 
the  change  in  inlet  mass  flow  ratio  produced  by  opening  the  inlet  is  greater  for  the 
powered  case  than  for  the  unpowered  case.; 

Further  examination  of  the  inpowered  data,  as  presented  in  Figure  1-15, 
indicates  that  the  longitudinal  coefficients  of  the  inlet  open  configurations  are  functions 
of  nozzle  deflection  angle  while  those  of  the  inlet  dosed  configuration  are  not..  This 
observation  in  combination  with  the  lack  of  difference  in  powered  data  suggests  that 
the  unaccounted  mlet  increment  may  be  due  to  the  exit  interactions  of  the  inlet  mass 
flow  as  opposed  to  the  variation  of  the  inlet  condition  per  se.  This  conclusion  must 
be  considered  tentative,  however,  because  the  inlet  mass  flow  ratio  as  well  as  the 
exhaust  exit  angle  are  varied  by  nozzle  deflection.  An  examination  of  the  pressure 
data  for  this  configuration  will  give  more  insight  to  the  cause  of  the  open  inlet  effect 
on  power  off  data. 

Wing  pressure  differentials  are  shown  for  both  inlet  conditions  in  Figure  1-16. 
Opening  of  the  inlets  causes  a  pressure  change  similar  to  that  produced  by  a  positive 
angle  of  attack.  Examination  of  both  upper  and  lower  surface  pressures  at  WS  10..  55, 


Figure  1-17,  indi  ates  that  the  effect  is  present  on  both  surfaces,  confirming  that 
opening  of  the  ’nlets  induces  a  positive  angle  of  attack  change  of  one  to  two  degrees 
across  the  wing,  and  that  the  effect  is  not  limited  to  one  of  the  surfaces. 

Pressure  distributions  about  the  six  fuselage  stations  are  shown  in  Figure  1-18. 
The  data  are  shown  as  a  locus  of  pressure  normals  about  the  local  fuselage  contour. 
Lines  within  the  section  indicate  positive  pressures;  lines  outside  the  section  indicate 
negative  pressures.  A+  interior  corners  of  the  body  some  data  have  been  eliminated 
for  clarity.  The  effects  of  the  open  inlet  on  fuselage  pressures  are  noted  below. 

At  the  most  forward  station  FG  7. 3,  the  pressure  change  is  similar  to  that  which 
would  be  produced  by  a  positive  angle  of  attack  of  approximately  two  degrees. 
Changes  occur  on  both  upper  and  lower  surfaces. 

In  front  of  the  nacelles  at  FS  11. 8,  pressure  changes  on  the  fuselage  sides  are 
dominated  by  apparent  changes  in  nacelle  blockage. 

Lower  fuselage  and  nacelle  pressures  are  made  more  positive  in  front  of  the 
nozzles  at  stations  11.8  and  19.65.  Upper  surface  pressures  are  not  affected 
at  these  stations. 

Aft  of  the  nozzles  and  nacelles  at  1 3  26. 425,  opening  of  the  inlets  induces  a 
more  negative  pressure  on  the  lower  fuselage  sides;  other  areas  are  not 
strongly  affected.  No  data  are  available  on  the  upper  surface  of  the  fuselage  at 
this  station. 

At  FS  34.  30,  a  more  negative  pressure  is  induced  about  the  section  due  to 
opening  of  inlets.. 

Upper  nacelle  centerline  pressures,  shown  in  Figure  1-19,  are  unaffected 
except  in  the  vicinity  of  the  inlet  lip.:  Lower  surface  pressures  are  made  more 
positive  in  the  vicinity  of  the  exit. 

Lower  fuselage  centerline  data,  Figure  I-2C,  indicate  that  opening  of  the  inlets 
creates  more  positive  pressures  ahead  of  thi  exit  and  more  negative  pressures  aft  of 
the  exit.; 

Data  from  che  tail  flow  angularity  rake  indicate  increased  do"  nwash  when  the 
inlets  are  open,  Figure  1-21. 


In  summary,  opening  of  the  inlets  produces: 

(1)  An  induced  positive  angle  of  at  ack  on  the  wing. 

(2)  An  induced  positive  angle  of  n.iack  at  the  nose. 

(3)  Positive  pressure  increments  on  the  lower  fuselage  and  nacelle 
ahead  of  the  exits. 

(4)  Negative  pressure  increments  aft  of  the  exits. 

(5)  Little  effect  on  upper  fuselage  pressures. 

(6)  Increased  downwash  at  the  tail  station.. 

These  flow  changes  are  consistent  with  those  associated  with  jet  exit  interference 
effects  at  high  velocity  ratios  and  tend  to  confirm  the  initial  observation  that  the 
"inlet"  effect  is  largely  caused  by  the  exit  of  the  inlet  flow. 

Further  confirmation  can  be  gained  by  observation  of  lower  fuselage  centerline 
pressures  for  the  configuration  equipped  with  the  large  aft  nozzles.  It  can  be  seen  in 
Figure  1-22  that  the  largest  effect  of  the  open  inlet  occurs  in  the  vicinity  of  the  exit 
not  the  inlet.  If  the  differential  lower  centerline  pressures  due  to  the  opening  of  the 
inlet  are  plotted  against  distance  from  the  exit  as  in  Figure  1-23,  it  is  seen  that  data 
for  the  forward  small  90  degree  nozzle  and  the  aft  large  90  degree  nozzle  show  a 
surprisingly  stiong  correlation.. 

It  is  concluded  that  the  proper  unpowered  data  for  the  vectored  thrust  configura¬ 
tions  are  that  taken  with  the  inlets  closed.  Data  taken  with  the  inlets  open,  though  the 
model  was  unpowered  in  the  sense  that  no  drive  pressure  was  supplied  to  the  ejectors, 
are  in  fact  power  on  data  at  a  very  high  velocity  ratio,  v/V.  =  2.  9  for  the  small  forward 
90  degree  nozzle  which  has  been  discussed. 


TABLE  I-I.  MODEL  PROPERTIES  AND  DIMENSIONS 


Wing: 


Span,  b . , 

Area  (projected  ),  S . . 

Mean  aerodynamic  chord,  c 
Location  of  25%  c 

F.  S . 

W.L . 

B.  L . 

Taper  ratio,  X . 

Aspect  ratio,  AR . 

Airfoil  section . . . 

Leading  edge  sweep,  AL£  , 
Quarter  chord  sweep,  A 

°/4 

Horizontal  Tail: 


Span,  bR . 

Area  (projected),  SR . 

Mean  aerodynamic  chord,  cR 
Location  of  25%  cR 

F.S . 

W.L.  ...... 


B.L 


Taper  ratio,  A  R . 

Aspect  ratio,  ARR. . . . 

Airfoil  section . 

Leading  edge  sweep,  A 
Quarter  chord  sweep  A 


leh 

CH/4 


Vertical  Tail: 


Span  (centerline),  by . . 

Span  (exposed),  b  . . 

ve 

Area  (exposed),  Sv  . . 

Maan  aerodynamic  chord  (exposed),  cv 


40.25  in 
324.00  in2 
8.35  in 


23.83  in 
3.555  in 
±  8.94  in 
0.50 
5.00 
63A010 
9.76  deg 
5.96  deg 


22.40  in 

2 

110.72  in 
5.14  in 

44.22  in 
6.62  in 
5.00  in 
.50 
4.50 
63A008 

13.40  deg 
8.86  deg 


11.00  in 
9.00  in 

71. 10  in  2 
£.66  :n 


TABLE  I-L  MODEL  PROPERTIES  AND  DIMENSIONS  (CONTINUED) 


Location  25%  Cy 

F.S .  43.93  in 

W.L .  7.25  in 

B.L .  0.00  in 

Taper  ratio  (exposed),  Xy . 26 

Aspect  ratio  (exposed),  ARy  .  4.56 

Airfoil  section  . . 63A008 

Leading  edge  sweep,  ALE  .  42.9  deg 

Quarter  chord  sweep,  A  „  .  34.0  deg 

/a 


Span ,  bp  . . .  0.65  wing  span 

Chord, o_ .  0.25  local  wing 

"  chord 


Deflection ,  <5^ 
Airfoil  Section 


45.  0  deg 


f  FLAP  ORDINATES 

x/c 

yu/c 

y,/c 

O 

.0350 

-.( 0350 

.0125 

.0544 

-. 0194 

0250 

.0650 

-.0150 

.0500 

.0788 

-0094 

.0750 

.0888 

-.0063 

.IOOO 

.0950 

-.0013 

.1500 

.1069 

-.0003 

.2000 

.1/38 

O 

.3000 

.1169 

O 

.4000 

.1138 

O 

.5COO 

.1050 

O 

.6000 

.0913 

O 

.7000 

.0738 

O 

.8000 

.0525 

O 

.9000 

.0281 

O 

.9500 

.0/50 

0 

1.0000 

.0013 

O 

Overlap 

Gap 

r  •••  >««•■•••••• 

Location  of  momt  -  •  •:  nter : 

F.S. 
W.L 
B.  L. 


0.01  c 
0.02  c 

23.82  in 
0.00  in 
0.00  in 
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Nozzle  Exit  Location* 


TABLE  I-m.  LOCATION  OF  PRESSURE  PORTS  ON  LOWER  FUSELAGE 


TABLE  MV.  WING  PRESSURE  POSTS 


Chord 

% 

Port  Numbers 

Surface 

Sjpanwise  Location  U/(b/2J 

0.25* 

.30 

.52 

.80 

0.0 

— 

1 

24 

47 

70 

U 

2 

25 

1.0 

L 

13 

36 

U 

48 

1.5 

L 

59 

2.5 

U 

3 

26 

49 

71 

u 

14 

37 

60 

81 

5.0 

u 

4 

27 

50 

72 

L 

15 

38 

61 

82 

10.0 

u 

5 

28 

51 

73 

L 

16 

39 

62 

83 

15.0 

u 

6 

29 

52 

74 

84 

L 

17 

40 

63 

25.0 

U 

7 

30 

53 

75 

L 

18 

41 

64 

85 

U 

8 

31 

54 

76 

40.0 

L 

19 

42 

65 

86 

U 

9 

32 

55 

77 

55.0 

L 

20 

43 

66 

87 

U 

10 

33 

56 

78 

70.0 

L 

21 

44 

67 

88 

U 

11 

34 

57 

79 

85.0 

L 

22 

45 

68 

8P 

U 

12 

35 

58 

80 

95.0 

23 

46 

69 

90 

Note:  U  —  Upper  Surface 
L  —  Lower  Surface 
*  —  Nacelle  Centerline 
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Velocity 

Ratio 

Unpowered 

0.3 

0.2 

0.1 

0  to  0.3/0.5 

0  to  0.3/0.5 

Unpowered 

0.3 

0.2 

0.1 

0  to  0.3/0.5 

0  to  0.3/0.5 


7.  RUN  SEQUENCE 


Angle  of 

Angle  of 

Attack 

Sideslip 

Degrees 

Degrees 

0  to  20 

0 

0  to  20 

0 

0  to  20 

0 

0  to  20 

0 

0 

0 

10 

0 

0 

-12  to  12 

0 

-12  to  12 

0 

-12  to  12 

0 

-J2  to  12 

0 

-8 

0 

8 

[*TT 


TABLE  l-VL  TEST  OOMZTOTXS 


EJECTORS 
3.67-tnch  D  Kozzlet 


EJECTORS 
4.S-incfc  D  Nozzles 


LIFT  JET 
2.2S-inch  D  Nozrle 


HIGH  THRUST  LOW  THRUST 


LIFT  JET  WITH  EJECTORS,  3. 67-inch  D  Nozzles 


TABLE  I-Vn.  TEST  PROGRAM  SUMMARY 
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BODY*  WING- TAIL 
&  FLAP 


FIGURE 


FIGURE  1-2.  WING -BODY- NACELLE  WITH  TAIL  FLOW  ANGULARITY 
RAKE,  POWERED  BY  SMALL  FORWARD  NOZZLES 
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FIGURE  1-3.  WING-BODY- NACELLE  TAIL  WITH  FLAP,  POWERED  BY 

AFT  SMALL  NOZZLES 


FIGURE  1-4.  WING-BODY-TAIL,  POWERED  BY  LI7T  JET 


sec  w  lb/ sec 


96 


lb/ sec  w.  lb/sec 


Lift  Jot 
T'=  72.4  lbs 


Small  Ejector 
T'=  77.1  lbs 


FIGURE  1-6.  TYPICAL  EXIT  PROFILES  ALONG  LONGITUDINAL 

NOZZLE  CENTERLINE 


FS  26.425 


FS  34.30 


FIGURE  1-7.  CIRCUMFERENTIAL  FUSELAGE  PRESSURE  PORTS 


143 


142 


FS  18.55 


FIGURE  1-7  (concluded) 


201 


FIGURE  1-8.  NACELLE  CENTERLINE  PRESSURE  PORTS 
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12  14  16 

h/D 

FFECT  OF  MODEL  HEIGHT,  a  -6° 
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Nozzle  Deflection  ~ Degrees 


FIGURE  1-11,  EFFECT  OF  FREE  FLOW  THROUGH  INLET  ON 
MODEL  AERODYNAMIC  LOADS  AT  a  -  0  DEG 
WING-BODY-NACELLE  WIT  !  SMALL  NOZZLES 


FIGURE  1-12.  INLET  WEIGHT  FLOW  FOR  UNPOWERED  EJECTOR 

SMALL  NOZZLES 
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Thrust  Coofflciont~C  T*>tal  Nor  isle  Thruit* 


0  20  40  60  80  100 

Nozzle  Deflection  ~6.  ~  Degrees 


FIGURE  1-15..  UNPOWERED  AERODYNAMIC  COEFFICIENTS 
WING-BODY-NACELLE,  SMALL  NOZZLES  a  =  0 


E 


»9 


FIGURE  1-18.  EFFECT  OF  FREE  FLOW  THROUGH  INLET  ON  FUSELAGE  PRESSURES 
WING-BODY-NACELLE,  FORWARD  SMALL  NOZZLES,  6.  =  90° 


FIGURE  1-18  (concluded) 


I - 1 - 1 - 1 _ L _ I _ L _ i  ■ 

14  15  16  17  18  19  20  21  22 

Fuselage  Station 


FIGURE  1-19.  EFFECT  OF  FREE  FLOW  THROUGH  INLET  ON  NACELLE 

CENTERLINE  PRESSURES 


5 


Wing  Body  Nacelle,  Forward  Small  Nozzles, 


FIGURE  1-20.  FREE  FLOW  INLET  EFFECT  ON  LOWER  FUSELAGE  CENTERLINE  PRESSURES 


Wing  Body  Nacelle,  Forward  Small  Nozzles,  6j  =  90* 
Unpowered 


Wing  Height 

_ I _ 


10 


12 


14 


W.  L. 


FIGURE  1-21., 


EFFECT  OF  FREE  FLOW  THROUGH  INLET  ON  DOWNWASH 
AT  THE  TAIL  FLOW  ANGULARITY  RAKE 
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Nacelle,  Aft  Large  Nozzles 
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FIGURE  1-22.  FREE  FLOW  INLET  EFFECT  ON  LOWER  FUSELAGE  CENTERLINE  PRESSURES 
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FIGURE  1-23.  EFFECT  OF  OPEN  INLET  ON  LOWER  FUSELAGE  CENTERLINE 
PRESSURES  WITH  FORWARD  AND  AFT  NOZZLES 


APPENDIX  n 


BALANCE  CORRECTIONS  AND  CALIBRATION 


Calibration  of  the  basic  balance  was  performed  by  NASA-  Langley  using  well 
established  procedures.  The  effects  of  tLe  air  supply  line,  which  spanned  the  balance, 
were  determined  fay  loading  of  the  balance  both  with  and  without  the  airline  using  the 
em?  3eries  of  loads.  A  matrix  of  correction  coefficients  was  then  derived  which 
caused  the  corrected  airline-on  loads  to  match  the  measured  airi ice-off  or  applied 
loads.  The  loads  were  applied  through  the  use  of  a  complex  device  utilizing  a  system 
of  weight  pans,  fnlcrmns,  and  levers.  The  resulting  correction  factors  are  shown 
below.  The  6x6  matrix  indicated  in  Appendix  I  is  shown  here  as  two  3x3  matrices 
as  longitudinal  and  lateral  components  were  not  found  to  interact. 


{Corrected 
Loads 


1  Correction  ( 

(Measured  Loads 

1  Matrix  j 

j  Airline-On 

. 

N 

1.0 

0 

-.0021 

N 

A 

= 

-.0019 

1.002 

-.0014 

A 

M 

.056 

-.0115 

A  .107 

M 

y 

C1 

1 

y 

M 

1.025 

-.0181 

0 

M 

X 

X 

M 

= 

-.0592 

1.085 

0 

M 

z 

z 

Y 

C 

.0008 

-.0014 

.9980 

Y 

Bal 


Bal 


Because  of  initial  difficulties  in  obtaining  the  longitudinal  corrections  shown 
above,  a  second  loading  of  normal  force  and  pitching  moment  was  made  with  the  use 
of  a  bar  and  moving  weight  pan.  This  resulted  in  a  different  set  of  longitudinal 
corrections.  Interactions  due  to  axial  forte  were  taken  from  the  previous  loading. 


N 

.9945 

0 

-.00198 

A 

.01167 

1.002 

-.00316 

M 

tr 

.27 

0 

1.109 

N 

A 

M 

y 
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o 


Bal 


This  matrix  dfiSers  sigstficwtty  from  the  original  in  the  eaUbratioe  of  corrected 
axial  force  and  pitching  moment.  The  pitching  moment  difiereace  wndd  cause  normal 
forces  to  ajpear  approximately  .2  inch  farther  forward  if  foe  second  matr  ix  is  need  in 
preference  to  foe  first.  The  ase  of  foe  second  matr  ix  woald  also  indicate  increased 
drag  with  positive  normal  force  and  decreased  drag  with  positive  pitching  moment  when 
compared  with  data  reduced  by  foe  first  matrix. 

The  second  matrix  has  been  used  to  reduce  foe  data  presented  in  this  report 
because  more  experience  was  available  with  the  method  of  loading.  In  addition,  the 
second  matrix  gate  better  results  on  the  position  of  a  small  test  weight  and  it 
indicated  that  file  lift  jet  center  of  pres  sore,  which  was  not  constant,  approached  the 
geometric  center  of  foe  exit  at  higher  thrusts  as  opposed  to  moving  away.  However, 


the  results  of  these  checks  were  not  conclusive. 

The  matrix  below  Indicates  the  final 

*0 

matrix  oi  correction  coefficients  K. 

.9945 

0 

-.00198 

0 

0 

0 

.01167 

1.002 

-.00316 

0 

0 

0 

K  = 

.27 

C 

1.109 

0 

0 

0 

0 

0 

0 

1.025 

-.0181 

0 

0 

0 

0 

-.0592 

1.085 

0 

0 

0 

0 

.0008 

-.0014 

.9960 

The  matrix  shown  below  will  convert  the  longitudinal  data  presented  in  this  report 
to  the  form  it  would  have  had  if  the  original  corrections  wer  used. 


N 

1.0055  0  -.0001 

N 

A 

= 

-.0142  1.000  .0016 

A 

M 

-.213  -.0115  .9975 

M 

y 

r. 

V 

Nominal  accuracy  of  the  basic  balance  is  indicated  below. 


Component 

N 

A 

Y 

Mx 

D1 

Mz 

Accuracy  —  lb  or  in. -lb 

±2.5 

±1.0 

±1.5 

±5 

±15 

±10 
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appenux  ra 


EFFECT  OF  FORWARD  SPEED  ON 
EJECTOR  JET  ENGINE  SIMULATORS 

fa  order  to  investigate  the  effect  of  forward  speed  on  foe  thrust  of  foe  ejector  jet 
engine  simulators,  foe  beUmooth  used  daring  static  calibration  was  attached  to  the 
ejector  unit  mounted  in  foe  test  section.  The  inlet  mass  flow  variation  with  forward 
speed  was  then  determined.  The  effect  of  forward  speed  on  inlet  mass  flow  is  shown 
in  Figure  ID-1.  The  changes  shown  are  small,  less  than  six  percent,  and  may  reflect 
changes  in  the  effective  area  of  the  bellmouth  which  was  calibrated  under  static  condi¬ 
tions. 

If  foe  thrust  of  foe  ejectors  is  assumed  to  vary  with  foe  square  of  the  total  mass 
flow,  the  change  in  thrust  may  be  computed. 


(pAjUj)  U, 


The  computed  change  in  thrust  du~  to  forward  speed  is  shown  below.  The  maxi¬ 
mum  thrust  change  occurs  at  the  highest  forward  speed. 


Nozzle 

T0.U> 

T 

Max 

o 

Q,  psf 

136 

1.026 

Small  Forward 

105 

1.026 

■89 

6i  =  90° 

30 

1.059 

0-70 

60 

1.061 

0-60 

Large  Aft 

_ 

178 

1.001 

0-60 

A  greater  change  occurs  at  lower  thrust  levels  because  of  the  larger  ratio  of 
secondary  to  primary  flow. 
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The  maximum  thrust  errors  are  of  a  magnitude  less  than  twice  the  nominal 
accuracy  of  the  normal  force  beam,  and  foe  maximum  pitching  moment  errors  are  of 
foe  same  relative  magnitude.  Correction  for  foe  errors  would  be  significant  only  at 
foe  higher  velocity  ratios  tested.  At  maximum  dynamic  pressures  correction  would 
result  in  a  three  percent  reduction  in  velocity  ratio  and  a  six  percent  redaction  in 
thrust  nondtmensionalized  data.  As  foe  majority  of  data  show  a  positive  slope  at  high 
velocity  ratios,  the  effect  is  minimized,  because  the  corrected  datum  tends  to  move 
along  the  established  curve. 
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APPENDIX  IV 


NORMAL  FORCE  AND  PITCHING  MOMENT  IN  THE 
LIFT  JET  WAKE 


The  jet  model  includes  no  mechanism  which  will  account  for  the  wake  region 
behind  the  jet.  Consequently,  calculations  utilizing  the  jet  flow  field  program  are  not 
in  good  agreement  with  test  data  in  this  region.  As  indicated  in  Section  HI  the  wake 
region  contributes  the  major  part  of  the  loads  on  the  fuselage  so  that  it  is  desirable  to 
have  some  method  for  estimating  the  integrated  force  and  pitching  moment  for  this 
region. 

To  enable  estimates  of  the  fuselage  loads  in  the  wake  behind  the  lift  jet  to  be 
made,  it  has  been  assumed  that  the  difference  between  test  data  and  calculations  in 
this  region  is  the  same  as  the  difference  between  test  data  and  calculations  for  the 
component  model  of  Reference  99. 

Thus,  for  example,  if  the  fuselage  is  at  zero  degrees  angle  of  attack  and  sideslip, 
the  wake  contribution  to  the  interference  lift,  AL.,  is  given  by 


in  which  T  is  the  jet  thrust.  This  double  integral  has  been  evaluated  for  a  range  of 
position  along  the  test  fuselage  and  is  shown  in  Figure  IV- 1. 

The  wake  contributions  to  pitching  moment  may  be  determined  in  a  similar  man¬ 
ner.  The  results,  for  the  three  velocity  ratios  .125,  .2,  .3  are  shown  in  Figure  IV-2. 
The  moment  axis  has  been  taken  through  the  center  of  the  jet. 

In  Section  IH  adjustments  to  the  calculations  have  been  made  to  account  for  the 
jet  wake  effect.  The  jet  wake  contributions  were  obtained  from  Figure  IV-1  for  a 
wake  length  of  13. 5  jet  diameters. 


225 


FIGURE  IV-1.  JET  WAKE  CONTRIBUTIONS  TO  FUSELAGE  LIFT 


IGURE  IV-2.  JET  WAKE  CONTRIBUTIONS  TO  FUSELAGE  PITCHING  MOMENT 
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